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Neue kritische Drehzahlen von einfach besetzten Wellen 
Von W. Keller 


1. Einleitung. Wie H. Melan in einer Arbeit iiber kritische Drehzahlen von Maschinenwellen 
mit Langsbelastung! gezeigt hat, werden die Biegesteifigkeit und damit auch die kritischen Dreh- 
zahlen von Wellen durch eine Zugbelastung hinauf- und durch eine Druckbelastung herabgesetzt. 
Kine ahnliche Erscheinung fand R. Grammel? bei Wellen, die ein Drehmoment zu iibertragen haben. 
Ein solches verkleinert namlich unabhangig davon, welchen Drehsinn es besitzt, stets die Biege- 
steifigkeit der Welle und daher auch ihre kritischen Drehzahlen. Beiden Untersuchungen liegt eine 
zeitliche konstante Belastung zu Grunde. 

Indessen kommt es im Maschinenbau auch haufig vor, da auf eine Welle neben konstanten 
noch pulsierende Torsionsmomente und pulsierende Axialkrafte wirken, wodurch sich ihre Biege- 
steifigkeit mit der Zeit periodisch andert. Da diese in den Bewegungsgleichungen einer auf der 
Welle sitzenden Scheibe als Koeffizient auftritt, wird in diesem Fall das Verhalten der Scheibe 
durch Differentialgleichungen mit periodischen K oeffizienten bestimmt. Solche Gleichungen besitzen 
unter gewissen Umstanden instabile Lésungen, die also mechanisch gesehen ein unbegrenztes Aus- 
weichen der Scheibe beschreiben. Hangt die Frequenz der pulsierenden Belastung von der Drehzahl 
ab, so sind kritische Drehzahlen zu erwarten, die sich von den seither bekannten zumindest durch 
eine véllig andere Ursache unterscheiden. 

Wie die durch diese Uberlegungen angedeuteten Zusammenhange im einzelnen sind, wird in der 
vorliegenden Arbeit untersucht. Dabei zeigt sich, daB es bei einem konstanten oder auch harmonisch 
schwankenden Antriebsmoment neben der kritischen Drehzahl erster Art unendlich viele Dubletts 
von kritischen Drehzahlen zweiter Art gibt, wenn die 
Axialkraft oder das Torsionsmoment harmonisch schwankt. 
AuBerdem gibt es noch kritische Drehzahlbereiche, deren 
Breite mit dem pulsierenden Anteil der Axialkraft oder 
des Torsionsmomentes zunimmt, und deren Anzahl von 
der GréBe der stabilisierend wirkenden Dampfungskrafte 
abhangt. 


2. Die Bewegungsgleichungen. Gegeben sei eine irgend- 
wie gelagerte Welle von iiberall kreisférmigem, sonst je- 
doch beliebig veranderlichem Querschnitt, die um ihre im 
unverformten Zustand gerade Achse rotieren médge. Sie 
trage eine einzige Scheibe, deren Masse man sich in ihrem 
Schwerpunkt (Massenmittelpunkt) konzentriert denkt. Von 
dem sekundaren Einflu8 einer etwaigen Kreiselwirkung 
wird also abgesehen. AuBerdem soll vorausgesetzt werden, 
daB die Wellenmasse gegeniiber der Scheibenmasse ver- pp. 1. Bezeichnungen an der rotierenden Scheibe. 
nachlassigt werden kann. 

In einem raumfesten Bezugssystem, dessen Ursprung O auf der Verbindungsgeraden g der Lager- 
mitten liegt, und dessen (y, z)-Ebene auf der Geraden g senkrecht steht und aufserdem noch den 
Schwerpunkt S enthalt, habe dieser die Koordinaten y und z (Abb. 1). Die Lage des Durchstof- 
punktes D der gebogenen Wellenachse mit dieser Scheibenmittelebene (y,z), der wegen der stets 
vorhandenen Zentrierungsfehler der Scheiben einen Abstand e vom Schwerpunkt S besitzt, ist dann 
mittels des Winkels py, den die Gerade DS mit der positiven y-Achse bildet, durch die Koordinaten 
y —ecosp, z—e sing, bestimmt. 

Die wegen der Elastizitat der Welle auf die Scheibe wirkende Kraft & = (0, K,, K,) kann in 
den Koordinaten des Schwerpunkts ausgedriickt werden. Mit der EinfluBzahl x, d.h. der Aus- 


1 H. Melan, Z. ost. Ing.- u. Arch.-Ver. 69 (1917) S. 610. 
2 R. Grammel, Stodola-Festschrift S. 180, Ziirich, Leipzig 1929. 
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lenkung des Punktes D infolge einer in D senkrecht auf g stehenden Einheitskraft von gleicher Richtung 
wie die Auslenkung!, werden die Komponenten der Kraft ®% in —y- und —z-Richtung 


K, =—(y—e cos), (1) 

K, =—(—esing), (2) 
und somit lauten mit der Scheibenmasse m die Bewegungsgleichungen des Schwerpunkts? 

my ey =— cosy, 

mi+—2=<sing o 


Dabei sollen alle Einfliisse auf die Wellenbiegung durch die EinfluBzahl « beriicksichtigt werden. 

Da bei Wellen, die gedriickt und tordiert werden, die EinfluBzahlen stark davon abhangen, wie 
das Torsionsmoment in die Welle eingeleitet wird, soll zunachst festgelegt werden, welcher Art die 
an der Welle angreifenden Krafte und Momente sein sollen. Von der Zug- bzw. Druckkraft 4b wird 
angenommen, da sie axial ist, d. h. daB ihre Wirkungslinie stets mit der Verbindungsgeraden (g) 
der Lagermitten zusammenfallt. In einem (x, y, z)-System, dessen x-Achse in der urspriinglich 
geraden Wellenachse liegt, besitzt also die Kraft S$ vom Betrag P die Komponenten 


B = (P, 0, 0). (4) 
Sie sei positiv, wenn sie auf die Welle driickt und negativ, wenn sie an ihr zieht. Das Torsions- 
moment ¥$ vom Betrag W wird dagegen als semitangential vorausgesetzt, d. h. der Momentvektor YY 
halbiere immer den Winkel zwischen der unverbogenen Wellenachse und der Tangente an die 
verbogene Wellenachse. Mit den Auslenkungen v(x) und w(x) der verbogenen Wellenachse in y- 
bzw. z-Richtung ergibt sich somit die Komponentendarstellung des Torsionsmomentes in erster 
Naherung zu 


1 
W =(W, 5 u,W, wi Ww). (5) 


wobei Striche Ableitungen nach x bedeuten und der Index 1 die Stelle x, kennzeichnet, an der das 
Moment Y angreift. Ein solches Torsionsmoment wird z. B. von mehr als zwei gleich groBen und 
iiber den Scheibenumfang gleichmaBig verteilten Tangentialkraften erzeugt’, wie sie etwa bei mehr- 
fach oder gleichmaBig beaufschlagten Turbinenscheiben auftreten. Belastungen von der Art (4) 
und (5) fihren bekanntlich zu Verschiebungsvektoren«, die dieselbe Richtung haben wie die Einheits- 
krafte, die sie erzeugen.* 

In alteren Untersuchungen wurden oft axiale Momentvektoren YS = (W, 0,0) beniitzt, die jedoch 
nur selten vorkommen und im allgemeinen keine konservative Belastung darstellen®. Wenn daher 
im folgenden von einem Torsionsmoment die Rede ist, so soll es stets semitangential sein. 

| Betrachtet man nun eine stationdre Drehung der Scheibe mit der Winkelgeschwindigkeit P=O—= 
konst., so fiihren bei festem « die Gleichungen (3) in bekannter Weise zu dem kritischen Zustand 


erster Art®. Nimmt man an, daf sich infolge eines periodisch schwankenden Antriebsmomentes der 
Drehwinkel der Scheibe gemaB 


yg =ot + Ht) (6) 
andert, wobei w die konstante mittlere Drehgeschwindigkeit und ((¢) eine periodische Funktion der 


Zeit ist, so treten neben den kritischen Drehzahlen erster Art auch solche zweiter Art auf7. Beide 
sind Resonanzerscheinungen mit ihren bekannten Merkmalen. 


Bei einer pulsierenden Beanspruchung der Welle auf Druck und Torsion wird jedoch auch die in 
den Bewegungsgleichungen auftretende EinfluBzahl « (P, W) eine periodische Funktion der Zeit t 


' Es sind auch bei rotationssymmetrischen Wellen Belastungen méglich, bei denen die Verschiebung «@ nicht 
dieselbe Richtung hat wie die Einheitskraft [siehe Ziegler, Ing.-Arch. 20 (1952) S. 384], jedoch werden solche 
Falle im weiteren nicht behandelt. 

? Siehe C. B. Biezeno und R. Grammel, Technische Dynamik, Bd. 2, S. 156, 2. Aufl., Berlin-Géttingen-Hei- 
delberg 1953. 

3 Siehe H, Ziegler, Ing.-Arch. 20 (1952) S. 377. 

SEibdasusoa. 

° H. Ziegler, Z. angew. Math. Phys. 3 (1952) S. 96. 

6° C. B. Biezeno und R. Grammel, a. a. O. S. 155. 

’ C. B. Biezeno und R. Grammel, a. a. O. S. 180. 
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sein. In diesem Fall hat man mit (3) zwei voneinander unabhangige Hillsche Differentialgleichungen. 
Um sie explizit anschreiben zu kénnen, muf zunichst « als Funktion der periodisch schwankenden 
Belastungen P und W bestimmt werden. Damit diese Funktion « (P, W) eine iibersichtliche Form 
annimmt, ist es zweckmafig, sie in eine Taylorsche Reihe zu entwickeln, was dadurch gerechtfertigt 
ist, da} die Amplituden der Belastungen P und W praktisch wohl immer weit unter den entspre- 
chenden Knickwerten einer statischen Axialkraft und eines statischen Torsionsmomentes liegen. 
Dazu empfiehlt es sich, die dimensionslosen Abkiirzungen 


TP Ue W 1\2 
Cag “= (F5| (7) 
einzuftthren, wobei J die Lange der Welle, E der Elastizitatsmodul und J das axiale Flachentrag- 
heitsmoment des Wellenquerschnitts ist. Beachtet man, da® fiir die Wellenbiegung und damit 
auch fiir die EinfluBzahl (0, ~) zwar das Vorzeichen von P, nicht aber das von W wesentlich ist, 
so stellt man leicht fest, daB «(0, ) eine gerade Funktion von g und y sein muB und demnach in 
erster Naherung die Form! 


Ba(o, 0 aa(0, 
a(Q, ) = (0, 0) + ¢ oe + s SNe Es (8) 
e e=0 M u=0 
besitzt. Mit den stets positiven GréBen 
= 1 da(o, 0) oe 1 0a(0, ) = 
«(0, 0) =a » Go ao? ee > % one tas Xe (9) 
kann man daher in derselben Naherung auch 

1 ul 

pee ea ty 0? — Hy pl? — >>) (10) 


setzen. 
Nun sollen sich die Axialkraft P und das Torsionsmoment W je aus einem konstanten und einem 
mit der Kreisfrequenz Q pulsierenden Anteil zusammensetzen und zwar so, daB 


Pi Pecos.) t, W=W,+ W,cos Qt (11) 
ist. Dabei kénnen die GréBen P,, W, und W, als positiv vorausgesetzt werden, ohne da dadurch 
die Allgemeinheit des Ansatzes (11) eingeschrankt wird. Dann bildet man analog zu (7) die Abkiir- 


zungen 
PoP 2__ PP 2 (FE 7) pbs G 5) (12) 


O=E’ Oy (Hie Ho = Tag) es a | 
Mit ihnen und den Abkiirzungen (7) findet man aus (11) leicht 
o = 02 + e? cos Qt, 
ee +o + 2 My fy cos Qt + 5 Hi cos 2 OT. 


womit die Entwicklung (10) in 
i) 1 

3 == 1 Hy 03 — Hy U3 — > Ha Mi — (41 OF + 2 Hp fly fy) C08 Qt— -—F x_ujcos2Mt—---| (13) 

es 0 

iibergeht. Daraus ist ersichtlich, daB8 nun die Differentialgleichungen (3) periodische Koeffizienten 


besitzen. 
Mit den dimensionslosen Veranderlichen 


t= Ot, © n=>, C=> (14) 


und den dimensionslosen Festwerten 


2 © © 
= 1m e— mm — za) 
2 
1 = — Fa ei + 2 2 Mot)» (15) 
we 1 
Y2=— Ge 7 el 


1 C. B. Biezeno und R.Grammel, a.a.O. S.179. 
6* 
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worin 
J 


Veo m 


die kritische Drehgeschwindigkeit der unbelasteten, querschwingenden Welle ist, lassen sich die 
Bewegungsgleichungen (3) in der Form 


(16) 


Mo — 


a (A+, cost + y,cos2t) 7 =r + y_ cost + yz, cos 27) cos, 
eS (17) 


aa t (A+, cost + y2 cos 27) $= (1+; cost + 72 cos 21) sin @ 


schreiben. 

Da bei praktischen Problemen meist W, << W, und demnach yj > pj ist, kann man in (15) 
die Glieder mit yu? vernachlassigen. In den Differentialgleichungen (17) verschwinden dann ins- 
besondere die Glieder mit cos 2 T, so daB das Verhalten der Scheibe durch zwei unhomogene Mathieu- 
sche Differentialgleichungen beschrieben wird. Um diesen eine fiir die weitere Untersuchung giin- 
stigere Gestalt zu geben, fiihrt man mit 


of = 0} (L—%4 08 — m8) (18) 


das Quadrat der kritischen Drehgeschwindigkeit der mit P) und W, vorbelasteten We'le’ ein und 
kommt damit zu den Bewegungsgleichungen 


22, 07 § 
t+ |L— et + 2 me Mo ey) = cost 
0) 


F cos P 
2 2 : 
Wr OT a is 


e w§ 
7 =T[— ei + 2 Ht Ug [My) = COST 
Of 


(19) 


$ we 
eas L (% 01 + 2 Ha Ho Ha) og COST 


2 
C= [LG 08 +2 Ha Ho Ma) 53 0087] sin > 
kk 


Dabei wird die Frequenz 2 der auBeren Belastung sehr haufig mit der Drehgeschwindigkeit w der 
Scheibe gekoppelt sein, etwa so, daB auf jeden Umlauf der Welle vy Impulse fallen, also 


Q=rve (20) 
ist. Setzt man in (19) zur Abkiirzung 
oes 2 2 OF e 
open ane Mee Bata 2 Mao Pa cece SG a se (21) 


so erhalt man bei Beriicksichtigung von (6) und (20) die endgiiltige Fassung der Bewegungsglei- 
chungen: 


@ bal + (1—ecost) 7 =y (1 — € cosT) cos E + Hz)], (22a) | 
are 
he a + (1—ecost) € =y(1—€ cost) sin eS + H(t)|. (22b) 


Dabei sind die Parameter q?, ¢ und y stets positiv. 
Fiir die Berechnung der bis jetzt noch nicht naher angegebenen Funktion #(t) wird nun ein 


Antriebsmoment zugrunde gelegt, das mit derselben Frequenz 2 = yw pulsiert wie die auBere 


Belastung P und W. Dann gilt in erster Naherung? 


At) = — a r5 008 (t + B,) , (23) 


worin @ das Tragheitsmoment der Scheibe beziiglich einer auf ihrer Mittelebene senkrecht stehenden 
Schwerpunktachse, a, die Amplitude des Antriebsmomentes und f, ein konstanter Phasenwinkel ist. | 


Die Gleichungen (22) lassen sich, ohne bestimmte Wellenmafe festzulegen, weiterbehandeln, da 


in ihnen keine benannten GréBen mehr auftreten. Zuniichst soll jedoch im nachsten Abschnitt an | 


einem Beispiel gezeigt werden, wie man die Koeffizienten x,, x, und damit auch wf und ¢ be- 
rechnen kann, wobei wir uns auf eine Welle von iiberall gleichem Querschnitt beschranken. 


3. Die Berechnung der Einflufzahlen. Um die Einflu8zahl a einer Welle explizit zu berechnen, 
geht man von ihren Biegegleichungen aus. Diese lauten? im Rahmen einer Theorie erster Ordnung 
— man darf sich auf kleine Auslenkungen v(x) und w(x) beschranken — fiir eine mit $$ und YW be- 


1 C, B. Biezeno und R. Grammel, a. a. O. S. 179. 
2 C. B. Biezeno und R. Grammel, a. a. O. S. 181. 
* Siehe H. Ziegler, Z. angew. Math. Phys. 2 (1951) S. 281. 
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lastete Welle von iiberall gleicher Biegesteifigkeit EF J 


ae (24) 
EJwY—Wo'" 4+ Pw" =0. 
Es empfiehlt sich, diese beiden Gleichungen mit der komplexen Veranderlichen 
r(x) = v(x) + i w(x) (29) 
in die eine Differentialgleichung 
EJrY—iwWr'' + Pr'-=0 (26) 


zusammenfassen. Deren allgemeinste Lésung kann man zwar leicht angeben, die Bestimmung der 
darin auftretenden Integrationskonstanten fiihrt jedoch schon bei einfachen Randbedingungen zu 
umstandlichen Rechnungen, nach denen die Aus- 
driicke fiir «(P,W) meist sehr undurchsichtig a 
werden. Hs ist daher zweckmaBig, die Gleichung 
(26) mit einem Naherungsverfahren zu lésen. 
Betrachtet man beispielsweise eine vor der 
Biegung kreiszylindrische, beiderseits frei drehbar 
gelagerte Welle (Abb. 2), die an der Stelle x = a 
(der Ursprung des beniitzten Koordinatensystems 
soll dabei im linken Lager der Welle liegen) mit 
einer in y-Richtung wirkenden Kraft Q belastet ist 
und nur in ihrem linken Teil (0 < x <a) ein 
semitangentiales Moment Y, jedoch keine Axial- 
kraft ‘8 iibertragt, so erkennt man sofort, daB 
in diesem Fall die elastische Linie des verdrillten jo) 
Wellenteils eine Raumkurve ist, deren Differen- 


tialgleichung zufolge (26) die Form 
rt Wr! = 0 (04s a) (27) 
besitzt. Fiir den restlichen Teil der Biegelinie gilt dagegen 


Abb. 2. Auf- und GrundriB der durch Q und W zu einer 
aumkurve verbogenen Wellenachse. 


TOP g fire 28) (Gx =), (28) 

Nach einer ersten Integration erhalt man 
EJr''—iWr'=—(A,+iA,) fir O0S4Sa, (29) 
E jr” = B, +1 B, Dar uae « (30) 


Die reellen Integrationskonstanten 4,, A, und B,, B, stellen dabei die Reaktionskrafte des linken 
bzw. des rechten Lagers dar. Man kann sie in den GréBen W und Q ausdriicken, indem man die aus 


Abb. 2 leicht abzulesenden! Momentengleichgewichtsbedingungen der Welle (mit b = / — a) 


bQ , Ww'(0) W w’(a) fie W v‘(a) W v’(0) 
ee ay et 2 21 Mr 2 (31) 
aQ , Ww'(a) W w’(0) Ree Wv'(0)  Wrv'(a) | 
eg rca 8 oy oa ae oy 21 
anschreibt. Fiir die weitere Rechnung ist es vorteilhaft, die dimensionslosen GréBen 
z 2 ee 32 
f= 5 on ae o3 =F ( ) 
und die EinfluBzahl 
3 
ie aye. (33) 


der nur mit Q belasteten Welle neben der bereits verwendeten Abkiirzung = W IE J einzufibren 
und die Gleichungen (29) und (30) mit Riicksicht auf (31) in 


mO—tpr'G)—rO) fr 05656, (34) 


tr : a 
r Sa CISL Pe Bae 
1 32 


Be Die in der Abb. 2 eingetragenen Pfeile weisen jeweils in der Richtung, die fiir die betreffende Kraft- bzw. 
Momentkomponente als positiv gewahlt wurde [siehe Gleichung (5)]. 
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und 
"= Pam Q—FHIG)—r'O)] fir &SES1 (35) 


umzuformen, wobei jetzt die Ableitungen nach € mit Strichen bezeichnet werden. 
Die Randbedingungen verlangen, daB fiir € = 0 die Auslenkungen v = w = 0 und die Kom- 
ponenten des Biegemomentes 


EJ i» W 


2 ar, Ui eee tie) WW, 
2 Srrig tcc Tete 1 


9 Bp w = 


(siehe Abb. 2) sind, ferner daB fiir € = 1 sowohl v = w = 0 als auch v" = w” = 0 ist. AuBerdem 
mu an der Stelle € = &, die Biegelinie und ihre erste Ableitung stetig sein. Diese Bedingungen, die 
man einfacher in der Form 


> (36) 
(Gy rte 0) (Ey 0) re at (Gaal) eee lode) 

schreibt, reichen aus, um die bei der Lésung der Gleichungen (34) und (35) auftretenden Integra- 

tionskonstanten zu bestimmen. 


Die Voraussetzung u < 1 legt nahe, die Lésung r(€) in eine Potenzreihe nach y zu entwickeln. 
Man setzt daher 


r(é) = Sn r,(8 (37) 


und fihrt, um das unbekannte Funktionensystem r,(€) zu ermitteln, die Summe (37) in die Diffe- 
rentialgleichungen (34) und (35) sowie in die Randbedingungen (36) ein. Dann verlangt man, daB 
die Summen der Glieder mit gleichen Potenzen von 4 verschwinden und findet so fir 0 < é < &, 
die Differentialgleichungen 


“tt 3 
ri tam =0, (38,0) 
$1 $2 
et ene ape Hee (EO enue = (bes 54 (38,n) 
sowie die Randbedingungen 
ro(0) —— 0 9 ro (0) == 0 > (39,0) 
rm(0) = 0, Fn (0) re 10) = 0 Gn, (39,m) 
Entsprechend gewinnt man fiir &; < € <1 die Differentialgleichungen 
ts 3 
Li = Ee Xo Q = 0, (40,0) 
im +> [Pm —1(6:) — Tm —1(0)] = 0 (el 2 ee) (40,m) 
und die Rand- bzw. AnschluBbedingungen 
Tm( 1) =0, rn (1) = (I), 


=.0,.15 2,5 2.) a(41, 
eee eo ee ig Tanai as ). (41am) 


Die Differentialgleichungen (38) und (40) lassen sich sukzessiv lésen; dabei sind die auftreten- 
den Integrationskonstanten durch die Randbedingungen (39) und (41) festzulegen. Man gelangt 
so zu dem Funktionensystem r,(&), womit sich dann nach der Vorschrift (37) die Durchbiegung r(é) 
der Welle ermitteln 1aBt. Falls x = a der Ort der Scheibe ist, geniigt es, fiir die Berechnung der 
Einflu8zahl « die Durchbiegung an der Stelle € = &, zu kennen, was dann der Fall ist, wenn die 
Werte r,(&,) bekannt sind. 

Anhand der Differentialgleichungen (38,0), (40,0) und der dazugehérenden Randbedingungen 
(39,0), (41,0) stellt man fest, daB ro(€) die Auslenkung einer mit Q belasteten Welle ist, weshalb 


To(F1) = %) Q (42) 
sein mu. Die GréBe r,(&,) verschwindet ebenso wie die Werte Ta(Ep)etG iste da die Durch- 
biegung nicht vom Vorzeichen des Torsionsmomentes abhangen darf. Um r,(£,) zu bestimmen, 
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wird das durch (38) und (39) baw. (40) und (41) vorgeschriebene Rekursionsverfahren angewandt. 
Man findet dabei nach einigen einfachen Zwischenrechnungen 


ry = — 32 fir O<é<&, (43) 
ri(0) = tee (EE + 26:62) (44) 
ri) = #E Gb — 2) (45) 
und danach 
=—GHeOP—3Hf—266—8) fir OSES4, (46) 
10) =—z ate +284), (47) 
rn) =—7ee G—H-B&), (48) 
womit die Differentialgleichungen (38,2) und (40,2) in 
Vt? 1 
V = GRE CP —WHEL SHO +8H—8E) fir OSESE (49) 
bad 1 B/E.) 
rrr __ 1 %Q : 
ete te, Oct Cent 45D (50) fe 
fiir Sy Ee SM 02 | 


iibergehen. Integriert man sie unter Beriicksichti- 
gung der Bedingungen (39,2) und (41,2), so er- 
halt man schlieBlich mit der in Abb.3 dargestellten ,,,| 


Funktion 
g 
D(E,) = s4p (20 — 45 & + 46 2 — 15) (51) : 
7 
den Wert 0 GF b2 83 G4 OF G6 47 OB 49 10 
r(E,) = D(E,) x) Q. (52) Abb. 3. Die Funktion &(£,). 
Nach (37) ergibt sich somit die Wellendurchbiegung an der Stelle € = &, zu 
Wy) = 05 0 OE) uF >) . (53) 
Setzt man darin Q = 1, so erhalt man die gesuchte Einflu®zahl 
& = % [1 + PE) we +--+] (54) 
und durch einen Vergleich dieser Formel mit (8) und (9) die GréBe 
2 = (6) ce) 


als Funktion von é,. 

Bei allen anderen Lagerungsarten von Wellen mit fehlender Axialkraft S$ laBt sich der Zahlen- 
faktor x, grundsatzlich nach derselben Methode berechnen. Analog gewinnt man die Werte %,, 
wenn man die Differentialgleichung (26) mit W = 0 fiir die einzelnen Falle integriert. Aber auch 
fiir Wellen, die sowohl mit {8 als auch mit Wjbelastet sind, kénnen die Faktoren %, und x, mit 
Hilfe der am vorhergehenden Beispiel angewandten Stérungsrechnung ermittelt werden, indem 
man mit dem Doppelreihenansatz ; 

(2) = 2 =o" Tn mS) 
in die Differentialgleichung (26) und die entsprechenden Randbedingungen eingeht und die Glieder 
mit den gleichen Potenzen von yw und 9¢ gleich Null setzt. Man erhalt so Differentialgleichungen 
und Randbedingungen, aus denen sich die Funktionen r,,,,() und damit auch die Durchbiegung 
r(E,) am Ort der Scheibe leicht berechnen lassen. 

Um eine etwas grobe Abschatzung der Faktoren ~, und %, vornehmen zu kénnen, betrachtet 
man die Formel (10), aus der ersichtlich ist, daB die Klammer (1 — x 0? — xu? —*+*) = 0 wird, 
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wenn «—> co geht. Das bedeutet, daB die Gleichung 
1—x,e?— *,w = 90 (56) 


in erster Naherung die Knickbedingung fiir eine mit einer Kraft P und einem Moment W belastete 
Welle ist. Allerdings wird die Bedingung (56) nur brauchbar sein, wenn die Entwicklung (10) rasch 
konvergiert; dann aber kénnen, sofern nur die Knickwerte von P und W bekannt sind, die Fak- 
toren x, und x, aus der Gleichung (56) berechnet werden. 

In Tabelle 1 sind fiir einige Lagerungsarten schon bekannte! Zahlenwerte von x, fiir WS 0 
und x, fiir 38 = 0 zusammengestellt. Darin sind die aus den Knicklasten gewonnenen Werte durch 
einen Stern gekennzeichnet. Ihre Abweichung von den exakt berechneten Faktoren kann als MaB 
fiir die Konvergenz der Entwicklung (10) angesehen werden. Ferner sind in Tabelle 1 neben den 
Faktoren x, fiir semitangentiale Torsionsmomente auch die entsprechenden Faktoren fiir axiale 
Torsionsmomente YW = (W, 0, 0) aufgefiihrt. 


Tabelle 1. Die Faktoren x, und x, fiir verschiedene Belastungsarten 


W (axial) 


Lagerung 


P (axial) W (semitang.) 


Belastung TSS Xo as 


ny 


0,0219 | 0,0253 


0,0250 0,0253 0,0124 0,0094 0,0124 


LZ LEE: 0,0063 0,0034 0 
P(E,) 
ida (Abb. 3) 


4, Eine partikulare Lésung der Bewegungsdifferentialgleichungen. Den folgenden Betrachtungen 
liegen die nicht gekoppelten, linearen Differentialgleichungen (22) zugrunde. Ihre allgemeinen 
Lésungen setzen sich jeweils aus der allgemeinen Lésung der entsprechenden verkiirzten Differen- 
tialgleichung [die aus (22) durch Nullsetzen der rechten Seiten hervorgeht] und irgendeiner Parti- 
kularlésung der nicht verkirzten Differentialgleichung additiv zusammen. Da in der vorliegenden 
Untersuchung nur kritische Zustande interessieren, geniigt es zu wissen, fiir welche Parameter- 
werte ¢, gq die Lésungen mit der Zeit unbeschrankt anwachsen, wozu die Kenntnis ihrer expliziten 
Form nicht unbedingt notwendig ist. 


Zunichst hat man also die Aufgabe, zwei Integrale 7(r), C(t) der nicht verkiirzten Differential- 
gleichungen zu finden. 


Unter der bereits getroffenen Voraussetzung, daB die Axialkraft P und das Torsionsmoment W 
im Vergleich zu ihren Knickwerten klein sind, wird stets ¢ < 1 sein. Daher liegt es nahe, die beiden 


Gleichungen (22) mit Hilfe der Stérungsrechnung zu lésen, indem man die Lésungen 7/(t) und C(t) 
nach Potenzen von ¢ entwickelt. Man setzt also, wenn man sich vorerst nur auf die erste der beiden 
Gleichungen (22) beschrankt, 


7i(t) =5 e, Mn(t) + (57) 


1 Siehe C. B. Biezeno und R. Grammel, a. a. O. S. 180; H. Ziegler, Z. . Math. Ph : 
H. Ziegler, Ing.-Arch, 20 (1952) S. 390. Sin SE er i 8S aloe 


XXV. Band 1957 Keller: Neue kritische Drehzahlen von einfach besetzten Wellen 719 


geht damit in die Gleichung (22a) fiir (rt) ein und verlangt, daB die Koeffizienten gleicher Potenzen 
von € verschwinden. Auf diese Weise erhalt man die Gleichungen 


Pe +10 =p cos |= + oe), (58) 
¢ oa ++ 1, = No COST — Wp cosT cos — + He Me (59) 
g oo = Hn—1 COST (9 282 Fee) (60) 


aus denen sich die 7,, prea berechnen lassen. Setzt man die Funktion (zr) (23) in die Gleichung 
(58) ein, so geht diese nach einer kurzen Zwischenrechnung, in der man sin } = 9 und cos } = 1 
setzt, weil voraussetzungsgemaB das schwankende Antriebsmoment und damit auch ?% selbst als 


klein anzusehen ist, in 
(L—s)r+8]} (60 


& ya 1 2 
? a + No = p cos as so Poe 6 fein ( ++ ah -}- p.| — sin 
iiber. Man erkennt nun sofort, daB 7, fiir y? = q? und »? = gq? (vy -- 1)? mit der Zeit unbeschrankt 
anwachst. Driickt man in diesen beiden Resonanzbedingungen die Gréfe g mittels der Abkiir- 
zungen (20) und (21) in m und @, aus, so erhalt man in 


oO = Oy, (62) 
und 
ees Ok 
o = (63) 


die allgemein bekannten kritischen Drehzahlen erster und zweiter Art. Sie werden durch ein mit 
der Zeit lineares Anwachsen der GréBe Mo, die ja in der allgemeinen Lésung 7 als Summand auftritt, 


hervorgerufen. 
Die allgemeine Lésung der Gleichung (61) findet man fiir y? = q? und »? + q? (vy + 1)? leicht zu 


iy = C, c08 = + Cy sin + P* ae Mla J ( +5)F +8, 
iL 
(1s) + A]: 


— Uh v 20°89 v®— q?(vy+ 1) 

yay | i 
20° »—q(v»—1)? 

Dabei kénnen die beiden Integrationskonstanten C, und C, willkiirlich festgelegt werden, da 1 eine 

beliebige Lésung der Gleichung (22a) sein darf. Man setzt daher in (64) C, = C, = 0, so daB die 

Gleichung (59) fiir 7, in 


d*n yp v? T W ay 1 ; 1 
py gn — pe ee ee moe 
dt? tM he heuer 207@ y—@g@ovtle a v t+ B, 


sin 


(64) 


sin 


q 


Y ay i 1 t 
Ber ara, gas (0 ry +B, —y c085 (65) 
—, 28 5° (si ( >t 7 B, sin ( sa + |) cos T 


tibergeht. Formt man die rechte Seite dieser Gleichung noch etwas um, so erhalt man 
q an ng ean 5 ta eos ( a ot + cos ( aE 
ate wae In (> +A) + om? tet 6,|| (66) 

ix 0 y2 fe D7] sin Ga — sin (2 ets st a | ; 


Hier zeigt sich nun eine neue Resonanz, falls y? = q? (2 v + 1)?, oder in den Drehgeschwindigkeiten 
ausgedriickt, 


os Ok 67 
Se re (67) 


ist. 
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Analog zur Berechnung von 7) kann auch die Berechnung der Lésungen 7, (n = 1, 2,. <<) gex= 
folgen. Es sind dabei stets Differentialgleichungen vom Typ (61) zu lésen, wobei sich Hiys eee um- 
fener Ausdriicke ergeben. Man iiberblickt jedoch, ohne die Lésungen 7, explizit anzuschrei- 
bear bei welchen Drehgeschwindigkeiten irgendeine von ihnen mit der Zeit unbegrenzt anwachst. 
Es treten namlich, wie man aus dem Aufbau der Gleichung (66) ersieht, in der Lésung 7, fiir y? == q P 
y2 +e g? (y+ 1)? und v? + q? (vy + 2)? Glieder mit sin [(2 + 1/y)t + ,] auf, weshalb nach (60) in 
der Differennaleiocuunes: fin ny auf der rechten Seite Terme mit dem Faktor sin [(2 + I/v)t + f,] 
cost vorkommen, von denen sich Glieder mit sin [(3 + 1/y)7 + 6 ] abspalten lassen. Demnach 
miissen auch wieder in der Lésung 7 mit sin [(3 + 1/»)t + B,] harmonisch pulsierende Anteile 
vorhanden sein. Setzt man diese Uberlegung fort, so erkennt man, da auf den rechten Seiten der 
Differentialgleichungen fir die 7,1 (n = 1, 2,3,...) immer Glieder mit sin [(n + 1/y)t + B,] 
stehen, die, wie ein Blick auf (60) zeigt, dann Resonate bewirken, wenn 7? = q? (nv + 1)? ist. Das 
besagt, daB die Drehgeschwindigkeiten 


o = — (n= 1, 2537063) (68) 


kritisch sind. Formal entsprechen diese den seither wohlbekannten kritischen Drehzahlen zweiter 
Art, die sich bei einem durch eine Fourier-Reihe mit der Grundfrequenz » m darstellbaren An- 
triebsmoment und nicht schwankenden P und W ergeben. Bei den hier abgeleiteten kritischen 
Drehzahlen (68) wurde jedoch ein mit der Frequenz y @ rein harmonisch pulsierendes Antriebs- 
moment und ein mit derselben Frequenz vy wm schwankendes P und W zugrunde gelegt. Der Unter- 
schied zwischen beiden Arten von kritischen Drehzahlen wird noch deutlicher, wenn man bedenkt, 
daB es bei nicht pulsierenden P und W zu jeder Frequenz » eines harmonisch schwankenden 
Antriebsmomentes nur ein Dublett von kritischen Drehzahlen zweiter Art gibt, bei ebenfalls mit 
der Frequenz yw schwankenden P und W jedoch unendlich viele. 
Ist das Antriebsmoment konstant, also ew = 0, so erhalt man mit a, = 0 aus (61) 


T 
q ot eal Ho aye COS oa (69) 
und aus (66) 


d*y yp gr l 1 
GP =) Se ai =e ae eos (1 +5)t +e08 (15 }z 4 


v 


(70) 


Die Lésung 77, besteht daher aus Gliedern mit den Faktoren cos (1 + 1/y)t, und man sieht aus dem 
Aufbau der Gleichungen (60), daB jetzt auf den rechten Seiten der Differentialgleichungen fiir die 
7, (n = 0,1, 2,...) immer Glieder mit cos (n + 1/v)7 auftreten, die ebenfalls fiir y? = q? (nv + 1)?, 
(n = 0,1, 2,...) Resonanz erzeugen. Daher sind die Drehgeschwindigkeiten (62) und (68) auch 
dann kritisch, wenn das Antriebsmoment konstant ist und nur eine periodische Langskraft oder 
ein periodisches Torsionsmoment auf die Welle wirkt. Trotzdem sollen im folgenden auch die Dreh- 
zahlen (68) als kritische Drehzahlen zweiter Art bezeichnet werden. 

Genau wie die Gleichung (22a) laBt sich auch die Gleichung (22b) behandeln. Man findet auch 
da, daB sich bei den Drehgeschwindigkeiten (62) und (68) Resonanz zeigt. 

Wenn das Antriebsmoment der Scheibe zwar periodisch, jedoch nicht, wie in der vorhergehen- 
den Untersuchung, harmonisch wire, stiinde auf der rechten Seite von (23) eine Fourier-Reihe mit 
der Grundfrequenz yw, und es wirden dann in allen Summanden ¢” m,, (= 0,1, oor. eer 
Summe (57), also auch bereits in 7, Glieder auftreten, die bei den Drehgeschwindigkeiten 


(ni OTS ae) (71) 


mit der Zeit unbegrenzt anwachsen wiirden. Dabei ware my die Anzahl der Schwankungen der 
héheren Fourier-Glieder je Umlauf. 

Vom praktischen Gesichtspunkt aus ist noch die Gefahrlichkeit eines kritischen Zustandes 
wichtig; sie 1aBt sich mit Hilfe der Ausweichgeschwindigkeit der Scheibe abschatzen. Um das fiir 
die Drehgeschwindigkeiten (68) zu tun, hat man die Gleichungen (58), (59) und (60) fiir 7,, und die 
entsprechenden Cieichenees fiir ¢,, bei den Drehgeschwindigkeiten (68) zu integrieren und findet 


dann, da®B der Fahrstrahl r, = LV 4+ des Scheibenschwerpunktes im wesentlichen mit einer 
Geschwindigkeit anwachst, die fiir 


: 1 
O == Oy, durch T,= 37 o> (72) 
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fiir 
ee 7, 
Die SNR == O, 
durch 
yo as 26 Gy a1 
eee 2PF02O0  4a,760 ee) 
und fiir 
2s Ok Ok 
< ois Sila) (=n Wad “1.7m = 12. 3,0...) 
durch 
ed, (vt ly fe\r : wp 
ce day ro a of cia anr ae 
oder auch * 
ae 
tie Se 
: @% @ +1), ew[@t) 45 (n+1)y +1 74 
is 4 wp v? O 2 n! (n+ m+1)y+2 (74) 
is 


festgelegt ist. Wie man leicht abschatzt, bleiben die Ausweichgeschwindigkeiten (74) auch fiir be- 
liebig gro8e n klein, und bei einem Vergleich der Geschwindigkeiten (72), (73) und (74) sieht man so- 
fort, daB die kritischen Drehgeschwindigkeiten erster Art weitaus gefahrlicher sind als die zweiter 
Art. Diese sind wohl kaum merkbar, da sie durch eine stets vorhandene Dampfung unterdriickt 
werden diirften. 


5. Die allgemeine Lésung der verkiirzten Bewegungsgleichungen. Nachdem die kritischen 
Zustande erster und zweiter Art durch die partikularen Integrale der Gleichungen (22) beschrieben 
wurden, bleiben nun noch die zu (22) gehérenden verkiirzten Differentialgleichungen vom Mathieu- 
schen Typ 

‘ d2n* 


a + (be cost) 9” = 0, (75a) 
2 at —_ 
Uae Sore (1 — ecost)6* = 0 (75b) 


zu untersuchen. Da beide Gleichungen identisch sind, geniigt es, sich auf eine der beiden, etwa die 
erste, zu beschranken. Die Lésungstheorie der Mathieuschen Differentialgleichungen ist voll- 
standig entwickelt!; es brauchen daher im folgenden nur einige bekannte Ergebnisse wiederholt 
zu werden. 

Die allgemeine Lésung der Gleichung (75a) laBt sich in der Form 

+00 +00 
je — C eet 5 is emt ae C, eet > €., enim (76) 

m = —oo m =—oco 
darstellen, wobei C, und C, beliebige Integrationskonstanten sind. Anhand von (76) kann man 
nun aussagen, unter welchen Umstanden die Lésung 7*(z) stabil ist, d. h. fiir unbegrenzt wachsen- 
des t beschrankt bleibt, und unter welchen Umstanden sie instabil ist, d. h. fiir unbegrenzt wachsen- 
des t gegen + co geht. Dafiir maBgebend ist der sogenannte charakteristische Exponent o, der 
reell, imaginar oder komplex sein kann und eine Funktion der beiden reellen Parameter ¢ und q ist. 

Es ist zweckmaBig, die folgenden drei Falle zu unterscheiden: 

a) Ist o eine von Null verschiedene, reelle Zahl a, also g = a == 0, oder ein komplexer Wert von 
der Form a + im/2 (m = 1, 2,...), so ist die Lésung instabil. 

b) Ist ¢ = 7} rein imaginar und 6b von Null und m/2 (m = 1, 2,.. .) verschieden, so ist die Lo- 
sung stabil. Sie ist auBerdem periodisch, wenn b rational, jedoch nicht periodisch, wenn b irratio- 
nal ist. 

c) Ist o rein imaginar und besitzt die in b) ausgeschlossenen Werte ¢ = 0 oder o =i m/2, so 
bilden die in (64) auftretenden partikularen Integrale kein Fundamentalsystem mehr. In diesem 
Grenzfall 1a8t sich noch eine weitere Lisung gewinnen, mit der sich die allgemeine Lésung in der 


Form 


n*(t) = Cy ar) + Colt malt) + x2(0)] (77) 


1 Siehe z. B. N. W. McLachlan, Theory and Application of Mathieu Functions, Oxford 1951; E. Mettler, 
Mitt. Forschungsanst. Gutehoffnungshiitte 8 (1940) S. 1. 
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schreiben la®t. Die Funktionen ,(t) und 7,(t) sind zwar periodisch mit der Periode 2 7 oder 4.7, 
je nachdem, ob m gerade oder ungerade ist, jedoch ruft der Faktor 7 ein mit der Zeit unbegrenztes 
Anwachsen von 7)* hervor, weshalb die Lésung (77) als instabil bezeichnet werden muB.t 

Weitere Werte von o brauchen nicht betrachtet zu werden, da, wie aus der noch folgenden Ent- 
wicklung (78) ersichtlich ist, der charakteristische Exponent o bei reellen Parametern ¢€ und q nur 
die unter a) bis b) aufgefiihrten Werte annehmen kann. 

Die in a) und b) gemachten Aussagen lassen sich nach entsprechenden Umformungen der Lé- 
sung (76) leicht bestatigen. Zu der Lésung (77) sei bemerkt, daB sie als Grenzfall — er wird spater 
noch etwas naher dargelegt — fiir die vorliegende Untersuchung kaum eine Bedeutung hat. 

Nachdem iiber das Stabilitatsverhalten der Lésungen ein Uberblick gegeben wurde, ist nun fest- 
zustellen, fiir welche Werte der Parameterpaare ¢, q der Exponent o die unter a) vorausgesetzten 
Werte annimmt. 

Bei geniigend kleinem ¢ gilt? fiir o 


PA 
7 sin — 
q 


20 
Co) 220 = cos t- 6% 
] q °° 2q(4—9') (78) 
2 
q (60 — 35 q+ 2q) sin = 4 (4 — q?) (1 — q2) x2 cos 7 

1a oat eee aE = SES 

ie 32 ¢ (4— ge (l— #") 
(Die Voraussetzung ¢ <1 war schon zum Auffinden einer partikularen Lésung der unverkiirzten 
Differentialgleichung nétig und stellt daher keine weitere Einschrankung des Problems dar.) Fir 


é 


| | | So) ee eee | | pena 
a OF O Bb WwW @ # 6: B.-20 a 2 a6 


Abb, 4. Die ersten fiinf Instabilitatsbereiche. 


die Werte q = 2/m (m = 1, 2,...) liefert die Entwicklung (78) zunachst unbestimmte Ausdriicke, 
die jedoch mit der de [’Hospital’schen Regel leicht bestimmt werden kénnen. So ergibt sich fiir 
q=2 


cas 7? , 4 2(1385—12n%) 

Co) 20 es ian oe | 
und fiir gq = 1 (79) 

ae 

ty = 4s 580: f 

Coj2xao=1+e Tiare 
Hieraus ist ersichtlich, daB in der Nihe der Stellen 

2 
g=—) (uel, 2) (80) 


der Betrag von Woj 27 gréBer als eins ist; denn die rechte Seite der Gleichung (78) wird haupt- 
sichlich durch das Glied cos 2 2/q bestimmt, wahrend die weiteren Glieder nur noch kleine Korrek- 
turen liefern, die allerdings gerade bewirken, daB ©oj 20 das Intervall (—1, +1) in der Um- 
gebung der Punkte (80) verlaBt. Wenn aber | Co] 270 | > 1 ist, kann o nur entweder eine von Null 
verschiedene reelle Zahl oder eine komplexe Zahl von der Formo = a + i m/2 (a = 0,m =1,2,.. 5 
sein, so daf} in diesem Fall nach a) die Lésung instabil ist. Dagegen ist die Lésung stabil, wenn 


1 In der Theorie der Differentialgleichungen wird die Lésung (77) meist als stabil bezeichnet, weil jede noch 


so kleine, geschwindigkeitsproportionale Dampfung dafiir sorgt, daB bei unbegrenzt wachsendem Tt die Losung 
7)* beschrankt bleibt. 
2 E. Mettler, a. a. O. S. 5. 
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il > Co} 220 > —1 ist; denn dann kann o uur die in b) vorausgesetzten Werte annehmen, Ist 
peblicBlich | Coj 2 ao | = 1, so liegt der in c) besprochene Grenzfall vor; die zu diesem gehorenden 
Wertepaare €, q bilden in einer (e, q)-Ebene sogenannte Grenzkurven, die gemaf Abb. 4 die Ebene 
in stabile und instabile Bereiche aufteilen, wobei die zuletzt genannten in den Punkten ¢ = 0, 
7 2/m (m = 1, 2,...) entspringen. Von den unendlich vielen Instabilitatsbereichen, die sich in 
q = 0 haufen, sind in Abb. 4 fiinf schraffiert gezeichnet. Dabei wurden beim Aufzeichnen der 
Kurven zahlreiche Punkte aus den Inceschen Tabellen! entnommen und weitere durch entspre- 
chende Reihenentwicklungen ermittelt. 

Fir die Geschwindigkeit, mit der die Lésung (76) innerhalb der Instabilitatsbereiche anwachst, 
ist der Realteil a des charakteristischen Exponenten o maBgebend. Seine GriBe kann man fir 
einen gegebenen Punkt (¢, q) mit Hilfe der Gleichung (78) bestimmen, indem man o — a + i m/2 
(m = 1, 2,...) setzt, wodurch die linke Seite der Gleichung in (—1)"Coj 27 a ttbergeht. FaBt man 
é als Parameter auf, so lat sich in einem Diagramm der Wert a als Funktion von q darstellen. Ein 
solches, von E. Mettler entworfenes?, Diagramm ist in Abb. 5 wiedergegeben. In ihm ist statt a die 
GroBe 22 aq als Ordinate gewahlt, was sich spater als vorteihaft erweisen wird. Abb. 5 zeigt, daB 
sich die halbovalen Kurven jeweils um die Punkte q = 2/m der Abszissenachse gruppieren und mit 


(ag) 


max 


Lag 


l 
19 40 a1 


a2 ads 
Abb. 5. Die GréBe 2 2a q, ees von den Parametern ¢ Abb. 6. Die GréBe (2 74 a) ee fiir die ersten drei 
und q. 1} 


Instabilitatsbereiche, abhangig von e. 


abnehmendem ¢ schlieBlich auf sie zusammenschrumpfen. Fiir die Schnittpunkte der Kurven mit 
der q-Achse trifft der Fall c) von Seite 81 zu; sie entsprechen also den auf den Grenzkurven der 
Abb. 4 liegenden Punkten. Bemerkenswert ist noch, da die Halbovale bei gleichbleibendem « 
um so kleiner werden, je naher sie an der 2 7 a q-Achse liegen. Das kommt auch zum Ausdruck, 
wenn man von den einzelnen Kurven nur die Maxima betrachtet. Zeichnet man sie in Abhangig- 
keit von ¢ auf, so entsteht Abb. 6. Dabei ist, um das Einmiinden der Kurven in den Ursprung 
deutlich herauszustellen, fiir die Ordinate eine quadratische, fiir die Abszisse jedoch eine lineare 
Teilung gewahlt. Auf diese Weise ist gut zu erkennen, daf die zum ersten Instabilitatsbereich ge- 
hérende Kurve schon fiir sehr kleine ¢ betrachtliche Werte fiir (2 7 a q),,,, aufweist. Die Kurven 
des zweiten, dritten und der nicht gezeichneten héheren Instabilitatsberiche verlaufen alle rechts 
von der des ersten Bereichs und schmiegen sich immer mehr der ¢e-Achse an. Damit ist die Frage 
nach der GréBe von a geniigend beantwortet. 

Nachdem so einige fiir die vorliegende Untersuchung wichtige Kigenschaften der Lésungen von 
(22) erlautert sind, soll noch einiges iiber die mathematische Behandlung der allgemeinen Diffe- 
rentialgleichungen (17) bemerkt werden. 

Die allgemeinen Lésungen von (17) kénnen auf dieselbe Weise gefunden werden wie die der 
Gleichungen (22), wenn nur die Parameter y,, y, <A sind, was in praktischen Fallen wohl immer 
zutreffen wird. Man hatte dann, um ein partikulares Integral zu finden, die Stérungsrechnung mit 
zwei Parametern durchzufihren, d. h. die Lésung ware zuerst nach dem einen Parameter zu ent- 
wickeln und die dabei auftretenden Glieder nach dem anderen. An dem Ergebnis wiirden dann 
wieder die kritischen Zustande erster und zweiter Art zu erkennen sein. Die zu (17) gehorenden 


1 N. W. McLachlan, a. a. O., Appendix II. 
2 E. Mettler, a. a. O.S. 6. 
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verkiirzten Gleichungen 
d?y* 
dt? 
acGe : 
qa + (A+ yy cost + 72 60s OPAC AU 


+ (A + y, cost + V2 cos 2t)4* =0, 


(81) 


sind wie die Gleichungen (75) vom Hillschen Typ. Im Gegensatz zu diesen ist die allgemeine Lé- 
sung der Gleichungen (81) nicht in allen Kinzelheiten bekannt. Es wurden jedoch schon Stabilitats- 
karten berechnet!, die gestatten, in einem (A, y,, 2)-Raum die instabilen Bereiche abzugrenzen. 


6. Der kritische Zustand dritter Art. Wie die Resonanzerscheinungen bei bestimmten Drehzahlen 
der Welle die kritischen Zustande erster und zweiter Art hervorrufen, so bewirkt auch die Instabi- 
litat einen kritischen Zustand. Dieser wird zum Unterschied von den beiden erstgenannten als 
kritischer Zustand dritter Art bezeichnet. Seine Ursache ist darin zu suchen, daB die pulsieren- 
den Belastungen P(t) und W(t) bei bestimmten, von der Drehgeschwindigkeit @ abhangigen Fre- 
quenzen Q nach einer von auBen kommenden Stérung — der Anla® zu einer solchen wird praktisch 
immer gegeben sein — Biegeschwingungen aufschaukeln, deren Amplituden theoretisch unbegrenzt 
anwachsen. Diese Erscheinung kann jedoch auch bei Wellen auftreten, die sich nicht drehen. Sieht 
man namlich von der Drehbewegung ab, so andern sich zwar die rechten Seiten der Gleichungen 
(22), die dazugehérenden verkiirzten Differentialgleichungen sind jedoch nach wie vor mit (75) 
identisch und fiihren, falls nur eine pulsierende Belastung P(t) und W(t) vorhanden ist, zu Instabili- 
tatserscheinungen. Im Gegensatz dazu treten die kritischen Zustande erster und zweiter Art nur 
bei sich drehenden Wellen und dem Vorhandensein einer Unwucht auf. Aber nicht nur in der Ur- 
sache unterscheiden sich die kritischen Zustande wesentlich, sondern auch in ihrer Gefahrlichkeit. 
Obwohl beide auf dem Anwachsen gewisser Biegeschwingungen beruhen, erfolgt dieses Anwachsen 
bei den kritischen Zustanden erster und zweiter Art proportional zur Zeit, wie bei allen Resonanz- 
erscheinungen, bei den kritischen Zustanden dritter Art jedoch exponentiell mit der Zeit. Das 
eine Mal sind also die Ausweichgeschwindigkeiten der Scheibe zeitunabhangige GréBen und kénnen 
somit direkt als Maf fiir die Gefahrlichkeit eines kritischen Zustandes angesehen werden, das andere 
Mal ist jedoch die Ausweichgeschwindigkeit eine Exponentialfunktion der Zeit; hier ist daher als 
MaB fiir die Gefahrlichkeit der Realteil a des charakteristischen Exponenten anzusehen, der sich 
bei bekanntem ¢ und q aus der Abb. 5 entnehmen oder auch mittels der Gleichung (78) berechnen 
JaBt. Es scheint demnach die Gefahrlichkeit eines kritischen Zustandes dritter Art erheblich gréBer 
zu sein als die der anderen kritischen Zustande, da ja die Exponentialfunktion schlieBlich steiler 
als jede lineare Funktion anwachst. Ist jedoch die Anfangsamplitude der Biegeschwingung sehr 
klein, so wird die Schwingung zu Beginn nur sehr langsam aufgeschaukelt werden. In diesem Fall 
diirfte sich z. B. das Durchfahren eines kritischen Zustandes dritter Art kaum als gefahrlich er- 
weisen. 

Weiter ist auch die Frage, inwieweit sich eine in Wirklichkeit stets vorhandene Dampfung auf 
die Instabilitat auswirkt, von Bedeutung. Um hier eine Antwort zu finden, miissen zunachst die 
Bewegungsgleichungen durch gewisse Dampfungsglieder erganzt werden. Wollte man dabei die 
Abhangigkeit der Dampfungskrafte einerseits von der 4uBeren, andererseits von der inneren Damp- 
fung (Werkstoffdampfung) genau angeben, so wiirde man auf erhebliche mathematische Schwierig- 
keiten stoBen. Die Rechnung wird einfacher, wenn man die Dampfungskrafte als proportional zu 
den Geschwindigkeiten 7 und ¢ ansieht. Man wird dabei die Proportionalitatskonstanten so be- 
stimmen, da durch die wirklichen und die angenommenen Dampfungskrafte je Schwingung die- 
selbe Energie verbraucht wird, und kann somit erwarten, daB der auf diese Weise beriicksichtigte 
Einflu8 der Dampfung zu Ergebnissen fiihrt, die wenigstens qualitativ mit der Wirklichkeit iiber- 
einstimmen., 


Es werden also die Gleichungen (22) durch Glieder von der Form at B dnp und 2g BaeD 


dt OCS ae 
bei 6 eine Dampfungskonstante ist, zu 
d*7 2qg° bd 
@ 7a e +2 + (L—e cost) np =y (1 — € cos) (co 5 +9), (82a) 
dl 2 q° B dl ; : 
ae ta ae t (le cost) ly = y (1 —¢ 087) (sin ral 8) (82b) 


' K. Klotter u. G. Kotowski, Z. angew. Math. Mech. 23 (1943) S. 154. 
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erganzt. Fiir die erste der beiden Gleichungen lat sich auch mittels des Ansatzes (57) ein parti- 
kulares Integral gewinnen, das wieder die kritischen Zustande erster und zweiter Art beschreibt. 
Die Ergebnisse weichen jedoch hier von denen des Abschnittes 4 etwas ab, weil dort bei Resonanz 
die Ausschlage unbegrenzt anwachsen, wahrend sie hier nur bis zu einem bestimmten endlichen 
Wert zunehmen. Ferner werden die kritischen Drehzahlen etwas kleiner, sobald eine Dampfungs- 
kraft wirkt. Jedoch diirfte dieses Absinken bei normalerweise vorkommenden Dampfungen kaum 
merkbar sein. Da der Unterschied zwischen Resonanzerscheinungen bei dae eelacayiten und ge- 
dampften Schwingungen wohlbekannt ist, eriibrigt es sich, diesen noch Si rekhondes zu betrachten. 
Um die bei kritischen Zustanden dritter Art durch die Dampfung Ter creeitneh Anderungen 
erkennen zu kénnen, muB noch die allgemeine Lésung der zu (82a) gehérenden verkiirzten Diffe- 
rentialgleichung . 


ody , 2478 dnp 


Tepe G7 7, 1 8 Ost) pp = 0 (83) 
aufgesucht werden. Diese geht mittels der Transformation 
ae 
nb(t) =e * d(z) (84) 
in die Gleichung 
of Ga + (scoot — ia |d = 0 (85) 


tiber, in der die erste Ableitung der unbekannten Funktion fehlt. 

Bevor die Gleichung (85) weiterbehandelt wird, soll noch der Dampfungsfaktor § zahlenmaBig 
abgeschatzt werden. Dazu bringt man ihn mit dem logarithmischen Dekrement A, dessen Groen- 
ordnungen eher bekannt sind, in Verbindung. Es ist definitionsgema8B 


e—Bbt 2 IG 
A=In ae == 8 (86) 
- (+5, 
und damit 
p a 
pe toeaas SD) 


Handelt es sich um Wellen aus Metall, so ist! immer A < 10-*. Folglich ist der Ausdruck (87) sehr 
viel kleiner als eins. Schreibt man mit 


iy ny ee und C tae Pipes 


die Gleichung (85) um in 
Che) 
Giza + (l— 24, cost) d = 0, 


so ist sie formal mit (75a) identisch. Nach (76), (84) und (87) lautet also die allgemeine Lésung der 
Gleichung (83) 


ee ak t +00 -+-0o 
{ipa eOn (c CE Tagg SN Real Aaa ek ON cage Gas AE aR in . (88) 
m= — oC m= —oo 
wobei sich der Wert des charakteristischen Exponenten o in (88) von dem in (76) wegen q, ~ q und 
€ & é, nur vernachlassigbar wenig unterscheidet. Fiir die Stabilitat ist also jetzt nicht mehr der 
Realteil a des charakteristischen Exponenten o sondern die GréBe a— A/2 aq mabgebend. Es 
herrscht namlich nur dann Instabilitat, falls a — A/2 2 q > 0 ist. Damit lassen sich nun fir jedes 
gegebene A die Instabilitatsbereiche eingrenzen. Man zeichnet in Abb. 5 die waagrechte Gerade 
2a aq = A cin und markiert ihre Schnittpunkte mit den bereits gezeichneten Halbovalen. Sodann 
iibertragt man diese Punkte in eine (¢, q)-Ebene und verbindet sie zu einer stetigen Kurve, die dann 
die Stabilitatsgebiete von den Instabilitatsgebieten trennt (siehe Abb. 7). Da sich die so gewon- 
nenen Grenzkurven fiir nicht zu kleine ¢ und kleine A kaum von den urspriinglichen Grenzkurven 
(A = 0) unterscheiden, erkennt man an Abb. 5. Fiir kleine ¢ ergibt sich jedoch eine merkwiirdige 
‘Erscheinung. Jede Gerade 2 aq = A wird in Abb. 5 von bestimmten Halbovalen beriihrt. Die 
au diesen gehérenden Parameter sind sogenannte Schwellwerte ¢), denn alle Halbovale mit Para- 


1 F, Forster und W. Kester, Z. f. Metallk. 29 (1937) S. 116. 
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metern ¢ < é liegen unter der Geraden 22 aq = A. Sobald also eine Dampfung vorhanden ist, 
reichen demnach die Instabilitatsbereiche nicht mehr bis zur g-Achse hinab, sondern nur noch bis 
in die Hihe ¢, (siehe Abb. 7 und 8). Fiir die ersten drei Instabilitatsbereiche lassen sich bei ge- 


gebenem A die dazugehdrenden ¢)-Werte aus Abb. 6 entnehmen. Man hat dazu nur die Gerade 
(22044) max = A mit den dort gezeich- 


vents neten Kurven zu schneiden und er- 
‘ \ halt dann die Schwellwerte als Ab- 
szissen der Schnittpunkte. 

In Abb. 7 und 8 sind fir A = 0 
und A = A* = 0,01 die Grenzkurven 
der ersten zwei Instabilitatsbereiche 
gezeichnet. Mit ihnen 1aBt sich der 
Einflu8B der Dampfung insofern ab- 
schaitzen, als man annehmen darf, 
daB bei Biegeschwingungen von Me- 
tallwellen immer A < A* ist und 
somit auch die zu A gehérenden Wer- 
te ¢ kleiner sind als der zu A* ge- 
hérende Schwellwert e§ (Abb. 6, 7 
und 8). Da sich auBerdem zwei Grenz- 
kurven nicht schneiden kénnen!, so 
miissen fiir Dampfungen mit 0 < A 
< A* alle Grenzkurven in den einfach 
schraffierten Gebieten von Abb. 7 

poaeeert l 7 bzw. 8 liegen. 

be oe ls dl a Berechnet man die Schwellwerte 
Abb. 7. Der erste Instabilitatsbereich. ex tir) disscrstenuires Instabilitates 

bereiche, so erhalt man nacheinander ef = 0,0064, e* = 0,11 und ¢e> = 0,26. Es wird also der 

Schwellwert im ersten Instabilitatsbereich sehr klein sein, wahrend fir die nachstfolgenden 

Instabilitatsbereiche ¢, stark anwachst, was auch in Abb. 6 zum Ausdruck kommt. 
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Abb. 8. Der zweite Instabilitatsbereich. Abb. 9. Rotor und Belastung der Welle 


Nachdem man nun iibersieht, wie sich die Dampfung auf die Stabilitatsverhaltnisse auswirkt, 
sei zusammenfassend festgestellt, daB bei einem gedampften Schwinger nur dann Instabilitat 
méglich ist, wenn € oberhalb eines bestimmten Schwellwertes ¢, liegt. Dieser ist bei dem ersten 
Instabilitatsbereich sehr Klein, weshalb hier schon eine verhaltnismaBig schwache pulsierende Be- | 
lastung ausreicht, um Biegeschwingungen aufzuschaukeln. Bei den ubrigen Instabilitatsbereichen 


1 N. W. McLachlan, a. a. O. S. 80—81. 
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wird €) um so gréRer, der Bereich also um so ungefahrlicher, je naher er an der e-Achse liegt. Fiir das 
_ Anwachsen der Lisung (tr) mit der Zeit ist nicht mehr wie im ungedampften Fall das Produkt 
a {2 maBgebend, sondern die GréBe (274 q—A) w,/2. Das bedeutet, daB die Ausweichge- 
schwindigkeit der Scheibe durch eine Dampfung stets herabgesetzt wird. 


7. Zwei Zahlenbeispiele. Die dargestellten Zusammenhinge sollen nun noch an zwei Zahlen- 
beispielen verdeutlicht werden. 

Kine beiderseits gelagerte Welle vom Durchmesser d = 3 mm tragt in ihrer Mitte einG = 0,1 kg 
schweres Antriebsrad (Abb. 9). Der gleichgroBe Abstand vom Gleitlager A bzw. B bis zur Ein- 
spannstelle am Rad S sei 1/2 = 120 mm. Das Wellenstiick AS wird durch ein pulsierendes Tor- 
sionsmoment W = W, + W, cos 3 wt beansprucht, herriihrend von einem mit derselben Frequenz 
harmonisch schwankenden Antriebsmoment der Scheibe S. Wie groB sind die kritischen Dreh- 
zahlen, wenn W, = 8 cmkg, W, = 0,2 emkg und der Elastizitatsmodul E = 2,15 + 106 kg/cm? ist ? 

Zur Beantwortung dieser Frage geht man von der hier durchaus berechtigten Annahme aus, 
daB die Masse der Welle gegeniiber derjenigen des Antriebsrades vernachlassigbar ist. Nach (16) 
und (18) ergibt sich die kritische Drehgeschwindigkeit «, der nur mit W, vorbelasteten Werte zu 


ee 
Ce 28 . (89) 
Da fiir die gegebenen Zahlenwerte x, u? << 1 ist, kann man in guter Naherung auch 

1 


O, = 0 —— (90) 
V mM % 
setzen. Dabei findet man die Einflu®zahl 
1 B cm 
Xo = 3a Edt — 0,084.22 ke > (91) 


wenn man in Anbetracht der geringen Wellendicke die Gleitlager A und B als starr und damit die 
Welle als beiderseits eingespannt ansieht*; fiir die Radmasse erhalt man m = 1,019 - 10-4 kgsek?/em 
und somit nach (90) den Zahlenwert 
@, = 341,2 sek. (92) 
Die kritischen Zustande erster und zweiter Art treten nun bei den durch die Gleichungen (62) 
und (68) mit »y = 3 bestimmten Winkelgeschwindigkeiten 
Ok 341,2 


Oe wy —= JA1,2 sek und ON Ereracrs Fea =a sek-1 i a eine an Oe) 


auf, aus denen die minutliche kritische Drehzahl erster Art 


Umi. (94) 


min ; 


N=, 3258 
JU 


und die nicht so gefahrlichen, in Tabelle 2 aufgefiihrten kritischen Drehzahlen zweiter Art folgen. 


Tabelle 2 
oo | ee eee ae 5 
1629 651 408 | 216 232 \ Um. 
Me 815 465 325 | 250 203 f min 


Um die kritischen Drehzahlen dritter Art bestimmen zu kénnen, berechnet man zunachst die 
GréBe e. Mit (21), (12) und dem aus der Tabelle 1 entnommenen Zahlenwert x, = 0,0034 erhalt man 


WW, al 5 10-5) (95) 


Da dieser Wert sehr klein ist, konnen nach Abb. 4 nur in unmittelbarer Nahe von q = 2/m (m = 1, 
2,...), also nach (21) bei den Winkelgeschwindigkeiten 


= O, (Se) (96) 


1 Hiitte, Bd. 1, S. 666, 27. Aufl. Berlin 1949. 
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kritische Zustinde dritter Art auftreten. Die den Winkelgeschwindigkeiten (96) entsprechenden 
minutlichen Drehzahlen ergeben sich fiir m = 1, 2, 3 zu 


0172 Um jog6 Ua, 794 Ut (97) 
min min 


min 

Ob bei diesen jedoch tatsichlich Biegeschwingungen aufgeschaukelt werden, hangt noch von der 
Dampfung ab. Bei einem logarithmischen Dekrement A = 0,01 wiirden sich z. B. keinerlei In- 
stabilitatserscheinungen zeigen; denn in diesem Fall ware bereits der Schwell- 
wert des ersten Instabilitatsbereichs erheblich gréBer als ¢ = 1,05- 10~° 
(Abb. 7). Es ist aber anzunehmen, da® fiir Schwingungen mit sehr kleiner 
Amplitude auch A auSerst klein ist und erst bei gréBeren Amplituden einen 
Wert von der GréBenordnung 10-* bis 10-2 annimmt. Es werden sich daher 
in einem kritischen Zustand die Amplituden der Biegeschwingungen nur so 
lange aufschaukeln, bis das logarithmische Dekrement den Wert A = (27 4q) max 
— 1,65: 10-5 erreicht hat, und dann konstant bleiben. Die Drehzahlen (97) 
werden sich demnach nur durch einen etwas unruhigen Lauf der Welle be- 
merkbar machen. 

Als zweites Beispiel wird eine beiderseits gelenkig gelagerte Welle von der 
Lange 1 = 40cm und der Biegefestigkeit EL J = 40 000 kgem? betrachtet, 
die in ihrer Mitte eine Scheibe vom Gewicht G = 4 kg tragt (Abb. 10) und 
mit einer axialen, harmonisch schwankenden Druckkraft P = P, cos4mt 
belastet ist, wobei P, = 10 kg ist. 

Mit der EinfluBzahl 


Xo 


ele 
~ 48 EJ 


und der Scheibenmasse m = 4,077 kg sek?/em erhalt man aus (62), (18) 
und (16) die kritische Drehgeschwindigkeit erster Art 


= 0,0333 = (98) 


Abb. 10, Rotor und 2S. 1 — 
Belastung der Welle. OO}, = WO = = = 85,78 sek 1, (99) 
ma 


die einer minutlichen Drehzahl von 
Umil. 


30 
: MS ed Nae 819 min (10 
entspricht. 
Nach Abschnitt 4 sind hier fiir die Drehgeschwindigkeiten 
(62) 
arpa (e212, 58) (101) 


kritische Zustande zweiter Art zu erwarten. Die zu diesen gehérenden kritischen Drehzahlen sind 
fiir n = 1, 2, 3 in Tabelle 3 angegeben. 

Die Instabilitatsverhaltnisse lassen sich iibersehen, sobald man ¢ berechnet hat. Man findet 
mit den gegebenen Zahlenwerten und dem aus Tabelle 1 entnommenen Wert x, = 0,10 


Tabelle 3 
: Beate ee c= Fy = 0,040. (102) 
n il 2 3 
Aus Abb.7 ist ersichtlich, da®B der Wert 
He 273 | 117 | 74 Um. (102) fiir A = 0,01 zwar den Schwell- 
164. 91 63 min wert ¢9 — 0,0064 des ersten Instabilitats- 


bereichs iibertrifft, jedoch bereits kleiner 

als der Schwellwert ef = 0,11 des zweiten Instabilitatsbereichs ist. Demnach werden also nur 
in einem sehr kleinen Bereich um q = 2 bzw. 

w@ = 29 — 42,89 sek (103) 

Biegeschwingungen aufgeschaukelt. Da aber A wohl stets kleiner als 0,01 sein wird, ist damit zu 

once daB auch noch der zweite Instabilitatsbereich mit dem Fu8punkt q = 1 gefahrlich wer- 

em kann. Mit A = 0,001 ergibt sich z. B. fiir den zweiten Instabilitatsbereich ein Schwellwert 

Ey = 0,03 < e = 0,04. In diesem Fall stellt sich also nicht nur bei der Drehgeschwindigkeit (103) 
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ein kritischer Zustand dritter Art ein, sondern auch in unmittelbarer Nahe der Drehgeschwindigkeit 


@ 


== a = 21,44 sek-!, (104) 
Es folgen daher aus (103) und (104) die kritischen Drehzahlen dritter Art 
N = 2020 _ 419 Um (105) 
min 
und 
30 @, Uni. 
NS NN rare (106) 


Uber die bei diesen zu erwartenden Ausweichgeschwindigkeiten gibt der Exponent 
[(2 % 4 q) max — A] ©,/2 Auskunft. Die Gleichung (79) liefert fiir ¢ = 0,04 und q = 2 den Wert 
(2% 4q)max = 9,063 und fiir ¢= 0,04 und gq = 0,99987 den Wert (2244) ma, = 90,0013. Fiir 
A = 0,001 erhalt man damit im ersten Instabilitatsbereich 


max 


(oe aq), — A] st = 0,846 (107) 
und im zweiten 
[(2 © @ q) maz — A] ~ = 0,00355 . (108) 


Die Amplitude der Biegeschwingung wachst also im ersten Instabilitatsbereich in jeder Sekunde 
auf das e°-4° — 2,32-fache ihres Betrages an und erreicht so bereits nach 8,2 Sekunden das tausend- 
fache ihres Anfangsausschlages. Im zweiten Instabilitatsbereich dagegen steigt der Betrag der 
Amplitude in jeder Sekunde nur auf das e° 355 — 1,00355-fache, so daB es hier etwa 1930 Sekun- 
den dauert, bis sie tausendmal gréfer ist als ihre Anfangsamplitude. 

Da im Maschinenbau kaum Faille vorkommen, die zu viel eréBeren e-Werten fiihren, als sie das 
letzte Zahlenbeispiel geliefert hat, darf man aus diesem wohl folgern, daB auBer dem ersten In- 
stabilitatsbereich nur noch der zweite beachtet werden muB und da wegen der kleinen ¢ beide 

Bereiche fast punktférmig sind und in unmittelbarer Nahe von q = 2 baw. q = 1 liegen. 


8. Zusammenfassung. Wird eine mit einer einzelnen Scheibe besetzte, masselose Welle mit einer 
periodischen Langskraft P, + P, cos Qt oder mit einem periodischen Torsionsmoment W, + W, 
cos {2 t belastet, so ergibt es unter der Annahme, daB bei umlaufenden Wellen die Frequenz 2 mit 
der Drehzahl wm der Welle gemaB der Beziehung 2 = y w gekoppelt ist, auch bei konstantem und 
um so mehr bei rein harmonisch periodischem Antriebsmoment der Scheibe neben der kritischen 
Drehgeschwindigkeit erster Art @ = w, unendlich viele kritische Zustande zweiter Art bei den 
Drehgeschwindigkeiten wm = w,/|nv—+1| (n = 1, 2, 3, ...), falls die Scheibe etwas exzentrisch 
auf die Welle gekeilt ist. Die Gefahrlichkeit der kritischen Drehzahlen zweiter Art ist jedoch viel 
geringer als die der kritischen Drehzahl erster Art. 

Auf8er den genannten kritischen Drehzahlen ergeben sich ohne Dampfung noch unendlich viele 
kritische Drehzahlbereiche in der Nahe des 2/my-fachen der kritischen Drehzahl w,, in denen 
die mit der Scheibe besetzte Welle exponentiell mit der Zeit anwachsende Biegeschwingungen aus- 
fiihrt. Die Breite dieser Bereiche und die Ausweichgeschwindigkeit der Scheibe nehmen mit wach- 
sendem P, bzw. W, zu. Da die Ausweichgeschwindigkeit exponentiell mit der Zeit anwachst, unter- 
scheiden sich diese neuen kritischen Drehzahlen dritter Art wesentlich von den kritischen Dreh- 

_zahlen erster und zweiter Art, bei denen die Scheibe ja mit konstanter Geschwindigkeit ausweicht. 
Dieser Unterschied kommt besonders dann zur Geltung, wenn man noch den Einfluf einer ge- 
_ schwindigkeitsproportionalen Dampfung beriicksichtigt. Eine solche begrenzt bei den kritischen 
_ Zustanden erster und zweiter Art das seitliche Ausweichen der Scheibe, wahrend bei den neu ge- 
- fundenen kritischen Drehzahlen dritter Art die Amplitude der Biegeschwingung nach wie vor un- 
begrenzt anwachst. Allerdings fallen infolge der Dampfung die Bereiche mit héheren Ordnungs- 
_ zahlen m weg, und die iibrigen werden um so kleiner, je gréBer die dazugehérende Ordnungszahl m 
ist. In der Praxis miissen — wie an zwei Zahlenbeispielen gezeigt wurde — ein erster Linie die 
Instabilitatsbereiche beachtet werden, die in der Nahe von 2 w,/y und ,/y liegen. 


| 
(Eingegangen am 1. Juni 1956) 
| Anschrift des Verfassers: Dipl.-Phys. W. Keller, Schwab.-Hall, Auf dem Klingenberg 5 
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Uber die Biegung des orthotropen Plattenstreifens durch Einzellasten 


Von S. Woinowsky-Krieger 


1. Einleitung. Die Theorie der Biegung orthotroper Platten hat in letzten Jahren im Hinblick 
auf ihre Anwendung auf die Berechnung von Hohlplatten, von Briicken-Fahrbahnplatten! ver- 
schiedener Art und auch von Tragerrosten? eine gréBere Bedeutung erlangt. Die im wesentlichen 
von M. Huber herrithrende Theorie orthotroper Rechteckplatten fiihrt auf Momenten- und Quer- 
kraftausdriicke in Form von einfachen Fourier-Reihen, die, wenn es sich um die nahere Umgebung 
stark konzentrierter Einzellasten handelt, nur sehr langsam konvergieren. Zweck der vorliegenden 
Arbeit ist es, die technisch wichtigsten Ergebnisse der Huberschen Theorie in geschlossener Form 
darzustellen, um auf diese Weise ihre Anwendung zu vereinfachen. 


2. Hubersche Lésung fir den freiaufliegenden Plattenstreifen. Die Differentialgleichung der 
Biegeflache w(x, y) einer orthotropen diinnen Platte lautet im Falle einer Einzellast 


7 4 
| = 1 
yy teat oem ae ay OS o 
Hierbei sind B,, B, und H gewisse, von den elastischen Eigenschaften der Platte bzw. des Trager- 
rostes abhangige Steifigkeiten*. Um alle méglichen Verhaltniswerte dieser Steifigkeiten zu er- 
fassen, setzen wir, im wesentlichen nach Huber, 


1a Se 
jek A=|/R Uaiee ) 


DEB 


und weiterhin® 


2 = ze Vu = yp? —1 : wenn fo 1S 


Vy a! bA 5 
Abb. 1. Plattenstreifen mit Einzellast (ay = 7a / io 2 wea fag alt 
im Koordinatensystem. p a6 fe 


(3) 


Die elastische Flache eines freiaufliegenden Plattenstreifens mit einer transversalen Einzellast in 


x = 0, y =7 (Abb. 1) driickt sich dann im Gebiete x > 0 folgendermafSen aus: 


(Bebe 1 Sy 1 _"* —"*\, OAR) 9. WOE SY 
Fall) “weet: a ae > 1 (pe P—— 6 (5 = | sin oe ae (4,1) 
n=1 
IP Te — 1 nux\ ———, WE Ue 3 (OE SE 
Pall 1s W= 9 5B A ale bi )e oA sin ages (4,11) 
n=l 
FHT, po 20s we NT eee oe eee 4,111 
SREB, ai AE F ps ee 


Der Fall If 1aBt sich im wesentlichen auf den Fall einer isotropen Platte zuriickfiihren® und im 
praktisch ebenso wichtigen Falle III erhalt man alle gewiinschten Ausdriicke aus den fiir den Fall I 
geltenden Ausdriicken durch die Substitution 


i 1 1 1 1 1 
a (be ca (5) 
: Siehe z. B. H. Olsen und F. Reinitzhuber, Die zweiseitig gelagerte Platte. Berlin 1950. 
_ = Yel. besonders Y. Guyon, Ann. Ponts et Chauss. (1946) S. 553 und Ch. Massonet, Abh. Intern. Vereinigung 
fiir Brickenbau u. Hochbau 10 (1950) S. 147. 
2 Soon Huber, Probleme der Statik technisch wichtiger orthotroper Platten, insbesondere S. 42, War- 
schau : 


* Naheres tiber ihre Berechnung in verschiedenen Einzelfallen siehe bei K. Girkmann, Flachentragwerke 
s. 475, Wien 1948. 

° Fur manche Umformungen erweist es sich als zweckmaBig im Falle I = oj 2@ und im Falle III 
t= cos 2@ zu setzen; dann hat man « = e— b A/a, B = ew bAla, «’ =b Ala cosa, fp’ =b A/x sin w. 

§ Vel. K. Girkmann, Flachentragwerke S. 487, Wien 1948. 
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mit = y—1. Im folgenden soll daher der formal einfachste Fall I stets als erster behandelt wer- 
den. Fiir die beiden Kriimmungen erhalten wir in diesem Fall die folgenden Reihenausdriicke: 


Bean Pb l en ft | Se ead at 
= ; ee TUMTUCY. 
Ox? = BA? B?— o? > | Aa ioe a a ieee be 


n=1 


Ow Pies 4 i —s aes 
es = FEN On LOB, cae TOMES 
dy? 7? B, p*—? - ae : pee a aed abs 


joel 


Die Biegemomente ergeben sich aus den Kriimmungen mittels der Ausdriicke 


(7) 


mit den Querdehnungszahlen vy, und 7,. Es geniigt somit an Stelle der Momente weiter die Kriim- 
mungen allein naher zu betrachten. 


3. Biegemomente und Schubkriifte in der Nahe einer Punktlast. Fall I. Die Entwicklung In (1 -+| 1 
=t—?/2 + 8/3 —..- gestattet uns, den Reihenausdriicken (6) eine geschlossene Form zu geben 
Man erhalt in der Tat durch wiederholte Anwendung der Summenformel! 


op eee es n 7 1] ny il Co} ea pe Uh 
=e “sine + sin aa = ; In. is (8) 
a Gye (y —7) 
a b 
folgende fiir den Fall I giiltige logarithmische Ausdriicke: 
Ow Pb i| 
axe = oe /B B, P—2# [p L(a) pale L(p)] ? 
bd il 2 
(9) 
ew Pb 1 
dy” =; 4x2 B, p— ae [p L(p) —= Ot L(x)] 
mit 
Coj a oy ”) 
L(x) = In - eee Se ee (10) 
Co} ad cos — oS ”) 
a 


und einem analogen Ausdruck fiir L(/). 


Fall II. Vergleicht man die elastische Flache (4,11) mit der elastischen Flache eines isotropen 
 Plattenstreifens im selben Belastungsfall, also mit der Flache 


PP wl [UGE eae ene CAT) enn IU 

W = 75R (1 + — Je sin —— sin = ; (11) 
n=l 

_wo B die Biegesteifigkeit der isotropen Platte sein soll, so sicht man, da} man von der Flache (11) 

direkt auf die Flache (4,11) schlieBen kann, indem man in (11) « durch «/A und B durch B, d ersetzt. 

Die wohlbekannten Ausdriicke fiir die Kriimmungen des isotropen Plattenstreifens? gehen dann 


 tiber in 


Cie 127 f Wi: A x S = 

ae 8 By SI” (12) 
Bei PS) BEA ee E 
dy? 82 BA 8 B,b 


1 Siehe z. B. W. Magnus und F. Oberhettinger, Formeln und Satze fiir die speziellen Funktionen der mathe- 


_ matischen Physik, S. 214, Berlin 1948. ; 
2 Vgl. z. B. A. Nadai, Die elastischen Platten, S. 86, 95. Berlin 1925. 
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mit den Abkiirzungen 


(On}) ieee 65 atl } 
ba b 
Lea Tb % nm (y —n) 
Co} 5 ae b 
hee i eee (13) 
: en ois. | 
d = Sy a(y +n) 
uD Gof - Loe e Bor s Mee 7 ) | 


Fall III. Mit Hilfe der Substitution (5) ergeben sich aus den Gleichungen (9) und (10) die 


folgenden Ausdriicke fiir die Kriimmungen 


2 12 / f. 
ae 3b B f In +8 (1 —#2))> 
i i : (14) 
at = a In e1 ‘ sees | 3 
oy? 8b B, 4, os B’ (G1 — Pe 
worin 
2 [ x x u(y + 7)? Spe Robby 22 
a = eae aes cos —— | + Gin ae (15) 
und 
Sin oa sin F 
Ba =arc tg (16) 
Pe oj — cos P cos — aS ”) 
sein soll. 


Die Ausdriicke (9), (12) und (14) in Verbindung mit den Gleichungen (7) erlauben die Biege- 


momente der Platte in beliebiger Nahe von der Punktlast zu berechnen. Zur Berechnung der 


Scherkrafte braucht man nur die jeweiligen Ausdriicke fiir 0?w/0x? und 0?w/dy* in die Formeln 


Fw B,v,— B,r,\ Bw 
qx =—Bis—(H+ 2 ax Oy?’ (17) 
Pw B,v,— B,»,\ ®w 
BoB VAD) 8 a pe | 
poe 18 dys (H 2 ax? dy 


einzusetzen und die Differentiation auszufihren. 


4, Biegemoment im Schwerpunkt einer belasteten Flache. In den technischen Anwendungen 
der Plattentheorie sind die Einzellasten stets auf eine Flache, meist einen Kreis oder ein Rechteck 
verteilt. Sind die Abmessungen dieser Druckflache klein 
im Verhaltnis zu der kleinsten Entfernung der Last vom 
Plattenrand, aber nicht klein im Verhaltnis zu der Platten- 
dicke!, so kénnen die nunmehr endlichen Werte der Biege- 
momente im Schwerpunkt der Lastflache auf Grund der 
Ergebnisse des vorhergehenden Abschnittes unschwer be- 
rechnet werden. Die weiter abgeleiteten Momentenformeln 
gelten zwar nicht streng, aber mit einer fiir technische 
Zwecke véllig ausreichenden Genauigkeit. 

Als Beispiel wahlen wir eine gleichférmige Verteilung 


Abb, 2. Lastverteilung auf eine Kreisflache. und beginnen mit dem 


Fall I, > 1. Hine elementare im Punkt (—x,, y,) wirkende Last von der GroBe dp = Pdxdy/xc? | 


erzeugt im Mittelpunkt des Kreises (0, y) Kriimmungen der elastischen Flache, die sich nach den 
Formeln (9), (12) oder (14) leicht ausdriicken lassen. 
Da nun die Entfernungen x,, und 


z=" (18) | 


‘In diesem Fall ware die Theorie der dicken orthotropen Platte zur Berechnung der Beanspruchungen 
heranzuziehen. 


der Last P auf eine Kreisflache vom Halbmesser ¢ (Abb.2) | 
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der Voraussetzung gema8 klein sind gegen b, so diirfen wir zunachst in der Gleichung (10) 


Sopets x 
Co} x 1+ tae? ) 
7 (¥1 —) m0 3 | 
cos ; ~1 apr? (19) 
Tey, 4 2 
cos ae ~ cos = ml ley) sin’ 


setzen. Mit diesen Vereinfachungen wird 


Ef): 2. In (2 sin — ‘ In( 


, | (20) 


und analog driickt sich L() aus. Zur Vereinfachung der nachfolgenden Integration gehen wir noch 
zu Polarkoordinaten iiber und setzen 


x =reos?, z=rsin?, dxdy =dxdz =rdrd}, (21) 
wonach sich 


Ip 
ap = or dr dd (22) 


€rgibt. Multiplizieren wir nun beide Teile der Gleichung (20) mit dp und integrieren iiber die Last- 
flache, so berechnet sich hieraus der neue We,t von L(x) zu 


ia eee coon 
Lia) = 2 In (2 sin 4 cere | | in & ae £ el eo “| r dr di 
066 


ae Oe b? 


oder, nach Ausfiihrung der Integration, zu 


iMG Sonn (2 Seg (23) 


Jetzt haben wir nur noch den Ausdruck L,(x) und den analogen Ausdruck fiir L,(8) an die Stelle 
von L(x) bzw. von L(f) in die Gleichung (9) einzusetzen. Nach einigen Umformungen kommt man 
dann zu den folgenden Ausdriicken fiir die Krimmungen im Mittelpunkt der Kreisflache 

1 


Tl eo acla b b } 
22 Ph 1 Abe? sin" «In (14+ 4) —pin(1 +.) 
= — | — = —— 
0x? 27?) B, B, poe cs It € fp? — a? i 
; b b = 
Pw Pb 1 Ake Bae «In (1 ! =.) pin(1 ! =3) | 
ay 27 B, |P+« 7 ¢ a e. Baa ie 


Praktisch ist es von Wichtigkeit zu bemerken, daB das zweite Glied in den eckigen Klammern von 
der Lastenstellung, sowie der Lastkonzentration véllig unabhangig ist und sich als Funktion der 
beiden Steifigkeitszahlen / und y ein fiir allemal berechnen laft. 

Fall II, ~ =1. Hier kann man entweder in der soeben dargelegten Weise vorgehen oder aber 
die nachstehenden Ausdriicke als Grenzwerte aus (24) gewinnen. Das Ergebnis lautet 


ee) } 
aw PA ec ehe 5 I 
ae 40 B® net] Be A+ 1’ 
(25) 
4bA sin 
du? a. SEE ee! 
dy? ss A BAA eC ly Be WL 


Mit 2 = 1 gehen die vorstehenden Formeln in die wohlbekannten von A. Nadat herriihrenden Aus- 
driicke fiir den isotropen Plattenstreifen iiber. Das von A. Nadai benutzte Verfahren, den Span- 
nungszustand des Plattenstreifens um die Einzellast herum mit dem einer zentrisch belasteten 
Kreisplatte zu vergleichen stot bei der anisotropen Platte begreiflicherweise auf Schwierigkeiten?. 
Diese umgeht man am einfachsten durch direkte Integration iiber die genaherten Ausdricke fiir 


1 Im vorliegenden Fall II ist das Problem gleichwertig mit dem Problem der Berechnung der Biegemo- 
mente eines isotropen Plattenstreifens im Mittelpunkt einer elliptischen Druckflache. 
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die Kriimmungen und gelangt auf diese Weise (im Falle der isotropen Platte) genau zu dem Er- 
gebnis von Nadav. 

Fall III, ~ <1. Hier erspart man sich die nochmalige Integration durch Substitution der 
Ausdriicke (5) fiir 1/x und 1/8 in die Gleichungen (24). Die beiden Kriimmungen driicken sich nun 


folgendermaBen aus 


1 
Den HEY 
20 Pp 4b e* sin b ; 
= aes alin —« no— py |. 
ihe Pp 4 be sin =~ 
—_— «’ n—_—_~—. — o Ino + fy }, 
dy? 4b) B, B, c 
worin 
b \2 b? 
ae San eked 
O =(" +3) ye 
p= ete 5 (27) 


sein soll. Die GréBen @ und g hangen hier wiederum nur von den Zahlen A und yw, d.h. von den 
elastischen Eigenschaften der Platte ab, nicht aber von der Stellung der Last und der GréBe des 
Lastkreises. 


5. Biegemomente einer Rippe. Soll eine orthotrope Platte naherungsweise einen Tragerrost oder 
ein System von Plattenbalken ersetzen, so interessiert in erster Linie das auf eine Rippe oder einen 
Trager anfallende Moment oder auch die Querkraft. Wir haben hier zwischen den Rippen der 
Richtung x (gegenseitiger Abstand u) und solchen der Richtung y (gegenseitiger Abstand v) zu 

unterscheiden. Zur Vereinfachung rechnen wir da- 
=. bei mit einer im Schnittpunkt zweier Rippenachsen 
oe 
: Ayes 
D 


angreifenden Punktlast (Abb. 3). Wir beginnen mit 


Sie 


der Bestimmung des Mittelwertes der Kriimmungen 
0?w/dx? und 0?w/dy? fiir eine Gesamtstrecke von v 
bzw. u in einer zum Angriffspunkt (0,7) der Last 
symmetrischen Lage (Abb. 3). 


FalllI, a) Rippen in der x-Richtung. Unter 


seal der Voraussetzung, daB die im Abschnitt 4 gemachten 
ie Vereinfachungen nach wie vor statthaft sind, haben 
wir in Gleichung (20) x, = 0 anzunehmen und als- 


dann den Ausdruck 


Abb. 3. Rippenplatte mit Einzellast. La) =— L(x) dz 


zu bilden. In ahnlicher Weise bestimmt sich L,(6). Ersetzt man jetzt in den Gleichungen (9) 
L(x) und L(f) durch L(x) bzw. durch L,(f), so erhalt man die folgenden gemittelten Kriimmungs- 


werte: 


dhesne 


Pw Pb 1 b 

(a) 2 = IID Dp. hoe In ) 

a 27?/B, B, P+ % RV (28) 
ote) 

aw 36 Pb 1 4b esin-— 

dy? } 27° Bo pra” ev : 


b) Rippen in der y-Richtung. In der Gleichung (20) wird jetzt z = 0 gesetzt, die Integration 


geht nach der Formel 
u/2 
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und das Ergebnis ist 


yw Pb 4e sin ~ 
( a eet een, Pe a lap 
Ox? * 202/B, B, \P + @ u B—a /’ 
29) 
Po Pb 4esin —/ 5 
w fe el ieee b | BmB—alne 
dy? , 2 7? B, Bra u B—a /* 


Fallil, «4 =1. Die nachstehenden Ergebnisse erhalt man am einfachsten aus den Aus- 
driicken (28), (29) durch entsprechenden Grenziibergang!. 
a) Rippen in der x-Richtung: 


4besin — 


Ow\ J err l b 
(ae : 47 B, sates hie 

a UT) ey 
Se Pp 4besin = 
OVE 42 BA = Tv 

b) Rippen in der y-Richtung: 

_ NH 
aw Bon Pi : 464 sin —— 
Ox? Pad 4x B, 3 TU : 

(31) 


EE 
Ow\ P sae 
(a3), a) a ee | 4) 
Fall III, ~ <1. Macht man wieder von der Substitution (5) Gebrauch, so liefern die Glei- 


chungspaare (28) und (29) die folgenden Ausdriicke. 
a) Rippen in der x-Richtung: 


aU 
4b e sin — 


fumes Pee b 
ax} | dB oe nee” 
(32) 
HLT 
= Ba ; 4besin>— 
= n —'. 
(a5). 46) B, B, mv 
b) Rippen in der y-Richtung: 
/ 
2w 5 P a 4besin—— Age | 
ya =e —fp are t8 
‘ (33) 
02w 12 aT eed cet t “| 
9 = if a n = ee arc S ay, 0) 
(a3), 4b |B, Bs mu 24 


Die Biegemomente in den Rippen berechnen sich jetzt nach den folgenden, den Ausdriicken (7) 


entspringenden Formeln 
[/O?w Ow 
va eG 


worin fiir die Kriimmungen, dem jeweiligen Fall entsprechend die Ausdriicke (28) bis (33) einzu- 
setzen sind. Es stiinde nichts im Wege, durch entsprechende Integration iiber die Ausdriicke (17) 
analoge Formeln auch fiir die Rippen-Querkrafte zu erhalten. 


6. Der in der Querseite eingespannte Halbstreifen mit einer Einzellast. Sehr cinfache Resultate 
ergeben sich auch fiir den in Abb. 4a dargestellten Begrenzungsfall einer orthotropen Platte; die 
Rander y = 0, b sind dabei nach wie vor frei aufliegend. Wir beginnen wieder mit dem 


1 Es ist dabei bequem, die in der FuBnote 5 S, 90 erwahnten Bezeichnungen einzufiihren und zum Schlub 
@ — 0 zu setzen. 
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Fall I, «> 1 und betrachten als erstes einen vollen Plattenstreifen, der zwei gleiche aber ent- 
gegengesetzt gerichtete Punktlasten P nach Abb. 4b tragen mége. Die betreffende elastische 
Flache ergibt sich sofort durch Superposition zweier Ausdriicke der Form(4, I), wobei in dem ersten 
x durch € —x, in dem zweiten, mit der negativen Einzellast, « durch € + x zu ersetzen sind. Aus 
der hieraus resultierenden Flache 


AYE ol = vo bo BOLE ee ne | WHEY. Woy 
oe Eee Se 0 ee eer, ee 35 
= RB Pow oo 4 e ° Gin gee Sin . )sin ee as (—é<%=&) (35) 


oder auch unmittelbar aus der Abb. 4b folgt, da8 die Durchbiegungen w, und die betreffenden Biege- 
momente auf x —0 verschwinden. Fiir den so erhaltenen freigestiitzten Halbstreifen x > 0 
berechnen wir noch die Neigung der Tangente 


oe) né& n& 
Ow, BPW? Selly aac 7) 3 NAA) 0, WOE SE 
Jee pe aes o P 6 
lea oe ae > pale e sin —[— sin — (36) 


n=1 


Um diese Neigung zum Verschwinden zu bringen, beniitzen wir eine Zusatzlésung w, (x, y), die der 
homogenen Differentialgleichung (1) geniigen und auBerdem eine Durchbiegung w, = 0 auf x = 0 
er ek soll. Eine solche Lésung ist 


CO mx nx 
. . n Jt 
CUP = C,, ( B__e « ) sin ~ ; 7 sin ; ~ 6P 


—_ 


Z 
Ys 


ie 
Abb. 4. Der auf der Seite x = 0 a) eingespannte und Abb. 5. Eingespannter Halbstreifen Abb. 6. Unendliche Kragplatte 
b) freigestiitzte Halbstreifen mit Einzellast. im neuen Koordinatensystem. mit Einzellast. 
mit gewissen willkirlichen Konstanten C, fiir jedes n = 1, 2,3,... Nun ist 


Ow, - — 1 1 5 OIG) c eaESY 
a> C(3 z)n sin 5 sin (38) 


und die Bedingung 0w,/0x + @w,/dx = 0 bei x = 0 liefert erst die Konstanten C, und alsdann die 
gesuchte Zusatzlésung 


né& n& nx 
DIP ap —_ 1 BI = * a | ni? no 
= ; Mies gee : tee am 
We TABS (6 ay Pe) e er @ ie C0 SU SUL aes (39) 


nl 


die elastische Flache des in x = 0 eingespannten Halbstreifens x > 0 ergibt sich hiernach zu 
w= w, + Ws. (40) 
Wir beschranken uns weiter auf die Betrachtung des Einspannmomentes lings x = 0, das aus- 


schlieBlich durch die Teillésung (39) geliefert wird. Die erste der Formeln (7) ergibt fiir das fragliche 
Moment 


Dh Ae 4 — 1 ( se *) 3, Lae ny 
1) boy, a SOS —{e #®—e *)sin Sih 
(maa => sin (41) 
Dieser Reihenausdruck ist wiederum gleichwertig dem geschlossenen Ausdruck 


(m,)en0 = — "5g —— [L(B, €) — Los €)] (42) 
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Der Logarithmus L(x, &) ist dabei durch den Ausdruck (10) gegeben, nur mit € an Stelle von x 
und L(f, ) ist die entsprechende Funktion von 8 und €. Um das Ergebnis (42) auch bei unendlich 
groBer Streifenbreite zu verwenden, verlegen wir die x-Achse in die Streifenmitte (Abb. 5) und nehmen 
uns vor, das Moment m, ebenfalls in Streifenmitte, d.h. in x — ¥y, = 9, zu berechnen. Im Ausdruck 


(42) haben wir dann 7 durch es ++ 4, und y durch 6/2 zu ersetzen. Das Moment wird dann gleich 


p ee b 
beg (43) 
[set F —cos aa [cot == Cos ae 
3 % 


(coi = + cos ah) (cot ie yee - 


b “ b 


Bisher haben wir uns an die Hubersche Bezeichnungsweise gehalten. Die weitere Betrachtung 
wird bedeutend vereinfacht, wenn man, wie bereits eingangs erwahnt, u = oj 2 w setzt und somit 


Lees Wad 
Ce aw a Pent (44:) 
erhalt. Setzt man dies in die Gleichung (43) ein und geht dann zu b> ow iiber, so gelangt man 
schlieBlich zum folgenden Ergebnis 
Oe tes yes ea Fy 
An Ginw (€e—%)?+(A7n,)? © 


(m,)x=y,=0 = (45) 
So groB ist also das Biegemoment im Punkt x = y = 0 des eingespannten Randes einer ortho- 
tropen Platten-Halbene, falls man eine Punktlast P in « = é, y, = 7, anbringt. Das Resultat (45) 
14Bt sich auch in der dimensionslosen Form 

PA eto + )2 tg? 8 
4x Gino eo + Ptg?d 


(m,)r=y,=0 = (46) 
schreiben, sofern man 7,/£ = tg # setzt (Abb. 6). Nimmt man insbesondere 7, = 0, d. h. tg 9 = 0 
an, so wird einfach 

Pio 
a Gin * 


(m,)x=y,=0 = (47) 
Das Einspannungsmoment einer parallel zur x-Achse verlaufenden Rippe lieBe sich schlieBlich durch 
Integration des Ausdruckes (45) zwischen den Grenzen 7, = —v/2 und 7, = v/2 bestimmen. 

Fall II, ~=1. Der unter ,,Fall I* beschriebene Vorgang laBt sich fast wértlich wiederholen. 
Durch Uberlagerung zweier Lésungen der Form (4 II) ergibt sich die elastische Flache des frei- 
gestiitzten Halbstreifens in der Gestalt 


oo nag 
IP iS 1 no Le et ED eR Be epg ey oe ee 48 
“1 78 BA a ne (4 bA )gin bse br pel b pee) 
oe (—E<xS8) 
Hieraus erhalt man 
Ce a ae ao 
Ow EE Sa gp) Poin — ein (49) 
be | vez a By n b b 
1) 
Die Zusatzlésung von der Form 
co niwx ; ee 
* pores ieee nN 7 . 7 
w= * > gare? sin “> “sin — (50) 
p=. E 
liefert 
CO 
Ba ES Cisne ain (51) 
Ox } <0 Me b b 
n=l 


Die Konstanten C, bestimmen sich aus der Bedingung 0w,/dx + @w,/dx = 0 bei x = 0, Die Sub- 
stitution der Werte C,, in die Gleichung (50) fiihrt auf die elastische F'lache 


= co aa 
ye ee a ee ee etic pe pamecas ; (52) 
n b b 
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98 
die sich auch in der geschlossenen Form 
Cor =o) aes my + ) 
aferee PAxsy, bia b ae (53) 
4 4x B, Coj T(% +> Bape 


darstellen laBt. 
Das Biegemoment langs des eingespannten Randes x = 0 ergibt sich nach den Gleichungen (7) zu 
/ 


ome ore Sa on I 

(m,)x=0 rz 1 | 5a | a 26 Coj GG eos my ae ”) Sof ma & COs may ahs ”) ( ) 

‘ bh b ba b 
Mit Bezugnahme auf Abb. 5 erhalt man insbesondere in der Randmitte 
ein wE EUS 
_ 2PéE sida sri cos is os) 
(m,)x-y,-0 = b o (ends : Paap be 2) 
er ae 

(56) 


) Pi & tf I A 
(m (/X=Y = —=—— = (2 + Br) — —— x(1 42 Se ) 


sofern wieder 7,/& = tg } gesetzt wird. 
Fall lll, uw <1. Die weiter folgenden Ergebnisse erhalt man entweder ausgehend von der 
Gleichung (4 ITI) oder auch direkt aus der Lésung im Fall I durch die bereits wiederholt benutzte 


Substitution (5). 
Die elastische Flache des freigestiitzten Halbstreifens hat im vorliegenden Fall die Gestalt 
(57) | 


Pb 1 GET) 5 WR 
i) — X,, in fl y ry 
| B, Bs n> b 
worin 
2 né 
nx nx n TUX TC nTe = 
oe (gin =~ C08 —- COS 2 + Gp) — sin — sin : aoe 
a B p B p 
b né& 
i, We? nx né nx nx SW Wo fae 
i (Ot = cos —— sin — O)) =. cine cose yea 
aan ame rs ai Caey all 
sein soll. Die Zusatzlésung 
as = 1 <) 
pel Peay 
Pb pr—a” e 5 US 6 WER og A) 9 EET) - 
li, = — ; ; sin —- sin — sin sin (58) 
zn? |B, B, a . n p p b b 
ph 


erfiillt die friihere Bedingung dw,/dx + dw,/dx = 0 auf x =0. Das Einspannmoment auf dem 
Rande x = 0 ergibt sich am einfachsten durch Transformation des Ausdruckes (42) unter Zuhilfe- 


(59) 


nahme der Relationen (5). Man erhalt 
Bp, 
m= — SE —) 
worin 
g Sin oe 
Pl al ™ 8B’ 
== are tg 2 2 (60) 
, Ss S my ate 11) 
P2 oj <7 Cos B cos ee 
bezeichnen soll. 

Wir wollen jetzt vom Koordinatensystem nach Abb.5 abermals Gebrauch machen und das 
Biegemoment in der Mitte des eingespannten Randes berechnen. Mit y = 6/2 und 7 = ie + 
ergibt sich dieses Moment wieder nach Gleichung (59), doch gilt diesmal 

van Sin Ss sin “: 
= arc tg : (61) 
, ar & é Te 1), 
YP, Co} 7 coe ae cos 
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Bevor man zu b —> o0 iibergeht, empfiehlt es sich (4 = cos 2 m zu setzen. Man hat dann 


een isda! : 
~ 71 c08 w’ B fo sin © (62) 


(o 5 . 
7 Sin W 


Mit b — oo verschwindet y’, und mit dem entsprechenden Grenzwert von gy’, lautet das Endergebnis 


& sin 2 w Pi sin 2 


= are tg (63) 


PBcos2m+/2 2 A2nsinw © cos 2m -+ A? te? h’ 


ord PY 
(m)x—y.—0 ~ Sasno” '8 


worin wiederum 7,/& = tg # sein soll. Der Ausdruck (56) ergibt sich hieraus als Sonderfall bei 
wo 0. 


7. Schlu®bemerkungen. In den Abschnitten 3 bis 5 wurden in einer Reihe von Fallen geschlossene 

Ausdriicke fiir die Biegemomente des freigestiitzten Plattenstreifens und des Halbstreifens ent- 
wickelt, deren Berechnung mit Hilfe der iiblichen einfachen Reihen sich als recht zeitraubend 
erweisen miBte. Die nur von den Plattensteifigkeiten, nicht aber von der Lastenstellung und der 
Lastkonzentration abhangigen Glieder lieBen sich leicht in Form von Tabellen angeben, wodurch 
eine weitere Vereinfachung im praktischen Gebrauch der Momentenformeln zu erreichen ware. 
_ Geht man von dem Plattenstreifen zu Rechteckplatten mit Einzelkraft-Belastung iiber, so diirfte 
es klar sein, daB der singulare Teil der Lésung durch die Abanderung der Randbedingungen unbe- 
rihrt bleiben miiBte. Bildet man also die Differenz der Biegemomente oder der Plattenkrimmungen 
des Plattenstreifens und der Rechteckplatte, so wiirden die betreffenden Fourierreihen auch an der 
Angriffsstelle der Einzelkraft gut konvergieren. Mit Hilfe von Spiegelbildern lieBe sich dabei der 
Spannungszustand der freigestiitzten Rechteckplatte aus demjenigen des Plattenstreifens 
besonders einfach entwickeln. 

Die betreffenden Zusatzmomente sind in diesem Falle bekanntlich durch Theta-Funktionen 
darstellbar; praktisch geniigt es bereits, den EKinflu8 von ganz wenigen positiven und negativen 
Einzelkraften, die auf den Plattenstreifen in passender Stellung aufzubringen sind, zu beriick- 
sichtigen. 

SchlieBlich ist zu bemerken, daB Probleme der Biegung von Tragerrosten oder Rippenplatten, 
deren Systemachsen einen Winkel von 45° mit den Plattenrandern einschlieBen, sich gleichfalls auf 
die oben diskutierten Probleme zuriickfiihren lassen. 


(Eingegangen am 5. Juni 1956). 
Anschrift des Verfassers: Prof. Dr. S. Woinowsky-Krieger, Quebec 6 (Canada), 890 Ave Marguerite Bourgeoys. 
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Resonance-Behavior of Non-Linear One-Dimensional Gas Vibrations 


Analyzed by the Ritz-Galerkin Method 
By Eyvind Frederiksen 


1. Introduction. An attempt is hereby made to develop an approximate analysis of the follow- 
ing problem: The gas column within a cylindrical pipe, closed at one end, is excited to vibrations 
by a harmonically oscillating piston, located at the other end. What is the resonance response, 
when non-linear gas properties are taken into account ? 

Analytical solutions to the problem are hardly known, but excellent experiments have been made 
by E. Lettau! and others®. Their results indicate the existence of periodical solutions. For all investi- 
gated orders of resonance (fundamental resonance, first overtone etc.) the vibration response was 
found to include shock waves traveling back and forth through the pipe. The vibrations are thus of 
a highly non-linear character. 

It will here be our aim to establish an approximate closed-form solution by means of a Ritz- 
Galerkin procedure, together with energy considerations. The general form of the approximate 
solution will be chosen with close agreement to the evidence of Lettau’s experiments, and our 
analysis will therefore be restricted to a quantitative evaluation of certain solution parameters. 

It might be mentioned that conventional analysis of non-linear flow (based for instance on the 
method of characteristics) involving successive calculations on wave propagation, is hardly appli- 
cable here, since the gradual approach through such numerical calculations towards a steady-state 
solution is likely to possess poor convergency. In this respect Lettau made an interesting separate 
experiment, which shows that not less than some ten cycles (of the piston motion) are needed to 
bring the gas column from the state of rest up to the fully developed steady-state vibration. Another 
drawback of such calculations is, that they require in each case numerically specified boundary 
conditions, and therefore they lack much in generality. 

In the following section we shall sketch the mathematical background of the problem to the 
extent necessary for the Ritz-Galerkin analysis, which is given in section 3. The supplementary 
energy considerations are given in section 4. The final evaluation of the solution parameters is 
found in section 5. The results are compared with Lettau’s experiments and are discussed in general. 
Some of the analytical details have been omitted from the text, but may be found in the appendix. 

2. Introductory analysis, Evidence of Lettaus’s experiments. Let us first consider the equations 
of motion for one-dimensional unsteady flow. The dynamical equation reads 


u,+uu, 4 =O (1) 
Q 
where u(x, t) is particle velocity, p(x, t) pressure and @(x, t) density of the gas. The length variable 
is x, the time variable is t. Subscripts denote partial differentiations. The linear term r u in (1) is 
included to account for flow friction, which we shall assume acts only in the immediate vicinity of 
the pipe wall; r is a friction coefficient. The friction or damping will be considered relatively small, 
since this will appear as a necessary condition for some essential steps in our analysis, but it might 
also be emphasized, that this restriction overrules the necessity for consideration of the non-linearities, 
which the friction term actually may involve. 
The continuity equation is 


We shall assume that changes in the state of gas are adiabatic, i. e. 
Gp CP, aa 
Lae a” , ®) 
Here a =a (x, t) is defined as a (positive) gas quantity. x is the ratio of the specific heats of the gas. 


In order to preserve the assumption (3) consistent with the friction term in (1), it is necessary to 
further assume that the developed friction heat radiates through the pipe wall to the surroundings, 


1 EK. Lettau, Deutsche Kraftfahrtforschung 39 (1939). 
* E. Schmidt, Z. VDI 79 (1935) S. 671 and VDI-Sonderheft ,,Dieselmaschine“ VI (1937) S. 125; C. Mayer- 
Schuchard, Schwingungen von Luftsiulen mit groBer Amplitude, Forschungsheft 376 (1936) S. 13—22, 
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so that no part of this heat remains to heat up the gas. This is a reasonable approximation when the 
thermal conductivity of the pipe wall is much larger than that of the gas. 

Equations (1), (2) and (3) in u, p and @ are sufficient to fully describe the dynamical behavior of 
the gas column. However, much simplification is obtained, when (3) is introduced in (1) and (2) 


to give a set of only two equations, where now the dependent variables are u and a. We get respec- 
tively 


E(u, a)=u,+uu,+kaa,+ru=0, (4) 
2 E,(u, a) =a,-+ ua, + S au, =0 (5) 
where 
9 
| iS a 5 for air) . (6) 


Our problem is then to solve equations (4) and (5) for u(x, t) and a(x, t), when the following end 
_ conditions are imposed 
a0, 4) =0', (7a) 
u(l,t) = U, sinwt, (7b) 
I being the pipe length, while @ is the angular frequency and U, the velocity amplitude of the piston 
motion. Actually, in Lettau's experiments the piston effects a displacement excitation rather than 
a velocity excitation, but the latter is much simpler to deal with in the present analysis, and the 
error involved is small, when the displacement amplitude U;/@ of the piston motion is small com- 


pared with the pipe length. 
To complete the boundary conditions we shall further make the reasonable assumption that the 


gas motion is periodic with period T = 2z/@ at any position x. Focusing our attention to the 
time interval 0 <t < T, we get the additional conditions 


ule, 0) == ule, T),, (8a) 
ax, 0) = a(x, T). (8b) 


Before we proceed with the non-linear problem, it is here worthwhile (for the sake of comparison) 
to briefly consider its linear counterpart. Linearizing equations (4) and (5) we get 


u,tkaa,t+ru=0, (9) 
a, aa, =a (10) 


where a, is an average of a in both x and t. (The corresponding average uy of u must be zero, since 
the total mass of the gas colum remains constant). 
In the case of fundamental resonance (« l/ay =) we have for sufficient small damping (r < w) 
with good approximation the steady-state solution 


Ue — U, fi sing x cos wt + (+) cos q x sin w ] h (11) 
Oe ep = Uy 4 = fa cosq xsin wt + (7) sin q x Cos W ] (12) 
| where 
(od) A 
=—S—- es il 
q dy 2 i) 


p=s~ <4, (14) 


q being the so-called wave number and / a damping coefficient. Note that other vibration com- 
ponents with frequencies different from @ cannot exist in this case. If they were somehow intro- 
duced, they would die out because of the damping. Condition (14) implies that bl<1, and it is 
thus seen that the terms in (11) and (12) having the factor (x/l) are quite insignificant, although, 
| strictly speaking, they are needed in order to have the end condition (7b) satisfied. 

With respect to the non-linear problem we shall, in this section, resort to a mere interpretation of 
| Lettau’s experiments, and we shall limit ourselves to consider the case of fundamental resonance. 
| Circumstances for higher order resonances are analogous, and they will be considered in the final 
4 analytical expressions only (in section 5). 
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The experiments are concerned with recordings of the velocity u and pressure pata number of 
stations equally spaced along the pipe. We have here chosen to analyze variations in a rather than 
in p. The relation between these quantities follows from (3) 


hie nk Falty rae for air) (15) 
Pp a a 


The relative pressure variations are found to be relatively small (< 0.16). We can, therefore, 
conclude from (15) that Lettau’s pressure diagrams may be interpreted without much error as 


recordings of variations in a, except for 


ss 


1) a is essentially an odd function 
of t for all x, and an even function of 


odd function of x for all t. 


cally independent of x. 


corresponds closely to a free and un- 
damped linear vibration in a cylindrical 


pipe of the same length, being rigidly 
closed at both ends. (Lettau proved 


on back and forth running waves). 
These observations enable us to ex- 
Saul press the response by an approximation 
of the general form 


0 u(x, t) = — (7) U,cosqxsinwt 


—k>A,sinngxcosnwt, (16) 


n=1 
6 facut) a = (7) Ursing xeos ot 
nm Fe iF 0 k\1 l 
’ 3 = . 
+> A,cosngxsinnwt (17) 
WZ im; 
Fig. 1. Recordings of pressure p and velocity u from equally spaced stations where n = Il. 2 Ohewen (The “tilda” . 


along a cylindrical pipe with piston excitation at one pipe end. Reproduction 


from E. Lettau, see ref. (1) on page 100. sign is here used to emphasize that we 


are dealing with an approximation). 
The insignificant terms with factor (x/l) are merely taken over from (11) and (12) just to have the 


end condition (7b) satisfied. The vibration is seen to contain infinitely many components each | 


having a frequency which is an integral multiple of the fundamental frequency w. 
However, probably the most striking feature of fig. 1 remains to be pointed out: 


4) The vibrations contain discontinuities in both a and u propagating back and forth through | 


the pipe as a shock wave. 


The shock wave pattern is apparently connected with the proper cooperation of the infinitely 
many vibration components and with the maximum amplitudes A, of these components. From (17) 
is seen that the As are merely the Fourier coefficients of the a-variation at x = 0. Inspection of 
fig. 1 shows that this periodical variation may be approximated by 


a(0, t) — ay = — A* sing (w t — 2) for Va tral (18) 


where A* is a maximum amplitude and « a ”form factor‘‘, which appears to be smaller than unity. 


at 

115 - a constant factor. 

110 Fig. 1 shows a reproduction of the | 
105 diagrams for fundamental resonance. 
cade From these we may summarize the _ 
aa | experimental evidence as follows: | 


2) a) (as a time average) is practi- | 


3) the recorded resonance response | 


085 
- x for all t. —u is essentially an even | 
function of t for all x (except for a slight | 
disturbance at the piston end), and an | 


this through a graphical, linear analysis 
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Fourier analysis of (18) yields 


MAT for n=1,2,3... (19) 


a (n?—@?) 
Incidentally, the discontinuities are associated with the circumstance that A, approaches pro- 
portionality with i/n for increasing values of n. 
With the various statements above we are sufficiently prepared for the actual analysis. 


3. Ritz-Galerkin analysis of the problem. So far, we have proposed to consider (16) and (17) as 
an approximate solution, and for the coefficients 4, we have proposed the relation (19), leaving the 
parameters A* and « to be determined. 

For the present analysis it is appropriate to start out from somewhat less specified relations. 
We shall write the approximate solution in the more general form 


u(x, t) = ox US PrlXz t) 2 (20) 
n=0 
a(x, t) aay ae Ss A, PnlXs t) (21) 
n=0 


where the “coordinate functions” gy, and y, are so chosen, that i and @ are periodic (with time 
period T) and satisfy the end conditions (7). In order to establish conditions determining the un- 
known coefficients U,, and A,, we shall require that the expressions FE, (i, @) and E, (u, @) — as 
found by substitution of (20) and (21) into (4) and (5) and thereby considered as functions of x 
and t — are orthogonal to any of the functions , and y,. Taking the periodicity into account we 
thus introduce the conditions 


iid B 

J [ Elis @ g(x, 0) dt dx = 0, (22) 
ie 

J [FG @ pale t) de dx = 0 (23) 


where i = 1, 2 andn = 0, 1, 2,3... These equations represent non-linear, algebraic relations be- 
tween the set of coefficients U,, and A,,. 

We shall not try to justify conditions (22) and (23) by a complete analysis originating in the 
calculus of variations. Reference is made to the ideas of W. Ritz! and B. G. Galerkin® and especially 
to the application of those ideas to non-linear vibration problems for systems of finite-degrees-of- 
freedom, as derived by K. Klotter?. For such systems corresponding conditions have been fully jus- 
tified. Equations (22) and (23) are thought to represent merely a formal extension hereof to 
cover our particular continuous system problem, and it is in this sense that the present treatment 
has been designated a “Ritz-Galerkin analysis’. The plausibility of the assumed conditions may be 
supported by the following remarks: 

1) If (20) and (21) should accidentally represent an exact solution (u, a), then the validity of (22) 
and (23) is evident, since in that case E, and E, vanish identically in x and t. 

2) On the other hand, if the two groups of functions :, and y,, each form a complete set of ortho- 
gonal functions, then conditions (22) and (23) express that E,(u, @) and E(u, @) are orthogonal to 
any function in x andt. Hence E,(u, 4) and E,(u, 4) must vanish for all x and t, i. e. (u, @) becomes 
identical with the exact solution (u, @). 

Thus if (22) and (23) cannot be strictly justified as a “direct method” for solving the underlying 


) variational problem, they at least may be considered as expressions for “orthogonalizing” * the error, 


e. g. functions E,(i, @) and E(u, 4), with respect to the coordinate functions @,, and »,. 
Now, in the present problem we have experimental evidence that the significant terms in an 
approximate solutions of form (20) and (21) contain the following coordinate functions [see (16) 


and (17)| 
y, (x, t) =sinngxcosnat, (24) 


y, (x, t) = cosngxsinnwt (25) 
1 W. Ritz, Crelles J. f. Reine und Angewandte Mathematik 135 (1908) S. 1. 
2 B. G. Galerkin, Vestnik Inghenerov | (1915) S. 879 (in Russian). 
3 K. Klotter, Non-linear vibration problems treated by the averaging method of W. Ritz, Proc. Ist Nat. 
Congr. Appl. Mech. (Chicago 1951) S. 125; K. Klotter, Trans. I. R. E. (1954) 8.13. 
4 For the “method of orthogonalization”, see L. Collatz, Numerische Behandlung von Differential- 
gleichungen, 2. Aufl., S. 30, Berlin 1955. 
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where n = 1, 2,3... These represent of course only a few of the functions needed to make up a 
complete set. But it is reasonable to assume that when we restrict ourselves to these functions, the 
effect is only to exclude those of the conditions (22) and (23) which are insignificant, and to slightly 
change the remaining conditions. The latter thus become an appropriate tool for approximate 
determination of the relevant coefficients U,, and A,,. 

Through the choice of the coordinate functions in (24) and (25) we have disposed of the first terms 
in (16) and (17), containing the factor (x/1), which terms correspond to Uy pq and Ay Yo in (20) and (21). 
This neglect is not necessary, but it effects without much influence a substantial simplification of 
our analysis. 

From expressions (20) through (25) we get the following four sets of orthogonality or Galerkin 
conditions (with n = 1, 2,3...): 


er 
i E(u, @ sinng xcosnwtdtdx =0, (26) 
36 

tan 

| [ E,(@, @ cosngxsinntdtdx=0, (27) 
00 

seh 
[ [ E\@ @ cosng«sinn aw tdtdx =0, (28) 
0 6 

Lv 

pi E,(u, @) sinn gx cosn mw tdtdx =0. (29) 
06 


In order to simplify the following we shall anticipate that each of the conditions (28) and (29) 
will be satisfied for all values of n, when 


U, =—kA, (30) | 


which is in close agreement with (16) and (17). We can thus continue with the simpler approxima- 
tion 


u(x,t) = —k > A snngxcosnwt, (31) 
n=1 
a(x, t) =a, + SA, cosngaxsinnot. | (32) 
n=1 
These expressions imply that 
u, ++ k a a, Se eee) (33) 
= 1 is 
a, ar k ao u, — ? (34) 
and we get from (4) and (5) | 
E\(u, @) =uu,+k(a—a)a,+ru, (35) 
ae ae Wcities S 
E(u, a) = wa, + os (a — a)) u, (36) 


which, with (31) and (32), become 


eo 
E,(u, 4) =e q>) 4,8 A, [sin (r+ s)qxcos(r+s)mt 


) 
+ sin (r — s) qx cos (r — s) wt} 
k(k—1 
— <==) 574, ¢ A, [sin (r + s)qx cos (r—s)awt 


+ sin (r — s) qx cos (r + s) w | 
—kr S14, sinng xcosn wt, | 


j 
! 
| 
q 
| 


f 
| 
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may k+1 
E(u, a) = — - q >) 4,5 A, [eos (r +s) gxsin (r +s) wt 


) 
| 
+ cos (r — s) qx sin (r — s) wt] | 
| 


k—] 
“4 D/A, 8 A, [eos (7 + 8) qxsin (rs) at 


+ cos (r—s) qx sin (r + s) wt] | 


where the summations are to be performed for all positive integers r, s and n. 
Substituting (37) and (38) into the remaining conditions (26) and (27) we get respectively 


De erate (39 
Sen (40) 


nm 


yy pa Mes Ae >. A.sA,— >; AsA., (41) 


n r+s=n r—s=n s—r=n 


(38) 


The two sets of equations (39) and (40) are obviously somewhat contradictory!, unless the fric- 
tion coefficient r vanishes. However, remembering that these equations are themselves the result 
of approximations, and considering only the case of small damping, it should not be too erroneous to 
continue with only one set of equations , 

Bap p 
2S k+ iq 7 (42) 
representing the “‘average”’ of (39) and (40). 
Elaborating the summations in (41) we get 


De lan Ans ee Ae SAGA Gv) An_4 A,| 
‘ + [4; 4,41 +2 A, Ani, +3 43 4,43 ++:°] (43) 
|r al) Ay Ane 2) As Ane, + (mn 3) A, Ane es]; 


Contracting the series expansions, substituting into (42) and remembering that r = 2 a, B 
and q = a/l we get 


4 pl — ~ 

= fet wi == >, 4,4,5—9 SAA (44) 
n p=1 p=1 

where p attains only positive integral values. 

To solve this infinite set of equations in the coefficients 4, appears to be a difficult task. At 
this point, therefore, it is suitable to introduce the experimental evidence we have available on these 
coefficients. Substituting (19) in (44) we get, after some intermediate steps (which are given in the 
appendix), 


2 (n? — a?) 


See siete) ee (45) 
where 
fe) =—1+420 >=, (46) 
ps ah 


1 This contradiction apparently results from the fact that in (26) through (29) we have twice as many 
equations as unknown coefficients U;, and A,. We could have avoided this redundancy, if instead of (20) 
and (21) had chosen the approximation 


co [o.e) 
u(x, t) = ~. (Wes Pn -- 2 Vin Wns 
n= n= 


lo-e} [o.e} 

a(x, t) = a De An Yn + » Bn Pn 
n=0 n=0 

where a new set of coefficients V,, and B,, but no new coordinate functions, are included. However, according 

to the experimental evidence, these new coefficients must be insignificant, and their corresponding solution 

terms may therefore be ignored without considerable error. 


8* 
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The function f(x) has been calculated from a series expansion (as derived in the appendix) and is 


plotted in fig. 2. pate | 
In order to have a unique relation between the parameters A* and q (i. e. independent of inte- 


gers n) we shall approximate the set of equations (45) by their “average”, which will be 


vy : 
a 1 == A * ‘ OXI > AT 
flos)tolee) eo pap Od a7) 
A ES considering that 
P 
F 1 (i a7) |e 
Jim E oe | 1. (48) 
8+ 4¥0r OF 
Equation (47) represents only one 
relation between the solution para- 
| meters, so we have not obtained full 
it tees information from our analysis. This 
is not the fault of the Ritz-Galerkin 
method as such, but is merely due 
aA to the fact that we have approxi- 
mated the whole set of equations (26) 
through (29) by a single expression 
A an Of (47). 
4, Energy considerations. Another 
relation between A* and « may be 
Ble AL ; o found by equating the vibration work 
95 06 a7 G8 Ga Hy exerted by the piston on the gas co- 
Fig. 2. Functions f(x), g(«) and h(«) as determined by expressions (46), (56) and (57). lumn to the energy dissipated through 
, friction. 
The mechanical work W>, transferred from the piston to the gas during one time period T is 
T 
Wp =—J pll,t) U,sinwt dt. (49) 
0 
As we consider pressure changes essentially proportional to changes in a, i. e. from (15) and (3) 
PCL, t) — Po = Po a 
ak a (50) 
we may rewrite (49) as 
7. 
Wp = —k oy a U, | [a(1, t) — ag] sin w t de. (51) 
6 
Substituting a/(I, t) from (32) we find 
k 
Wp = “9 0 % TU,A, (52) 


where A, is the amplitude of the fundamental component of (@ — a). 
The dissipation work W, per time period may be written as 


ier 
Wy =r fourdtdx. (53) 
00 
Introducing herein (31) and (14) and neglecting density changes we get 
ke ee 
Wy = 50% TBI > AD. (54) 
n=1 : 
Equating W, and W>, and introducing (19) we obtain 
Page, (UE > ; 
i FI = (1 — a?) g(a) A* sinaaz (55) 
where 
Tay ee 
ate) = Di gear) (56) 


The function g (x) has been calculated from a series expansion (as derived in the appendix) and is 
plotted in fig. 2. Expression (55) is the desired second relation between A* anda. 
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5. Discussion ot the solution parameters 4* and x. From (47) and 


dimensionless relations! : (55) we get the following two 


1 k a a 
een = al a5 = h(x) , (57) 
4n U * 

k (k 7 1 = = (1 —o?) f(x) g(a) ball sin? x 77 , (58) 


_ from which it is seen that the form factor « and the maximum amplitude A*/a, of a/ay are both 
functions of the relative excitation amplitude U,/a, and the damping measure f 1, i. e. 


=F (Bd, (59) 
= (co °F): (60) 


So far, we have considered the case of fundamental resonance only. Applying the same analysis 


to higher order resonances of order s (with s = 1 for fundamental resonance, s = 2 for first overtone 
etc.) we find 


U, Bl Ale 
oa [Pe Bae Zo 
: (; ay > hea a iE ye Vi 
RS QD; . f 

A* Uy pl make S © 6 

eee ele @ |S WW So | oh(t 

ay Sidi os es Jol | Ve / Ye, WA 

Va : = 

where the functions F and G are the same as ’ Los 


A 


in (59) and (60). The arguments U,/a, and f / / V 
3 5 oz 
003 oft Th 
5 LA 


have merely been replaced by U;,/s ay and f l/s. 

From relations (57) through (62) a diagram 
has been calculated (as presented in fig. 3), 
showing A*/a,) as function of U,/s a, for x and 
{las parameters. 

Let us here first compare this diagram 
with Lettau’s four experiments for resonance 
of orders s = 1, 2, 3 and 4. 001 

In the first three columns of Tablel are 
given the experimental data on the resonance 1a ¥ 
frequency f = /2a cycles per second, the I/: Siig 
maximum pressure amplitude P* atm. (corre- 0 gar G02 G03 Go4 005 


AME 


Vhae 
ie 


sponding to A*) and the form factor«a. The Fig. 3. Response plane, as determined by expressions (57) and (58), 
f - a KG proportional té the showing vibration amplitude A*/a, versus excitation amplitude Uj/s ay 
requencies are seen 

resonance order s. They thus determine uni- 
quely the (actual) sound velocity ay = 21 f/s 


= 344 m/s, considering that the pipe length 


with damping f l/s and form factor « as parameters. The points marked 
1, 2, 3 and 4 correspond to Lettau’s experiments. They are all located in 
the shock region above the parabola « = 1. Point la below the para- 
bola « = 1 corresponds approximately to the shock-free condition 
which — according to the present analysis — should exist for funda- 
mental resonance, if in Lettau’s experimental set-up the pipe had been 


twice as long, i. e. 24 m. 


is ] = 12,00 m. P* and « have found been from 
the pressure recordings obtained at the closed 
pipe end. Their values in Table 1 appear to be about 2 to 5 percent uncertain. The next two co- 
lumns give values for 


where R = 0.034 m is the displacement amplitude in the piston motion, and 
Ae 1 Pe iP 
dy xk po 7 


With these values are four “response points’ plotted in fig. 3, each marked with their corre- 
sponding resonance order s. From this diagram are then estimated the corresponding values of 


‘« and f las given in columns 6 and 7. 
The «-values obtained in this manner are seen to be in excellent agreement with the measured 


values. 


1 See fig. 2 for a plot of the function h(z). 
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Table 1. 


Measured quantities Calculated quantities 


i Pe 


atm, 


U A* é Bl Bol 
c. p.s 


1 14,35 0,12 0,75 0,0089 0,0170 0,75 0,086 0,040 
2 28,7 0,14 0,71 0,0089 0,0200 0,70 0,153 0,059 
3 43,1 0,15 0,68 0,0089 0,0170 0,67 0,213 0,073 
4 57,5 0,16 0,65 0,0089 0,0230 0,64 0,254 0,084. 


The corresponding f I-values are in reasonable agreement with earlier experiments made by 
K. O. Lehmann! and E. Frederiksen?. Both have found that the damping coefficient /, in case of 
(nearly) harmonic vibrations of air columns in cylindrical pipes, is 2 to 3 times as large as calculated 
from Kirchhoff’s® formula 


a2 Joy 

Boece a 

considering maximal amplitudes of the same magnitude as in Lettau’s experiments. When in (63) 

we substitute y = 15+ 10-® m?/s for the kinematic viscosity, d = 0,070 m for the inner diameter 
of the test pipe and a) = 344 m/s for the sound velocity, we find with the pipe length 1 = 12 m 


Bot = 0,0042 Yo . 


Substituting herein the values of w = 2 2 f we get fy 1 as given in Table 1. 
The agreement convinces us of the dependability of our analysis and 
justifies further investigations into the general character of the problem. 
The validity of our analysis is restricted to the interval 0,9 [a0 <1. 
The lower limit follows from (57) which shows that h(x) and f(x) cannot 
attain negative values, and according to fig.2 this requires « > 0,5. The 
upper limit follows from the fact that « > 1 would be associated with 
negative pressure jumps which are impossible for entropy reasons. Fig. 4 
shows the general character of a(0, t) — ay for selected values of x. 

Obviously « = 0,5 corresponds to the case of undamped vibrations 
(8 1=0). Here, according to (58), 4* becomes unlimited unless U, vani- 
shes. If U, = 0, we have undamped free vibrations in a cylindrical pipe 
closed at both ends. In that case we need proper initial conditions for the 
determination of A*, 

When « = 1 our approximation (31) and (32) with the restraint (18) 
or (19) becomes identical with a linear response, e.g. a pure harmonic 
vibration (without discontinuities). To check this we find — through a 
limiting process with the aid of the series expansions for f(x) and g(a) as 
given in the appendix, and again considering higher resonance orders by 


ry replacing U,/ay by U|/sao, resp. BI by Bl/s — that x = 1 specifies (57) 


(63) 


2(0t)-ay 


BX 


Fig. 4. Vibration response 
(a(0, t) —a,) at the closed pipe 
end, according to the assumption 
(18). The diagrams show the 
effect of the form factor, The 
top diagram (a = 0,5) corre- 
sponds to the case of undamped 
(free) vibration. The next di- 
agram (x = 0,75) shows a case, 
where the damping is not suffi- 
cient for eliminating the presence 
of shocks, such as in the bottom 

diagram (« = 1), where the 
response is essentially linear, 


resp. (58) to A 
2 
: (Sai) Te ie wee 
ee 1 A*2 

k(k +1) sa, =a(7| : 
and combining these relations we get 
A* ib oe Wh i Gh 
als - 


i) 


> -k Bl sa, baie 


(64) 


(65) 


(66) 


which is in agreement with the amplitude of the significant term in expression (12) for a(x, #) in the 


linearized case. 


In the “response plane” fig. 3 the “response points” (U,/sa), A*/a)) are located on parabolas, 


when a is assigned to constant values (see (58)). The interval 0,5 < « <1 is then characterized by 


Ve K. 0. Lehmann, Ann. Physik 21 (1934) S. 533. 


* E. Frederiksen, Luftpulsationer i Rorsystemer (Pulsating air flow in pipe systems), Diss. Techn. Univ. 
Copenhagen 1954. Chapter 4 (in Danish). nee oe ee 
3 G. Kirchhoff, Pogg. Ann. 134 (1868) S. 177. 
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the “response region” bounded by the ordinate axis and the parabola ~« = 1, and our analysis is 
valid for this region only. To account for response points below the parabola « = 1 (and above the 
abscissa axis) it is — with reference to (66) — reasonable to adopt the linear solution (11) and (12). 
Therefore — considering only the significant terms without the factor (x/l) — fig. 3 is, in this 
region, provided with straight lines radiating from the origin, corresponding to constant values of 
the parameter /1/s. Incidentally, this simplifie dlinear response could still be considered a solution of 
form (31) and (32) with restraint (18) or (19) in so far as this case corresponds to vanishing coeffi- 


- cients A, for n> 2, while 4, = A* anda = 1. 


Within this interpretation of our analysis we can make the following remarkable statement: 
Response points above the parabolaa = 1 correspond to essentially non-linear 
vibrations including the presence of shocks. 


Response points below the parabola « =1 correspond to essentially linear 
vibrations where shocks are impossible. 


With the aid of (64) and (65) we can elaborate this statement as follows: In a pipe with given 
“damping measure” I, shocks will be present 
only when 


U; 1 2k /pi2 0,05 
ead eo) 
and ¢ tl 
n onsequen Ay 004 
oi ee 
s et (68) 


ee wtk+l1s° 


Obviously, when the excitation amplitude 9% 
U,/a) is given, shocks can be present only 
when f 1 is less than a critical value specified 
by (67). 

Unfortunately, none of Lettau’s closed-pipe 
experiments exhibit absence of shocks. Ac- 
cording to our analysis a shock-free response 
could have been obtained if the pipe had been, 
for instance, twice as long, whereby f 1 would 
have been approximately doubled, considering 
the case of fundamental resonance. This would 0 ai Due We WE Wa 
correspond approximately to the rScep sotdehaahe Fig. 5. Response plane for the shock strength [@/a,] as determined 
point la in fig. 3. However, Lettau has made from (69), (57) and (58). 
another series of experiments, where the pipe 
at x = 0 is changed into an open end. Although the resonance responses in these experiments 
are far more complicated (presumably because of imperfect wave reflections at the open pipe 
end), they do disclose the general tendency that the pipe damping counteracts the shock for- 
mation. Thus no shocks appeared in the case of fundamental resonance, but were present at 
higher order resonances having magnitudes increasing with the resonance order. 

As fig. 1 indicates, two pressure jumps occur within the time period at any position x, and it is 
easy to show that the magnitude or strength of these jumps is maintained constant throughout the 
pipe. At the pipe ends, only one shock of double strength appears within the period. In terms of 
the quantity a we find from (18) the shock strength 


[a] = A* sina x (69) 


Fig. 5 shows a diagram (analogous to fig. 3) for the evaluation of [a]. 


G02 


607 


6. Summary. The present paper concerns an approximate analysis of damped resonance vibrations 
in a closed cylindrical pipe with harmonic excitation at one pipe end. Earlier experiments made 
by Lettau have disclosed that such vibrations may be highly non-linear and may even include shock 
waves travelling back and forth. The analysis is based essentially on the socalled Ritz-Galerkin 
method. Emphasis is put on the mechanism of the analysis rather than on a strict mathematical 
justification hereof. The theoretical results are in close agreement with the experimental evidence. 
Further a criterion for the presence and absence of shocks is derived. 

The present investigation was undertaken while the author, on leave from the Technical 
University of Denmark, was working as a Research Associate ar Stanford University in California 
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ona Research Fellowship from International Cooperation Administration, Washington D. C. Thanks 
are due to all these institutions. 

In particular the author wishes to express his sincere thanks to Professor K. Klotter, Division of 
Engineering Mechanics, Stanford University, for his kind encouragement of the work and valuable 
criticism. 


Appendix 1. Evaluation of >’. 


From (44) we have 


n—l 


D1 = 7 Oy Ae Ane” 2) Ay Ane (70) 
= jaz 


nm 


where the coefficients 4, from (19) are to be substituted. With 


K = Ax Sine (71) 
to 
we can write 
1 1 1 1 
Ae 2 
ApAL, 5 ee errs Sere (72) 
il il 1 1 
= 2 
Bp Give a lores pal eee | eel (72) 


Substituting (72) and (73) in (70) and performing the summations we shall need the following sub- 
summations 


n—l n—l n—l 
I 1 il 1 v i 
— | a 
210+ 50—P FH Goro. ! SEE? GE ee (74) 
n—l l 1 n—l 1 1 9 n—l l 
= | = aan ©, 
peer) Geel oe res) G—IH 2 p= ‘ (75 
~ 1 Tey 1 
— | 
ace) leer el 
n—lL 
cid: 1 1 if 1 
Sas: eee == jo a> ee) (76) 
a ip eee 1 1 
n Pre Or eee eels 
n—l n— 
ieee lie iro Es 1 1 } 
S| AUD eA PH) Oe = Pa : | 
n—l 
al 1 1 il 1 
4 =a. | 
Pe Baa aaB oe rye | ‘a 
n—1 nr 
a 2 1 1 1 il 
Speen a ess wenee | 
S = =i >(3 | 1 
ot (p+a)(n+p+a) i pte | nd (pata) (78) 


7 zt f i il 1 1 e il 
Paes n eer le ae (79) 
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> 1 a eel > hea eg Meal | 
pepe) Ps) eee meer 
= — 1 1 1 
Siory ae nt+pta eee) tay) 
< l 
co ea aH poles Tone 
OS 1 a 1| = 1 | 
—(p—a)(n+pt+e2) (n+ 22) p—¢ ae) | 
ea esi 
1 i wl 1 l 1 
Seat peat= n-+p—a eres | (81) 


1 4 1 1 = 1 1 


With these sub-summations, straightforward computations lead to 


Dee aa | L+20t a): (82) 


p=1l 
and substitution of (82), (71) and (19) into (44) yields immediately (45) and (46). 


Appendix 2. Series expansions of functions f(x) and g(x). Making first the following trans- 
formations of (46) and (56): 


fo-—ss2eSigtg iss +2 >| > (5) (83) 
0-2 letel* 2) yl Sige oe ae | 
P 


“Ziel Gy Gl +46) + eZ Gl 


pa 


then introducing the Bernoulli numbers B,, defined by 


ces ny! = ‘L\en 85 
B,, =2 aay (5) ( ) 
we get 
eee eV eeh geo, (86) 
n=1 tc n) 
gla) 2 >) nh Bot (87) 
fale 


These series converge rather slowly, but noting that 


ling ( We 2, <2 (88) 


a eepe mPa 
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we can considerably increase the rate of convergency as follows. Adding and substracting in (86) 
the series 


a 2 O77 = i (89) 
and in (87) the series 


CO r 4 
2 Da? n Xn” n—-2 — Gd — 2p (90) 
iat 


we get 


{== ay en) "B,—2] a2", (91) 
n=1 


G9) Gare ye 2 ye nen aes — 2latn2. (92) 
Substituting herein Bernoulli numbers 
we find for the first important terms of (91) and (92) 
5) =a hes | 1,28987 a2 + 0,16465 a4 + 0,03469 v8 on 
+ 0,00815 x8 + 0,00199 x19 + ..., 


mee 4 2 4 6 
(0) = q—gp t 2.57974 + 0,65860 a? + 0,20814 a4 + 0,06520 o |! an 


+ 0,01990 a§ + 0,00600 «1° + . 
Equations (93) and (94) specify f(«) and g(x) in the interval 0 < « < 0,5 to the accuracy of 6 ciphers. 
To obtain the same accuracy in the interval 0,5 < « < 1 would require more terms than considered 


above, especially in the neighborhood of ~« = 1. However, from (46) and (56) one can derive the 
formulas 


fl—vj=— 2-2 Gy, (95) 
g(1—a) =e) +0 =9.. (96) 


Thus the values of f(x) and g(x) in the interval 0,5 <« < 1 may readily be calculated from 
known function values for the interval 0 < « < 0,5. 


Fig. 2 shows a plot of the calculated functions f(x) and g(x) together with a plot of h(x) defined 
by (57). 


(Eingegangen am 20. Juni 1956.) 
Anschrift des Verfassers: Dr. Eyvind Frederiksen, Gentofte (Danemark) Lyngbyvej 315 A. 
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Knickung von Schraubenfedern unter Druck und konservativer Torsion 


Von H. Maag 


1. Einleitung. Die Knicklasten einer durch eine axiale Druckkraft beanspruchten Schraubenfeder 
wurden erstmals von E. Hurlbrink! und in einem allgemeineren Rahmen 1924 von R. Grammel? 
berechnet, indem die Feder in geeigneter Weise durch einen weichen Stab ersetzt wurde. C. B. Bie- 
zeno und J. J. Koch® zeigten, da die Resultate durch Beriicksichtigung des Kinflusses der Querkraft 
speziell fiir nichtschlanke Federn erheblich verbessert werden kénnen, und J. A. Haringx* * machte 
neuerdings darauf aufmerksam, daf8 man dann auch die von der Querkraft erzeugte Schiebung 
beriicksichtigen mu8. In dieser Arbeit berechnete Haringx die kritischen Lasten fiir eine Belastung, 
die sich aus einer axialen Druckkraft und einem axialen Torsionsmoment zusammensetzt. 


Diese Resultate wurden alle mit Hilfe des statischen Stabilitatskriteriums gefunden, wonach die 
Stabilitatsgrenze durch das erstmalige Auftreten von nichttrivialen, d.h. ausgebogenen Gleich- 
gewichtslagen bedingt ist. Nun hat aber H. Ziegler® gezeigt, daB® das statische Stabilitatskriterium 
nur fiir konservative Systeme legitim ist, also insbesondere auf Systeme, in denen axiale Torsions- 
momente auftreten, nicht angewendet werden darf, da diese ja im allgemeinen nicht konservativ 
sind. Die Resultate von Haringx hangen damit 
teilweise in der Luft und bediirften einer Uber- 
priifung durch das kinetische Stabilitatskriterium 
(Untersuchung kleiner Schwingungen um die 
Gleichgewichtslage), das in allen Fallen gii]tig ist. 

Anderseits stellte H. Ziegler’, indem er auf 
die den Momentvektor erzeugenden Einzelkrafte 
zuriickging, fest, da} die durch die einfachsten SSN 
Kraftverteilungen erzeugten Momentvektoren nur 
bei Einspannungen axial bleiben, bei gelenkig 
gelagerten oder freien Enden aber ihre Richtung 
andern. Die so je nach dem Kraftangriff ents I. die beidseitig horizontal unverschieblich eingespannte Feder, 
stehenden semi-, quasi- und pseudotangentialen _ IL. die beidseitig eingespannte Feder mit gegeneinander horizontal 
Momente sind alle konservativ. Fir dieBelastung 44; eater sertsbs oe cares ie 

to} . die elnseitig eingespannte eder, 
mit solchen Momenten darf also wieder das sta-  !V: Siete Spice eine o- paiiatone oye cde AE acenie 
tische oder auch das energetische Stabilitats-  V. die beidseitig gelenkig gelagerte Feder. 
kriterium angewendet werden, wonach die Stabili- 
tatsgrenze erreicht ist, wenn die potentielle Energie desSystemes erstmals nicht mehr positiv definit ist. 


1 =; 


AN 


wf IT IT IV eens 


Abb. 1. Die untersuchten Lagerungen. 


In der vorliegenden Arbeit wird nun die Knickung einer Feder untersucht, deren Belastung sich 
aus einer axialen Druckkraft und einem solchen konservativen Torsionsmoment zusammensetzt. 
Im Abschnitt 2 werden die Differentialgleichungen der elastischen Linie der Feder bzw. des elastisch- 
aquivalenten Ersatzstabes hergeleitet. Unter Beriicksichtigung der Randbedingungen werden dann 
in Abschnitt 3 die Knickgleichungen aufgestellt und zwar fir fiinf Lagerungsarten (Abb. 1). 


Selbstverstandlich miissen in allen fiinf Fallen die Lager so ausgebildet sein, daB sich die Feder- 
enden gegeneinander um die urspriingliche Federachse verdrehen und vertikal verschieben kénnen. 
Im Abschnitt 4 schlieBlich werden die Knickgleichungen diskutiert und in den einfacheren Fallen 
die kritischen Dynamen in Diagrammen dargestellt. 


2. Grundlagen. a) Der Ersatzstab. Eine Schraubenfeder, deren Lange im ungespannten Zu- 
stand J, und deren Wicklungsradius R betragt, sei durch die Belastung auf die Lange / zusammen- 
gedriickt. Die Feder sei so gewickelt, daB die Haupttragheitsachsen des Drahtquerschnittes parallel 


1 F. Hurlbrink, Z. VDI, 54 (1910) S. 133 u. 181. 

2 R. Grammel, Z. angew. Math. Mech, 4 (1924) S. 384. 

3 C. B. Biezeno u. J. J. Koch, Z. angew. Math. Mech. 5 (1925) S. 379. 
4 J. A. Haringx, Proc. Nederl Akad. Wet. 45 (1942) S. 533 u, 650. 

5 J. A. Haringx, Philips Res. Rep. 3 (1948) S. 401 u. 4 (1949) S. 49. 

6 H. Ziegler, Z. angew. Math. Phys. 4 (1953) S. 89. 

7 H. Ziegler, Z. angew. Math. Phys. 3 (1952) S. 96. 
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und senkrecht zur Schraubenachse gerichtet sind; die Haupttragheitsmomente beziiglich der Nor- 
malen bzw. Parallelen zur Schraubenachse durch den Schwerpunkt des Drahtquerschnittes seien 
mit I,, baw. Iz, das entsprechende polare Tragheitsmoment mit I) bezeichnet. 

Nach Biezeno-Grammel! laBt sich eine solche Feder unter den iiblichen vereinfachenden Annah- 
men (geringe Steigung, Vernachlassigung der Kontakteffekte) als weicher Stab behandeln. Aus der 
Verformung einer Federwindung unter verschiedenen Belastungen lassen sich die Biegesteifigkeit a, 
die Schubsteifigkeit B, die Drucksteifigkeit y und die Torsionssteifigkeit 0 berechnen. Diese vier 
GréBen, die das elastische Verhalten des Ersatzstabes beschreiben, hangen im Gegensatz zu den 
Steifigkeiten eines realen Stabes wesentlich von der relativen Lange I/I, ab. Es gilt namlich 


& =X I/Iy, B= Bollh, ¥=YVd/b, 6 = 69 1/ly, (2.1) 


wobei X», fy, Yo, 09 Konstanten sind und die Steifigkeiten der ungespannten Feder darstellen. Man 
findet fiir ihre Kehrwerte 


1 n m+1 J,\2R 1 n a R8 

ee! (1 2k Ee aes os 

a Ik die TE Po lL £1, (2.2) 
1 n,m+147R 1 _on 2aR | j 
le ie ae (bal eed Sn 


Dabei bezeichnet E den Elastizitatsmodul, m die Querzahl des Materials, n die Windungszahl der 
Feder und k einen Profilbeiwert, der fiir kreisformige Drahtquerschnitte gleich und fiir andere Quer- 
schnittsformen gréBer als eins ist. 
Bei der Diskussion der Knickgleichungen in Abschnitt 4 treten noch die GréBen 
1 Wo ie the 


—, PO 
Gow He re) 8 ae 


(2.3) 


~ Yo 1 m 5 
breed row 4Akm-+1 I, 


auf. Im folgenden beschranken wir uns immer auf Federdraht mit der Querzahl m = 10/3 und haupt- 
sichlich auf kreisférmige Drahtquerschnitte. Dann ergeben sich fir die Konstanten a und 6 die 


Werte 


— 23/26 = 0,8846 , 
ieee (2.4) 


b = 8/13 =0,6154. 


Fiir rechteckige Drahtprofile mit den Rechteckseiten c bzw. yc parallel zur ersten bzw. zweiten 
Hauptrichtung (Abb. 2) erhalt man fiir k, a und b die Werte der Tabelle 1. 


Tabelle 1 
lu k a b 
2 ———. 
PLA) 0 fore) 1,2692 1 
1 


1,820 1,0283 0,8679 


/ 
é j i 1,182 0,8254 0,6746 
Abb, 2. Eine halbe rwindung mit r ige 2 ? ? 
vie eakeaacsechiitael © desea 2 | 1,820 |) 0,6321 | 0,4717 


foe) lore) 0,5 0,2308 
Im folgenden werden nun die Knickformeln fiir einen durch Druck und konservative Torsion 
belasteten und durch die vier Steifigkeiten «, f, y, 0 elastisch gekennzeichneten Stab unter Beriick- 
sichtigung der Querkraft und der Schiebung hergeleitet, und erst bei der Diskussion dieser Glei- 
chungen in Abschnitt 4 wird mittels der Formeln (2.1) bis (2.4) wieder auf die Feder zuriickgegriffen. 
b) Die Differentialgleichungen. Nach dem statischen Stabilitatskriterium wird die kri- 


tische Belastung fiir einen geraden Stab, der durch eine axiale Druckkraft P und durch ein vor dem 


Knicken axiales, konservatives Torsionsmoment W beansprucht ist, dann erreicht, wenn zum ersten- 
mal nicht nur fiir die urspriinglich gerade, sondern auch fiir eine ausgebogene elastische Linie Gleich- 
gewicht méglich ist. 

Zur Aufstellung der Differentialgleichungen der elastischen Linie fithren wir zuerst das raumfeste 
Koordinatensystem (x, y, z) (Abb. 3) ein, dessen z-Achse mit der urspriinglich geraden Stabachse zu- 
sammenfallt. Da fiir den Ersatzstab in jedem Querschnitt alle Achsen durch den Schwerpunkt 
Mauptachsen sind, brauchen die x- und die y-Achse nicht speziell festgelegt zu werden, wahrend der 


wo 


1 C, B. Biezeno u. R. Grammel, Technische Dynamik, 2. Aufl. Bd. I, S. 616, Berlin, 1953. 
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Ursprung fiir jeden Knickfall geeignet gewahlt werden soll. Denkt man sich dieses System im 
undeformierten Stab in einen beliebigen Querschnittsschwerpunkt parallel verschoben und dann 
die Deformation des Stabes vollzogen, wobei das System mit dem Stabquerschnitt kérperfest ver- 
bunden bleibt, so gelangt man zu einem zweiten Koordinatensystem (€, y,¢). Seine €-Achse steht 
zum Querschnitt normal, fallt aber wegen der Schiebung im allgemeinen nicht mit der Tangente 
der elastischen Linie in diesem Schnitt zusammen. 

Um die Knickgleichungen zu finden, kann man sich auf kleine Ausbiegungen und Schiebungen 
beschranken und somit annehmen, daf die ¢-Achse nur wenig gegentiber der z-Achse geneigt sei. Im 
Rahmen einer linearisierten Theorie laBt sich dann die gegenseitige Stellung der beiden K oordinaten- 
systeme, abgesehen von einer Translation, durch den Verdrehungswinkel x des Querschnittes und 
durch die x- und y-Komponenten eines Einheitsvektors 7 (n,,n,, 1) in der Richtung der Quer- 
schnittsnormalen festlegen (Abb. 4). Der Querschnitt selbst wird durch die von 0 aus gemessene 
Bogenlange t bestimmt. Unsere Aufgabe besteht in diesem Abschnitt darin, fiir n,(t), ny(t) und 
x(t) Differentialgleichungen und deren Lésungen zu finden, denn es wird sich zeigen, daB durch diese 
Funktionen im wesentlichen auch die elastische Linie, d. h. x(t), y(t), bestimmt ist. 


Abb. 3. Die beiden Koordinatensysteme Abb. 4. Die gegenseitige Stellung der Abb. 5. Der unter einer Druckkraft und 
(x, y, z) und (&, 7, ¢). beiden Koordinatensysteme. einem Torsionsmoment ausknickende 
36 cy Ersatzstab. 


Wie man nach H. Ziegler! aus den Formeln der Kreiseltheorie herleiten kann, gelten fiir einen 
beliebigen Vektor a, dessen Komponenten a,, a,, az, a, von erster Ordnung klein sind gegen a, und 
az, die Transformationen 


a,, = a, cos ¥ — a, sin ¥ + a,n,, 
a, = agsin ¥ + a, cos x + any, (2.5) 
a, — ar. 


Bezeichnet man die Ableitungen nach t mit ', so folgen aus (2.5) und der Voraussetzung kleiner 
Schiebungen speziell fiir den Tangenteneinheitsvektor v die Gleichungen 


. : 
Vv, =x = vz,C0OS¥Y — Uv, sINY + N,, 


vy =y' =v;esin ¥ + v, Cos Y + ny, (2.6) 
v, = il = Ue . 

Zur Beschreibung der elastischen Deformation des Stabes wird noch der Vektor % gebraucht, 
der als Ableitung nach t des vektoriellen Drehwinkels des (&, 7, ¢)-Systems gegeniiber dem (x, ¥, 2)- 
System definiert ist. Nach dieser Definition ist seine Komponente x; gleich dem spezifischen Ver- 
drehungswinkel des Stabes und nach (2.5) darf dafiir in erster Naherung auch x, gesetzt werden, so 
daB gilt 

: py (2.7) 
© Durchlauft das kérperfeste System die elastische Linie gleichférmig, so stellt % seine Winkelgeschwin- 
i digkeit dar; somit gilt fiir die wirkliche Ableitung des in diesem System konstanten Vektors n die 
- Gleichung 
| ON Ts (2.8) 


1 H. Ziegler, Z. angew. Math. Phys. 2 (1951) S. 265. 
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Zerlegt man (2.8) in Komponenten, so erhalt man 
Ny = Xy — 4, Ny, (2.9) 


wiahrend die dritte Gleichung in erster Naherung identisch erfillt ist. 
Der Stab sei am Stabende B (t =6) durch ein Moment W (W,, W,, W) und durch eine Hinzelkraft 


K (K,, K,, —P) belastet (Abb.5). Dabei sind K, und K, Reaktionen, bei Einspannungen sind es 
auch W, and W,,. Diese Dyname erzeugt im Querechnitt F eine Beanspruchung, die durch die 


Dyname R Cie R, eke )s M (M,,M,,M,) dargestellt wird (Abb. 3). In bekannter Weise ergeben 
sich dann die Glecheevihehedincuiee! fiir ein Stabelement 


dR dM = 
a ee (2.10) 
Die erste dieser Gleichungen kann mit der Randbedingung r = K sofort integriert werden, in 
Komponentendarstellung erhalt man 
R,=K,, R,=K,, R,=—P. (2.11) 
Zerlegen wir auch die zweite Gleichung von (2.10) unter Verwendung von (2.11) in Komponenten, 
so finden wir : 
M,=K,+y’' P, 
M, = — K,—«’' P, (2.42) 
M, = K,y' — K,«' ~ 0. 
Die Terme auf der rechten Seite der dritten Gleichung sind von zweiter Ordnung klein und kénnen 
also durch Null ersetzt werden. Aus der dritten Gleichung folgt mit (2.5) und der Randbedingung 
M(b) = W 
M,=M,=W. (2.13) 
Zur Herleitung der gesuchten Differentialgleichungen benétigen wir noch den Zusammenhang 
zwischen der Beanspruchung und der Deformation. Im Rahmen der Festigkeitslehre gilt 


M;=ax:, M,=«x,, M;=0x,, 


" 
1—] 2.14 
R, =f vg, R, = 8 v4, Re =/( ak ( ) 
Setzt man in der dritten Gleichung die Werte aus (2.7) und (2.13) ein, so folgt 
WHO, 
also 
yomt mit mx, =%, (2.15) 


da die Integrationskonstante durch geeignete Wahl der x-Richtung zu Null gemacht werden kann. 


Mit Hilfe der Formeln (2.5), angewendet auf M und x, lassen sich die ersten zwei Gleichungen 
von (2.14) ins (x, y, z)-System iibersetzen. Man findet leicht 


M, =« (x, — %-n,) + Men,, 
M, =o (%, — xz ny) + M;n, 


und, indem man x, und %, aus (2.9) entnimmt, schlieBlich mit (2.13) und (2.7) 


M,=—an,+ Wn,, an 
M = on, + Wn, 4 Cas 


In der gleichen Weise transformiert man auch die vierte und fiinfte Gleichung von (2. 14), Durch 
Anwendung von (2.5) auf Rund Verwendung von (2.6) und (2.11) ergibt sich 


=F ta(1 +5), 
pe eee ae (2.17) 
ay B T ny ( | 3) . 
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Aus diesen Beziehungen ersieht man, daf es berechtigt ist, n, und n, als die eigentlichen Unbe- 
kannten zu wahlen, denn sind diese und die Konstanten K,, und K, bekannt, so wiirde sich die 
elastische Linie selbst durch Integration von (2.17) bestimmen lassen. 

Setzt man jetzt (2.17) und die Ableitungen von (2.16) in (2.12) ein und ordnet, dann findet man 
die gesuchten Differentialgleichungen fiir n, und ny 


aon, + Wn, + P(l +3) =—K,(1 +3): 
2.18) 
any — Wr, + P(l+5)n,=—K,(t +5). ( 


Diese Gleichungen, sowie auch die Beziehungen (2.16), die noch fiir die Aufstellung der Randbedin- 
gungen gebraucht werden, kénnen mit den komplexen GréBen 


N=n,+in,, K=K,+ik,, M=M,+iM, (2.19) 
komplex zusammengefaBt werden und lauten dann 
; B : 
EN IW NE PLN K (147, (2.18a) 
M=aiN'+Wund«. (2.16a) 


Zur Vereinfachung der Schreibweise und um die Resultate besser diskutieren zu kénnen, fiihren 
wir dimensionslose GréBen ein. Wir ersetzen die unabhangige Variable t durch 


f= 7 (2.20) 


und versehen jede von 7 abhangige Funktion auch mit einem Querstrich. Fiir eine beliebige reelle 
oder komplexe Funktion f gelten dann die Beziehungen 


fO =f, 


- eum ldh ae le oe 
u d°f 1 af 1 Fis 
f(t) Gita Age h at (t). 
Dividieren wir die Gleichung (2.18 a) durch «//? und beniitzen noch die dimensionslosen GréBen 
PP Wi _ KP _ MI Wel V1 uke 
= 7a a? L= ore Lk aaa SR) ee a Ta 'y ae E BE (2.22) 


so lautet die gesuchte Differentialgleichung schlieBlich 


N"—2iwN +p(ltep)N=—k(l+ep), (2.18b) 
wahrend die Beziehung (2.16 a) in 


= i 


m=;N'+wN (2.16b) 
tibergeht. 
Mit den komplexen Integrationskonstanten 
A=4,+i4; B=B,+1B; C=C,+41¢, (2.23) 
findet man als allgemeine Lésung der Differentialgleichung (2.18 b) 
N=Acit' + Best +C. (2.24) 
Dabei ist 
Cae (2.25) 
Pp 
und 4,, A, sind die Wurzeln der charakteristischen Gleichung 
—A = 2ws | pUlzeep)— 0. (2.26) 
Mit der Abkiirzung 


q =)u? + p (1 + ep) (2.27) 
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ergeben sich dafiir die Werte 
Ap Weg - (2.28) 


Die partikulare Lésung (2.25) wird ungiiltig fiir den Grenzfall p = 0 (reine Torsionsbelastung)- 
Dann erhalt man als Lésung von (2.18 b) 


N = | e2iwt 1 B+Ct (2.24 a) 
mit 
es (2.25) 
2 w 


Damit wir die den verschiedenen Knickfillen zugehérigen Randbedingungen aufstellen konnen, 
miissen wir zuerst die verschiedenen konservativen Torsionsmomente einfiihren. 


c) Konservative Torsionsmomente. Um die konservativen Torsionsmomente zu finden, 
berechnen wir die virtuelle Arbeit 6 4,, der an den Enden B und A des Stabes angreifenden Torsions- 
momente W, (Wie Voy Wi. = +W) und W,(W,,, Way, Wi,=—W). Man findet dafir nach 
H. Ziegler}, indem man dort itberall wegen der Schiebung den Tangentenvektor v durch den Normalen- 


vektor n ersetzt, 
b 
6A, = W f (n, 6 n,— n, 0 n, +6 xz) dt + (W,, 0 n,, — W,, fs) fe) str ay. 0 gz — Wz 0 We,)- 


Formt man diesen Ausdruck noch etwas um, so kommt man auf 


b 
, re 1 1 
os W | (1, ny — ny Mz + 2 x) a se (May 3 Me Muy) Mp x ee mua "us O My 


a 


ae (ay ae = We "as ) a es a Way "es 6 Nay + (2.29) 


Die Torsionsmomente Wz und W,, sind aber dann und nur dann konservativ, wenn sich auch die 
Randterme als Variation eines Potentials darstellen lassen. Hierfiir ergibt sich die Bedingung 


0 1 Om/al 
any (Vs — a Me s| = Goa Wn, W.) 


die an beiden Stabenden gelten mu. (Die Zeiger a und b sind hier und im folgenden der Einfach- 
heit halber weggelassen.) Wertet man diese Bedingung gemafB dem Vorgehen von M. Beck? mit emem 
linearen Ansatz 


WW AA* n, + BU n,)5 2) We (Cre cd aorta) (2.30) 
aus, so erhalt man fiir konservative Momente die Beziehung 
AS DY =a, (2.31) 


Von diesen allgemeinen konservativen Momenten (2.30) und (2.31) greifen wir drei Typen heraus, 
die sich verhaltnismafig einfach realisieren lassen, namlich 


1. das semitangentiale Moment, mit 


x z x 


1 
A*=D* ==, B*=C*=0; alo Wi25W,n,, Wy Wine 


2. das quasitangentiale Moment, mit 
A* = 1, BY— C= D*=0: aloo = Wine.) Wo 0r (2.32 q) 
3. das pseudotangentiale Moment, mit 
A* =|, B*=tgz, C* = D* =) 0:2 ealso W, = W,(n, + tgtn,), W,=0. (2.32p) 
Nach H. Ziegler’ (wobei man auch hier in unserem Falle wegen der Schiebung den Tangenten- 
vektor v durch den Normalenvektor n zu ersetzen hat) lassen sich diese Momente so realisieren, dab 


man in den Fallen 1) und 2) den Endquerschnitt des Stabes starr mit einer Kreisscheibe verbindet 
und an dieser gleichmafig verteilte, tangentiale Krafte von konstantem und gleichem Betrag K 


* Vel. die in FuBnote 7 von Seite 113 zitierte Arbeit S. 104—105. 
® M. Beck, Ing.-Arch. 23 (1955) S. 231. 
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und konstanter Richtung angreifen laBt, und zwar fiir das quasitangentiale Moment zwei und fiir das 
semitangentiale Moment drei oder mehr (Abb. 6). 


Um das pseudotangentiale Moment zu erzeugen, ist der Stabquerschnitt starr mit einem Quer- 
balken zu verbinden, an dessen Enden ein Kraftepaar mit Kraften von konstantem Betrag und kon- 
stanter Richtung angreift. t bedeutet dabei den Drehwinkel des Querbalkens, die vertikale Kom- 
ponente W des Momentes hat dann 
den Betrag Zz Zz 

W—=2aK cost. — (2.33) 


Je nachdem, ob der Stab am 


zweiten Ende a) st i P 
) starr eingespann nat PRK 


ist, oder b) sich axial drehen kann, 


K 


ist t gleich dem ganzen oder halben K 
totalen Verdrehungswinkel 

oO= x, 1 (2.34) 
des Stabes. Mit (2.15) erhalt man 

te y z Y 

also 

ee ear et semitangentiales guasitangentiales pseudotangentiales 

i ee: (2.35 a) Moment Moment Moment 

ae fi o]2 x = (2.35 b) Abb. 6. Die drei konservativen Momente. 


Weil aber t nicht gréBer als 7/2 werden kann, folgt daraus, daB auch W beschrankt ist, namlich 


) 
Wi 1° (2.36 a) 
bzw. 
) 
SchlieBlich ist noch zu bemerken, daB wegen (2.30) im allgemeinen weder das axiale (A* = B* = 
Cc* = D* = 0), noch das tangentiale (A* = D* = 1, B* = C* = 0) Moment konservativ sind. 
Sind freilich beide Stabenden eingespannt (Lagerungen I und II), so sind dn,, = 6 ny =O m, = 
6 n,, = 0 und die Randterme in (2.29) fallen weg, das belastende Moment ist dann immer konser- 
vativ, und es 1a8t sich nicht mehr zwischen axialen, semitangentialen und quasitangentialen Mo- 
menten unterscheiden; nur das pseudotangentiale Moment wahrt wegen (2.36a,b) noch eine 
gewisse Sonderstellung. 


3. Die Knickgleichungen. Wie jedes Knickproblem ist auch dieses mathematisch ein Eigenwert- 
problem. Es handelt sich darum, die Kigenwerte p und w so zu wahlen, daB die Lésung (2.24) unter 
Beriicksichtigung der Randbedingungen nicht identisch verschwindet. Zu jedem Knickfall gehéren 
drei komplexe oder sechs reelle Randbedingungen, die fiir die drei komplexen Unbekannten A, B 
und (€ oder ihre sechs Real- und Imaginarteile ein lineares, homogenes Gleichungssystem liefern. 
Ein solches Gleichungssystem mit gleich viel Unbekannten wie Gleichungen besitzt dann und nur 
dann eine von Null verschiedene Lésung, wenn seine Determinante verschwindet. Durch diese 
Forderung erhalt man die gesuchte Knickgleichung, die eine Beziehung zwischen den kritischen 
Werten p und w darstellt. 


a) Knickfall I. Fir den beidseitig eingespannten Stab wahlen wir den Koordinantenursprung 
des (x, y, z)-Systems in der Stabmitte, so.daf den Stabenden A und B die dimensionslosen Bogen- 
langen t = —1/, und t= + 1/, zugeordnet sind. Die Einspannungen an den Enden des Stabes er- 
geben auf Grund der Definition des Vektors n und mit den Bezichungen (2.19) und (2.21) 


a(t)-a(-4)0 an 


PA; Ag AY . As 


Ae®tBe*tO=0, Ae ?7iBe 7? +0=0. (3.2) 
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Ferner ist zu fordern, daB die Verschiebung zwischen den beiden Stabenden in der x- und der 
y-Richtung Null ist, d. h. daB 
1/2 1/2 
fxidt= fy dt=0 (335) 
—1/2 1/9 
gilt. Indem man wieder auf komplexe und dimensionslose GréBen iibergeht, kann man mit (2.17), 


(2.19), (2.21) und (2.22) fiir (3.3) 


1/2 1/2 7 
f(@e +iy)dt = f [N(Q_+ep)+ek]dt =0 (3.4) 
—1/2 —1/2 
schreiben, woraus sich schlieBlich durch Einsetzen von (2.24), Ausfiihren des Integrals und Multi- 
plizieren der Gleichung mit 2, A, = — p (1 + € p) (wegen (2.28) und (2.27)) die Beziehung 


if Seeds Par apes 
glee 2) 4 pale? —e Capen (3:5) 
ergibt. Die gleich Null gesetzte Determinante des Gleichungssystems (3.2) (3.5) liefert jetzt unter 
Verwendung von (2.27) und (2.28) die Knickgleichung 

2 q (cos w — cos q) — psing = 0. (1) 
b) Knickfall II. Sind die beiden Stabenden gegeneinander horizontal verschiebbar, so treten 

keine Querkrafte auf. Wegen (2.19), (2.22) und (2.25) folgt sofort 
C=0. (3.6) 
Wahlt man hier den Koordinantenursprung des (x, y, z)-Systems im unteren Stabende, so ergeben 
die wie im Fall I hergeleiteten Einspannbedingungen N(0) = N(1) = 0 mit (2.24) und (3.6) die 


Gleichungen 


As B=0, Ae*-— Bes = 05 (3.7) 
Setzt man die Determinante dieses Systems gleich Null, so folgt mit (2.27) und (2.28) 
w?>tp(l+tep)=h2, (al See) 
und wenn wir uns hier, wie auch im folgenden, auf die kleinsten Knickwerte beschranken 
w+ p(l+ep)=277. (IT) 
c) Knickfall III. Wie im Fall II gilt auch hier wegen des Wegfallens der Querkrafte die Glei- 
chung (3.6) und damit wegen der Einspannung am unteren Stabende — das auch hier Ursprung 
des (x, y, 2)-Systems ist — die Gleichung 
A+B=0. (3.8) 


Das belastende Torsionsmoment liefert die dritte Gleichung, doch haben wir zwischen den 
verschiedenen Belastungsarten zu unterscheiden. 

x) Semitangentiales Moment. Setzt man in der ersten und zweiten Komponentengleichung 
der Randbedingung 


M(l) = W (3.9) 
die Werte fiir das semitangentiale Moment (2.32 s) ein und geht mit (2.19), (2.20), (2.21) und (2.22) 


wieder auf komplexe und dimensionslose GréBen iiber, so entsteht die Beziehung 


wo — 


m(1) = 3 N(1). (3.10) 
Mit Hilfe von (2.16 b) entsteht aus (3.10) die Gleichung 
; N'(1) —iw M1) =0. (3.11) 


Setzt man jetzt die Lésung (2.24) unter Beriicksichtigung von (3.6) und (3.8) in (3.11) ein, 
so erhalt man mit (2.28) und (2.27) die Knickgleichung 


ead 


wis P (Ler ep) =a (III s) 


6) Quasi- und pseudotangentiales Moment. Diese beiden Falle kénnen gemeinsam be- 
handelt werden, denn rein formal kann das quasitangentiale Moment (2.32 q) aus dem pseudotangen- 
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tialen (2.32 p) erhalten werden, indem man tgt durch Null ersetzt. Fiihrt man diese Substitution 
erst in der pseudotangentialen Knickformel aus, so erhalt man die Knickgleichung fiir das quasi- 
tangentiale Moment und eine direkte Berechnung eriibrigt sich. 

Fiir das pseudotangentiale Moment (2.32 p) ergeben sich in dimensionsloser Schreibweise mit 


(2.20), (2.21) und (2.22) die Komponenten 
w, =w[n(1) + tgtn,(1)], w,=0. (3.12) 


Da sich diese Komplex zusammengefaBt nicht durch N ausdriicken lassen, ist es nétig, die Gleichung 
(2.26 b) gemaf (2.19) in Real- und Imaginarteil zu zerlegen. Es ergibt sich 


es 


m0) =— > m(F) + wii,(t), ——my(?) = 4 nil) + waa,(2), (3.13) 


und indem man die Randbedingung (3.9) mit (3.12) und (3.13) auswertet, entstehen die Gleichungen 
ny(1).+ 2 wtgtn(1) =0, nl) + 2wn(1) = 0. (3.14) 


Zerlegt man die Lésung (2.24) mit (2.23) in Real- und Imaginarteil und setzt sie, sowie ihre Ableitung, 
unter Beriicksichtigung von (3.6) und (3.8) in (3.14) ein, so bekommt man mit (2.28) die beiden 
Gleichungen 


A, [A, cos A, — A, cos dA, + 4w tg t sing cos w] + A; [—A, sind, +A, sind, —4 w tgt sing sin w] = 0, 
A, [A, sind, —A, sin A,] + A, [A, cos A, —A, cos A,] = 0. 


Durch Nullsetzen der Determinante erhalt man mit (2.27) und (2.28) die Knickgleichung 


P(l+ep) + (w+ q)cos2qg+2wqtgtsn2g=0, 
die durch goniometrische Umformung und Abspaltung des Faktors sin? q auf die Form 
qg ctg?q +2wqtgtctgq—w? = 0 (IIT p) 
gebracht werden kann. | 


Indem man tgt durch Null ersetzt, erhalt man jetzt auch die Knickgleichung des quasitangen- 
tialen Momentes 


q? ctg?q@— uw? = 0. (IT q) 

d) Knickfall IV. Da hier infolge des Auftretens von Querkraften C nicht mehr Null ist, ergibt 
die Einspannbedingung am unteren Stabende (t = 0) statt (3.8) 

A+B+C=0. (3.15) 


Eine zweite Gleichung folgt aus der Bedingung, da die Stabenden in der x- und der y-Richtung 
nicht gegeneinander verschiebbar sind. Diese wird wie im Falle I hergeleitet, nur da® den Stab- 
enden jetzt die Bogenlangen 0 und 1 zugeordnet sind. Die Rechnung fiihrt auf 

AA, (e'4 — 1) + B(A, e'*= —1)—Cip=0. (3.16) 
Fir die dritte Randbedingung hat man wieder zwischen den belastenden Momenten zu unter- 
scheideh. 


x) Semitangentiales Moment. Genau gleich wie im Falle IIs fihrt die Belastung durch ein 
semitangentiales Moment auf die Randbedingung (3.11) oder ausgerechnet mit (2.24), (2.27) und 
(2.28) auf 

Aqes—Bqex—Cw=0. (3.17) 

Setzt man die Determinante der drei Gleichungen (3.15), (3.16) und (3.17) gleich Null, dann 

entsteht mit (2.27) und (2.28) die Knickgleichung 
pqcos q— (wv? + q’) sing +2wqsinw=0. (IVs) 


B) Quasi- und pseudotangentiales Moment. Von Fall IIIp kann die Randbedingung 
(3.14) tibernommen werden. Setzt man die zerlegte Lésung (2.24) und ihre Ableitung in (3.14) ein 
und driickt © gemaB (3.15) durch A und B aus, so ergeben sich mit (2.28) die beiden Gleichungen 


A_[A, cosa, + 2wtgr sind] + A; [—A, sind, + 2 w tgt (cos A, — 1)] 
+ B [A, cosd, + 2 w tgt sind] + B;[—A, sind, + 2 w tgt (cosA, — 1)] =0, (3.18) 
A, Ay sind, + A; [A, cos A, — 2 w] + B,A, sind, + B; [A, cos A, — 2 w] = 0. 
g* 
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Zwei weitere Gleichungen fir A,, A4;, B, und B; liefert die Beziehung (3.16), wenn man sie unter 
Beriicksichtigung von (3.15) in Real- und Imaginarteil zerlegt. Durch einige Umrechnungen dieser 
vier Gleichungen gelangt man schlieBlich zur Determinantenbedingung 


WY Ay : 

A, (1 — cos A,) p tA, sind, 

A, sind, + p A, (cos A, — 1) 

A, cosd, +2 w tgt sind, A, sind, + 2 w tgt (cos A, — 1) 

i, i Bik | 

A, (1 — cos A) p +A, sin A, E 

rf sind, + p A, (cos A, — 1) | 0. (IVp) 
Az cos Ay + 2 w tgt sind, —A, sind, + 2 w tgt cos (A, — 1) 


Die Knickgleichung fiir das quasitangentiale Moment erhalt man wie im Fall III, indem man 
tg t durch Null ersetzt, 


p Ay Pp hs 

Ay (1 — cos A) P +A, sind, A,(l—cosA,) p+tA,sinA, | _ 0 (IVq) 
A, sind, + p A, (cosA,—1) A, snd, + p A, (cosAg—1) | q 
A, cos A, —A, sind, Ay C08 Ay —/, sin A, 


e) Knickfall V. Sind beide Enden des Stabes gelenkig gelagert, so kénnen an beiden Enden 
konservative Momente, die nicht axial bleiben, angebracht werden. Wir beschranken uns hier auf 
Falle, wo an beiden Enden gleichartige Momente angebracht werden. Ferner sollen beim quasi- 
und beim pseudotangentialen Moment die erzeugenden Krafte an beiden Stabenden parallel sein. 
Mit dem Koordinatenursprung des (x, y, z)-Systems in der Stabmitte findet man als erste Rand- 
bedingung wie im Falle I die Gleichung (3.5). Fiir die weiteren Randbedingungen treffen wir die 
bekannte Fallunterscheidung. 

ax) Zwei semitangentiale Momente. Die Bedingungen fiir die Momente lauten an den 
beiden Stabenden allgemein 


a 1 is sth l = 
M(t) = + Wp, it(—1) -_W,, (3.19) 
wobei zu beachten ist, daf fiir die vertikalen Komponenten gilt W,, = W und W,, = — W. Durch 


Einsetzen des semitangentialen Momentes (2.32s), Ubergang auf dimensionslose und komplexe 
GréBen und Verwenden der Lésung (2.24) erhalt man vdllig analog zum Fall [Vs die beiden (3.17) 
entsprechenden Gleichungen 


Ay As Ay Ae 


INGO I HO. FOr al Age =— Bae 7>_Cw=0. 
Zusammen mit (3.5) fithren sie auf die Knickgleichung 
2 w® cos w — 2 w* cos gq — pqsing = 0. (Vs-s) 


B) Zwei quasi- oder zwei pseudotangentiale Momente. Wertet man die Randbe- 
dingungen (3.19) fiir das pseudotangentiale Moment (2.32 p) aus, so ergeben sich mit dem Vorgehen 
des Falles II p die zu (3.14) analogen vier Gleichungen 


na (tg) +2wi, (+ 3) =0, my (+g): 2w teri, (+ 5) =0. 


Durch Einsetzen der Lésung (2.24) und ihrer Ableitung, sowie durch geeignete Addition und Sub- 
traktion von Gleichungspaaren, erhalt man daraus mit (2.28) die vier Beziehungen 


A 
A, A, cos 2 + Bh cos 2 +-C,2w =0, 
aa 8 abe 
A, i, sin rate B.A, sin > 08 


a 


A Aes 
A, cos > + 2 w tgt sin > +B, |izcos F +2 utersin 3] = 0, 


ag. Uh A eed 
A, [2s sin — 2 tg cos 3 + B; sin —2w tgr cos 3]—C,2wtgr =O. 
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Dazu kommen noch die beiden Gleichungen, die durch die Zerlegung der Beziehung (3.5) entstehen. 
Die Determinante dieses Systems von sechs Gleichungen lat sich nach einigen elementaren Um- 
formungen in ein Produkt zerlegen. Man findet so als Knickgleichung fiir den pseudotangentialen 


Fall 

{(w> +- q?) sin g — 2 wq sin w + 2 w tgt q (cos w — cos q)} 

{2 w tgt [pq (cos q + cos w) + 4 wq sin w — (2 w? + 2 q?) sin q] (Vp-p) 

+ 8 w? q (cos w — cos q) — p(w? + q?) sing —2 pwqsinw} =0, 

und daraus fiir den quasitangentialen Fall 
{(w? + q?) sing —2 wqsin w} {8 w? q (cos w — cos q) — p (w? + q?) sing —2 pwqsinw} =0. (Vq-q) 
: Bevor wir die elf Knickgleichungen genauer diskutieren, bemerken wir noch, daB sie fiir unend- 
liche Schubsteifigkeit 6, also nach (2.22) und (2.27) fiir e = 0 und q = (w? + p)!/2 in die entsprechen- 
den Formeln fiir den schlanken Stab iibergehen, die von M. Beck! berechnet wurden. Ferner stim- 


men die Knickgleichungen in den Fallen I und II iiberein mit den Gleichungen, die Haringx® ge- 
funden hat. 


4, Auswertung und Diskussion. Um die Diskussion bequem durchfiihren zu kénnen, miissen 
in den Knickgleichungen Variable auftreten, die zu den eigentlichen kritischen Lasten P und W 
proportional sind. Die GréRe w erfiillt diese Bedingung und wird deshalb als die das Torsionsmoment 
W charakterisierende Variable verwendet, denn nach (2.22) und (2.1) gilt 
Wl, 


wWw=—-aza— - 
2 & 


(4.1) 


In gewissen Fallen wird W auch durch den totalen Verdrehungswinkel w ausgedriickt werden, denn 


mit (2.34), (2.15), (2.1) und (4.1) erhalt man 
? 2 W il, _ 2a 
i «ae, aah (4.2) 


Dagegen ist die Gré8e p nicht zur Druckkraft P proportional, weil in (2.22) der Beiwert [?/x 
selbst mit P variiert. ZweckmaBigerweise verwendet man zur Charakterisierung der Druckkraft P 
in den Knickgleichungen die negative Federdehnung 


g=1—{. (4.3) 


Die letzte Gleichung von (2.14) lautet namlich dann mit (2.11), (4.3) und (2.1) 


1 Sa a (4.4) 


also ist £ zu P proportional. Die in den Knickgleichungen auftretenden GréBen p und p (1 + ¢ p) 
lassen sich mit (4.3), (2.22) und (2.1) durch € und Konstanten ausdriicken. Man findet 


p=B4a—H8, p(ltep) =" o4 at 4 , 


Diese Formeln lassen sich unter Verwendung der nur von der Drahtquerschnittsform und dem 
Material abhangigen Parameter a und b (2.3) sowie durch Einfiihrung des Schlankheitsgrades 


(4.5) 


l 
S= oan (4.6) 
als 
p=H4as(1—éé, pit+ep)=4as(1—béjé (4.7) 
schreiben. 


Bevor wir uns der Diskussion der allgemeinen Knickformeln zuwenden, wollen wir noch die 
Spezialfalle untersuchen, die sich durch Weglassen des Torsionsmomentes oder der Druckkraft er- 
geben. 

a) Reine Druckbeanspruchung. Setzt man in den verschiedenen Knickformela iiberall w 
und gleich Null, so ergeben sich unter Verwendung der Substitution (4.7) die folgenden Gleichungen 
zwischen dem Schlankheitsgrad s und der negativen Federdehnung &. 


1 Siehe FuBnote 2 von Seite 118. 
2 Siehe FuBnoten 4 u. 5 von Seite 113. 
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Knickfall I: va(l—bépé =z’. 

Knickfall II: Asta(l—béEé =z". 

Knickfall II: 16s a(l1—bé)E =a’. (4.8) 
ae ari ites ays 

Knickfa g s Va ( E)E jada 


Knickfall V: 4Astta(l—bé)E=z7. 


Natiirlich kann man diese Gleichungen, die mit den von Haringx gefundenen Formeln identisch 
sind, sowohl aus der semi-, quasi- oder pseudotangentialen Knickformel des entsprechenden Falles 
herleiten, da mit verschwindendem Torsionsmoment diese Unterschiede hinfallig werden miissen; 
jaim Knickfall V, wo sich die quasi- und die pseudotangentiale Knickformel in zwei Faktoren zer- 
legen lassen, zeigt es sich, daB beide Faktoren auf die Gleichung (4.8V) fiihren. Der durch die 
Gleichungen (4.8) gegebene Zusammenhang zwischen é und s ist in Abb.9, 10, 11 und 12 durch die 
Kurven w = 0 graphisch dargestellt; die Diskussion wird im Rahmen des allgemeinen Knickpro- 
blems spater durchgefiihrt. 


b) Reine Torsionsbeanspruchung. Setzt man in den Knickgleichungen p gleich Null, 
dann erhalt man die kritischen Momente fiir reine Torsionsbeanspruchung, die zur Unterscheidung 
vom allgemeinen Fall mit dem Zeiger , versehen werden. Die Tatsache, daf die allgemeine Lésung 
(2.24) fiir diesen Spezialfall nicht richtig ist, fiihrt dazu, daB die meisten Knickgleichungen fir 
p = 0 zu Identitaten werden. Man kann die Knickwerte trotzdem daraus erhalten, indem man die 
Knickgleichungen fiir kleine p entwickelt oder aber, indem man die Rechnungen von Abschnitt 3 
nochmals durchfiihrt, jedoch die Lésung (2.24a) verwendet. Setzt man nach (2.22) 


e(p=0) =m = 5%: (4.9) 


so erhalt mit beiden Methoden fir das kritische Torsionsmoment wy, bei reiner Torsionsbeanspru- 
chung in den verschiedenen Knickfallen die folgenden Gleichungen 


Knickfall I: 


W, 
18 = Ta a? (4.101) 
Knickfall I: 
Wy = 2; (4.10 IT) 
Knickfall III semitangential: 
mats, (4.10 IIs) 
Knickfall II quasitangential: 
wm =+ 5, (4.10 III q) 
Knickfall III pseudotangential: 
—>=+4, (4.10 III p) 


Knickfall IV semitangential: 
1+ w? 


Cte, == esa ree (4.10 TVs) 
Knickfall IV quasitangential: 
sin? Wy — Wy sin 2 wy + we — wé cos 2 wy 
“+ & (4 w§ sin 2 wy — 8 w) + 8 w3 sin 2 wy) + 16 &2 w8 cos 2 wy = 0, (4.10IV q) 


Knickfall IV pseudotangential: 
ae ‘ : : 
SIN" Wy — Wo Sin 2 wy + wh — ué cos 2 wy + & (4 w3 sin 2 wy — 8 wo + 8 we sin 2 wy) 
+ 16 ¢5 wp cos 2 wy + tg t (2 wy sin® wy — 2 w2 sin 2 wy + 2 we — wi sin 2 wo) 
+b & tg t (8 wo sin? wy — 8 uf cos 2 wo) + 16 ek tg 7 wS sin2 wy = 0, (4.10IV p) 
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Knickfall V semitangential: 


a Wo 
tg Wo = 5 reas (4.10 Vs) 
Knickfall V quasitangential: 
erster Faktor Wy = + = : 
zweiter Faktor Wy = + ee ‘ (4.10 V q) 
= &y 
Knickfall V pseudotangential: 
erster Faktor Wy — T= S > 


Se : 
zweiter Faktor tgw, = —* Abas 6 80) 08 * (4.10 V p) 


Wy (1 — 4 w2 €9) + ter’ 

Um diese Gleichungen bequem diskutieren zu kénnen, fiihren wir noch die GréBe ¢, nach (4.9) 

mittels der Beziehungen (2.3) und (4.6) auf den Schlankheitsgrad s zuriick. Man findet 
1—b 


éy = >: 
0 4as* 


(4.11) 


Fiir kreisférmige Drahtquerschnitte (2.4) sind die Funktionen w,(s) fiir semi- und quasi- 
tangentiale Momente in Abb. 7 graphisch dargestellt. In den Knickfallen IT und III, wo keine 
Querkrafte und somit keine Schiebungen auftreten, sind die Werte w, konstant und stimmen er- 
wartungsgema8 mit den Knickwerten fiir den schlanken Stab iiberein. In allen andern Fallen 
stimmen nur die Grenzwerte fiir s = oo mit den Knickmomenten fiir den schlanken Stab itiberein, 
wahrend fiir endliches s die KnickgréBen wy kleiner sind und mit zunehmendem s monoton zunehmen. 


70 
if 
Y 
IVs 
| 
w/G 
3 
| i—Vs+ 7 a =a | 
Wq 
aA A ____ 
} Is, Vg (1.faktor) 
Al 
‘ A [Ve (z. Fakior) 
LY = 
/ ; 
0 2 Y 6 8 10 


ie Kni ui ine Torsionsbeanspruchun i i- sitangentialem M : 
Abb. 7, Die Knickkurven I bis V fiir reine Torsionsl pruchung bei semi- und quasitangentialem Moment 


Dabei ist freilich zu beachten, daB der Proportionalitatsfaktor 1/2 Xo zwischen dem eigentlichen 
Torsionsmoment W, und der GréBe wy selbst mit dem Schlankheitsgrad variieren kann, a ie 
ist 1,/2 xp nach (2.2) bei festem Drahtquerschnitt zum. Wicklungsradius R Papen anes und una ss 
hangig von /,. Beriicksichtigt man diese Abhangigkeit, so ergibt sich, dab W bei peers naere 
in allen Fallen mit zunehmendem R abnimmt, wahrend bei konstantem Rin den Fallen II un 
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W,, fiir alle 1, fest bleibt, in den andern Fallen aber mit zunehmendem /, zunimmt. Andert man R 
und [, proportional, so bleiben s und wy fest, Wo aber nimmt mit zunehmendem R ab. Durch nicht- 
proportionales Verandern von J, und R kann schlieBlich auch erreicht werden, da W, mit zuneh- 


mendem s abnimmt. . 
Die Werte, zwischen denen wy variieren kann, sind in der folgenden Tabelle 2 zusammenge- 


stellt. (Dabei ist selbstverstandlich auch der Wert s = 0 nur ein mathematischer Grenzwert, der 
sich mechanisch nur niherungsweise realisieren liBt.) 


Tabelle 2. 


Fall Wy min Wo max Bestimmungsgleichung fiir Wo ee 
(s = 0) (s = 0) 
u 7 4,494 tg wy = wy 
Il 1 PS 
IIIs /2 a/2 
Ill q m/4 m/4 
1 + wé 

IVs m/2 3,406 clg wy = ara hal 

0 
IVq | 2/4 2,476 sin? wy — wy sin 2 wy + wi — w§ cos 2 wy = 0 
Vs | 0 2,456 tz w) = — * 
Vg i 0 n/2 


Besondere Beachtung verdient noch der Fall Vq. Die Knickgleichung zerfallt in zwei Faktoren, 
die beide beriicksichtigt werden miissen. Im Bereich s > 1,036 wird die Knickkurve durch die 
Gerade wy = 2/2 dargestellt, die schon im Fall IIIs auftrat; fiir s < 1,036 fallen die kritischen 
Werte dann aber mit s linear gegen Null ab. Daraus folgt aus Stetigkeitsgriinden, daB auch bei der 
allgemeinen Druck- und Torsionsbeanspruchung mindestens fiir kleine Druckkrafte keiner der 
beiden Faktoren zum vorneherein weggelassen werden darf. 

Fiir Federn mit andern als kreisférmigen Drahtquerschnitten, also mit andern Werten a und b 
kénnen die Kurven der Abb. 7 auch verwendet werden, wenn man nur die Abszisse s durch den 
Wert s ee *” ersetzt. Man findet damit z. B., daB fiir in Richtung der Schraubenachse lang- 
gestreckte Rechteckprofile (a und b kleiner als beim Kreisquerschnitt, Tabelle 1) die kritische Tor- 
sionsgréBe wy fiir einen festen Schlankheitsgrad im allgemeinen kleiner ist als beim Kreisquerschnitt. 

Komplizierter liegen die Verhiltnisse bei den pseudotangentialen Momenten, da zu den 
durch die Feder bestimmten Parametern a, b und s auch noch das Verhaltnis 


ay 21, 


ee zat panes 4.12) 
Oo ( - m-+1 I, ( 
I (1+ 2h" 2 i.) 
hinzutritt, denn gemaé8 (4.2) und (2.35a und b) gelten die Beziehungen 
2a 
Ve ee an Wo > (4.13a) 
— O = he 
tae oe (4.13b) 


je nachdem, ob die Feder am unteren Ende a) starr eingespannt ist oder b) sich um eine vertikale 
Achse drehen kann. Im Knickfall V kommt natiirlich nur b) in Frage, da dort nach Voraussetzung 
die Momente an beiden Federenden pseudotangential sind. Ferner miissen immer die Ungleichungen 
(2.36a und b) beachtet werden, die fiir wy und 7 obere Schranken abgeben. 

Fiir die Knickfalle I und IT andern sich die Knickgleichungen nicht, die Ungleichungen 
(2.36 a und b) ergeben aber mit (2.22) und (2.1), daB fiir o/d) = 1/4 (Fall a) bzw. a/d) => 1/2 
(Fall b) tiberhaupt kein Knicken eintritt. Speziell fiir den Kreisquerschnitt (9/69 = 0,8696) ist 
Knicken in beiden Fallen unméglich. 

Abb. 8 zeigt die Abhangigkeit wy (x9/d) fiir den Knickfall IIIp, die KnickgréBe wy ist hier 
unabhangig vom Schlankheitsgrad. Man findet durch Rechnung oder aus der Figur, daB fiir den 
Kreisquerschnitt nur im Falle b Knicken méglich ist und zwar bei einem Moment Wo —= 1,390; 

; Die Untersuchung der Gleichung (4.101 Vp) zeigt, daB im Knickfall IV p gleich wie im Fall IIIp 
Knicken nur méglich ist fiir x9/d) < 1/2 (Fall a), baw. 9/59 < 1 (Fall b), aber bei einem bestimmten 
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Querschnitt auch nur fiir Schlankheitsgrade s < s*(x9/d9, a, b). Dieser gréBte Schlankheitsgrad s* 
wird fiir o%/g S 0,231 (Fall a) bzw. o/5) < 0,461 (Fall b) unabhangig von a und b gleich unend- 
lich, so daB® unterhalb dieser Werte Knicken fiir jeden Schlankheitsgrad méglich ist. Fiir den Kreis- 
querschnitt ergibt sich daraus, da8 Knicken nur im Falle b méglich ist und zwar nur fiir s < s* = 
0,3928 und mit einem Moment wy < 1,806. 

Beim Knickfall Vp liegen die Verhaltnisse so, daB auch hier Knicken nur fiir »/d) < 1 még- 
lich ist, und zwar liefert fiir 1/2 < «,/d) < 1 der zweite Faktor allein das Knickmoment. Dabei ist 
_ Knicken wieder nur méglich fiir s < s**(ao/5p, a, 6). (Kreisquerschnitt s** = 0,655 und wy < 1,806). 
Wenn a/d9 < 1/2 ist, so ist der erste Faktor fiir groBe und der zweite fiir kleine Schlankheitsgrade 
maBgebend, Knicken ist dann immer méglich. 

c) Kombinierte Torsions- und Druckbean- 7% 
spruchung. Um die allgemeinen Knickgleichungen 
aus Abschnitt 3 fiir eine bestimmte Feder auszuwerten, 
hat man die Federgré8en a, b, s (bei pseudotangentialen 
Momenten auch ,/d9) in die Knickgleichungen einzu- 
fiihren und fiir variables € die entsprechenden w- Werte Falla: Fall b 
zu berechnen. Da die Knickgleichungen auBer in den 
Fallen II und IIIs transzendent sind, ist diese Arbeit 
bei den komplizierteren Gleichungen allerdings ziemlich 
miihsam. In manchen Fallen kann man die Rechnungen 
durch Einfiihren eines geeigneten Parameters stark ver- 
einfachen; diese Falle werden im folgenden genauer 
untersucht. ¥; 

Dabei wird zunachst nur die Abhangigkeit der 
kritischen Dynamen vom Schlankheitsgrad untersucht 
und graphisch dargestellt; die Abb. 9, 10, 11, 12 und 13 
beziehen sich alle auf Federn mit kreisférmigem Draht- 
querschnitt (2.4). Die Abhangigkeit von der Quer- 
schnittsform wird in Abschnitt 4d diskutiert werden. 

Ferner sind die Knickkurven nur fiir das praktisch 
allein wichtige Intervall 0 <€& < 1 gezeichnet. Selbst- 
verstandlich kénnen aber, falls nétig, mit den Knick- 
gleichungen auch die Knickkurven fiir € < 0,_also fiir 0 05 7 
Federn unter Zug und Torsion berechnet werden. Abb. 8. Die Knickkurven fiir reine pseudotangentiale 

Haringx+ hat ‘gezeigt, daB sich auch der Bereich Torsionsbeanspruchung im Knickfall III. 

€ => 1 realisieren ]aBt, indem man die Feder umstiilpt, 
d.h. den ersten Ring festhalt, den zweiten unter leichter Deformation iiber den ersten zieht, dann 
den dritten tiber die ersten beiden usw. Haringx zeigte fiir reine Druckbeanspruchung, dab die 
Knickgleichungen mindestens qualitativ die tatsachlichen Verhaltnisse auch hier richtig wieder- 
geben; und auch wir kénnten fiir kombinierte Druck- und Torsionsbeanspruchung die Knick- 
kurven chne weiteres in diesen Bereich => 1 fortsetzen. 

Knickfall I. Um die Knickkurven Punkt fiir Punkt auszuwerten, wahlen wir die GréBe q 
(2.27) als Parameter. Die Knickgleichung (I) lautet dann nach p aufgelést 


a 


a & 


O5 


AHR 


__ 2q(cos w—cosq) 
ame sin q 2 (4.14) 
und aus (2.27) ergibt sich 
pUsep)=P—v. (4.15) 
Dividiert man diese beiden letzten Gleichungen durcheinander, so kommt man auf 
(q? — w*) sin q 
ee ~ 2q(cosw—cosq)’ (tp) 
wahrend nach (4.7) 
1—bé 
: Se 4. 
(2 ESD Toot (4.17) 


1 J. A. Haringx, Appl. Sci. Res. A 1, (1949), S. 417. 
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ist. Mit (4.16) und (4.17) berechnet man leicht € in Funktion von w und q und daraus laBt sich 
mit (4.7) und (4.14) auch der Schlankheitsgrad 


a BEAL 4.18 
| aene Ce 
in Funktion von w und q darstellen. Variiert man jetzt q bei festem w, so erhalt man die Knick- 
kurven in einem (s, £)-Diagramm mit w als Scharparameter, wie sie Abb. 9 — die der Arbeit von 


Haringx! entnommen ist — fiir Kreisdrahtquerschnitt (2.4) zeigt. 

An Stelle von w kann gemaB (4.2) auch der — in der Abb. 9 auch angegebene — totale Ver- 
drehungswinkel @ als Scharparameter verwendet werden. Die Kurven sind so ausgewahlt, daf 
nicht w, sondern w besondere Werte annimmt; speziell erhalt man in allen Knickfallen fir w=w=0 
die Knickkurve fiir reine Drucklast (4.8). 


E 


Abb. 9. Die Knickkurven im Knickfall I. 
e 


Die Knickkurven zerfallen in zwei verschiedene Typen, die analog auch in den andern Knick- 
fallen auftreten und mit Kurven erster und zweiter Art bezeichnet werden sollen; je nachdem, ob 
w unterhalb oder zwischen den fiir reine Torsionsbeanspruchung méglichen Grenzknickmomenten 
Wo min UNA Wo max (Labelle 2) liegt. Die beiden Scharen werden getrennt durch die Kurve w = 
Womin = 7%, die selbst keiner Schar angehért, sondern als einzige Kurve die €-Achse in einem end- 
lichen Punkte schneidet. (Fiir w > wo,,,, Wirden noch Kurven einer dritten Art auftreten, die 
aber in allen Abbildungen weegfallen, weil sie vollstandig unterhalb der s-Achse verlaufen, in Uber- 
einstimmung mit der Tatsache, dafS eine solche Feder nur durch Zug stabilisiert werden kann.) 
Fir w zwischen wy,,;,, = 7% und Wong, = 4,494 besitzen die Knickkurven (zweiter Art) ein Maxi- 
mum, es existiert fiir jedes solches w ein bestimmter Schlankheitsgrad mit maximaler kritischer 
negativer Federdehnung, fiir wachsendes s strebt € gegen Null, fiir abnehmendes gegen — oo. 
Liegt w unterhalb az, (Kurven erster Art), so fallen die Knickkurven im allgemeinen monoton. 
(Dies gilt allerdings nur fiir € < 1.) Da sie die €-Achse nie schneiden, gibt es zu jedem solchen Tor- 
sionsmoment einen Schlankheitsgrad, unterhalb dessen die Feder fiir jede Druckkraft knickfest ist. 
Die Eindeutigkeit der Funktionen €(s) ist fiir € < 1 nur gestért bei sehr kleinen Torsionsmomenten 
und s in der unmittelbaren Umgebung von 5,38, wo es zu jedem s-Wert zwei oder drei €-Werte gibt. 
VergréBert man fiir festes s den zweiten &-Wert etwas, so nehmen auch die w zu, d. h. eine knickende 
Feder kann in diesem kleinen Bereich — und nur hier — bei fester Druckkraft durch VergréBern 
des Torsionsmomentes und bei festem Torsionsmoment durch VergréBern der Druckkraft wieder 
stabilisiert werden. Der Bereich ist allerdings so klein, daB sich dieses Verhalten fiir Kreisquer- 
schnitte kaum feststellen lassen diirfte (Vgl. Abschnitt 4d). 

Knickfall II. Die Tatsache, daB fiir eine Feder mit beidseitig gelenkig gelagerten Enden und 
einer aus Druck und axialer Torsion zusammengesetzten Belastung das — allerdings nicht legitime 


1 FuBnote 5 von Seite 113. 
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— statische Stabilitatskriterium auf die Knickgleichung IT fiihrt, hatte zur Folge, daB diese Knick- 
gleichung schon friiher eingehend diskutiert worden ist. So findet man das Knickdiagramm Abb. 10 
auch bei Haringx, wahrend bei Ziegler und Huber! oder Biezeno und Grammel? die Diskussion auch 
fiir gewisse Rechtecksquerschnitte durchgefiihrt ist. Es ist dort nur zu beachten, daB die Knick- 
kurven in einem (&, 7)-Koordinatensystem dargestellt sind, wobei € wie hier (4.3) und 7 als 7 = 
w/w, = w/z definiert ist und als Scharparameter der Schlankheitsgrad s verwendet wird, der dort 
aber als /,/2 7 R und nicht wie hier und bei Haringx als |,/2 R definiert ist. 
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Abb. 10. Die Knickkurven im Knickfall II. 
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Abb. 11. Die Knickkurven im Knickfall IVs. 


Die in der Abb. 10 gezeichneten Kurven sind alle zueinander affin; es handelt sich um die 
Kurven erster Art, wahrend die Kurven zweiter Art hier gar nicht auftreten, da im Falle IT wy,,;,, = 
Wo max ist. Im Gegensatz zum Falle I besteht die Méglichkeit, die Feder durch Erhéhung des Tor- 
sionsmomentes zu stabilisieren fiir alle s < 2,69 und die zugehérigen Torsionsmomente umfassen 
alle Werte w < x. 

Knickfall IIIs. Der Fall IIIs erfordert kein eigenes Knickdiagramm, da sich die Knick- 
gleichung (IIIs) durch die Substitution wy, = 1/2 wiz, sis = 1/2 sy auf die Knickgleichung (IT) 


1 H. Ziegler u. A. Huber, Z. angew. Math. Phys. 1 (1950) S. 189. 
2 FuGBnote 1 von Seite 114. 
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zuriickfiihren laBt. Man hat in der Abb. 10 nur die Abszisse und den Scharparameter gemaB obiger 
Substitution zu verindern; das Verhalten der Feder ist analog zum Verhalten im Fall I. 
Knickfall IVs. Wie im Falle I folgt aus der Knickgleichung (IVs) 
(w? + q?) sin q — 2 w q sin w 
x q cos q 


P 


mit (4.15) und (4.17) 
1—bé (q? — w*) q cos q 7 
ae (w? + q2) sing—2wqsinw’ 


woraus sich € und mit (4.18) s durch w und q ausdriicken lassen. Die Knickkurven in Abb. 11 ver- 
laufen ahnlich wie im FalleI; es sind beide Arten vorhanden, die durch die Kurve w = W pin = 7/2 
getrennt werden. Im Gegensatz zu Fall I sind aber die Funktionen &(s) im Bereich & < | fiir alle 
s eindeutig, so daB in diesem Knickfall die VergréBerung der kritischen Druckkraft durch Erhéhung 
des Torsionsmomentes fiir Kreisdrahtfedern tiberhaupt nicht méglich ist. 


g 


Abb. 12. Die Knickkurven im Knickfall Vs. 


Knickfall Vs-s. Gleich wie bisher erhalt man aus der Knickgleichung (Vs-s) mit dem Para- 
meter q die Beziehung 
I= be) ) | (gq? — 2) q sing 
1—é& 2 w? (cos w— cos q) ” 


woraus man & und mit (4.18) s berechnet. Da in diesem Knickfall wy,,,;,, = 0 ist, sind alle Kurven 
der Abb. 12 (auBer w = 0) Kurven zweiter Art. Die Funktionen é(s) sind iiberall eindeutig, auBer 
fiir sehr kleine w und fiir s in der unmittelbaren Umgebung von 2,65, wo es wie im Falle I zu jedem 
s zwei oder drei €-Werte gibt. Wieder wiirde sich hier und nur hier die Druckstabilitat durch Ver- 
gréBerung des Torsionsmomentes erhéhen lassen; doch ist auch hier das Gebiet so klein, daB sich 
diese Tatsache praktisch kaum feststellen lassen diirfte. 

Knickfall lg. Unter Verwendung des Parameters q ]aBt sich die Knickgleichung (IIIq) als 


w=——+qctgq 


schreiben, wahrend mit (4.15) 
p (1 + €p) = ¢@° (1 — ctg? q) 

wird. Mit Hilfe der zweiten Bedingung (4.7) lassen sich aus diesen beiden Gleichungen fiir festes s 
zu jedem q die Groen w und & berechnen. Dieser Fall mu® daher anders als die bisherigen darge- 
stellt werden, es lassen sich die Knickkurven nur in einem (6, w)-Koordinatensystem mit s als 
Scharparameter angeben. Abb. 13 zeigt diese Knickkurven fiir Kreisdrahtfedern. Selbstverstand- 
lich sind die Knickkurven symmetrisch in bezug auf die -Achse, da es wie in allen andern Knick- 
fallen auf das Vorzeichen von w nicht ankommt. Zu jedem Torsionsmoment w < m/4 gibt es auch 
hier einen Schlankheitsgrad, unterhalb dessen die Feder fiir jede Druckkraft knickfest ist. Im Be- 
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reich 0,8125 > & fallen die Knickkurven monoton, wahrend sie, falls tiberhaupt vorhanden, im 
Bereich € > 0,8125 wieder ansteigen. In diesem Gebiet ist also wieder Stabilisierung durch Er- 
héhung des Torsionsmomentes méglich und zwar nur fiir Schlankheitsgrade s < 1,346. 

Vergleicht man schlieBlich die verschiedenen untersuchten Knickfille noch miteinander, so 
stellt man fest, daB der Satz sich auch hier bestatigt, wonach eine Vermehrung der Bindungen eine 
Erhéhung der kritischen Belastung bewirkt. So schneiden sich z. B. entsprechende Kurven der 
Falle I und II oder Vs-s und IVs nicht, dagegen diirfen verschiedenartige Bindungen nicht mit- 
einander verglichen werden, wie z. B. die Schnittpunkte der beiden Kurven w = 1,81 in den Knick- 
fallen IT und IVs zeigen. 


d) Andere Drahtquerschnitte. Um die Abhangigkeit der kritischen Lasten der Feder von 
der Querschnittsform zu diskutieren, kann man die Querschnittsparameter a und } in den Formeln 
(4.7) beim Berechnen der Punkte der Knickkurven durch andere Werte als (2.4) ersetzen und die 
neuen Knickkurven wieder in einem (s, £)-Koordinatensystem darstellen. Einfacher ist es jedoch, 


TU 
a8 - 
06 
Oo 
G2 
S=302\ \2,593 
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Abb. 13. Die Knickkurven im Knickfall Iilq. 


die alten Knickkurven fiir den Kreisquerschnitt beizubehalten und dafiir ein krummliniges Koor- 
dinatennetz zu verwenden. Die Tatsache, dafS in den Knickgleichungen nur p und p (1 + ¢ p), 
nicht aber w von a und b abhangen, bewirkt, daB die Knickgleichungen erfiillt bleiben, solange nur 
diese beiden Ausdriicke ihren Wert nicht andern. Bezeichnet man die den KreisquerschnittsgréBen 
a = 23/26, b = 8/13, s und & beim neuen Profil entsprechenden GréBen mit a, b, s und €, so hat 
man gemaf (4.7) nur zu verlangen, daB 


p=4as(l—é)& =4a8(1—£)E 
und 
p(l+ep) =4a82(1—béE=—448(1—b HE 


ist. Daraus folgen durch Auflésen nach s und E die Beziehungen 


(L —b)é_ a pl bet Lop (4.19) 
(eye aes a (1 —b) (1 —b) 


Um die Knickwerte des neuen Profils zu bekommen, hat man jetzt nur das Kurvennetz s = kon- 


Wy | 


stant und é = konstant im (s, £)-Koordinatensystem iiber die alten Knickkurven zu legen und die 
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neuen Werte s und € aus diesem Netz abzulesen. Der Scharparameter w bleibt dabei ae 
dagegen muf der totale Verdrehungswinkel « mit (4.2) fiir jede Kurve aus w neu berechnet ti en, 
weil die GréBe 2 x»/d) nach (2.2) von der Querschnittsform abhangig ist. Natiirlich kénnen mit ieser 
Methode nur die Falle behandelt werden, bei denen die Knickkurven schon in einem (s, &)-Koordi- | 
natensystem vorliegen, also z. B. nicht der Fall IIIq. 
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Abb, 14. Das Koordinatennetz fir einen rechteckigen Drahtquerschnitt mit dem Seitenverhaltnis y = 1/2. 
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Abb. 15. Das Koordinatennetz fiir einen rechteckigen Drahtquerschnitt mit dem Seitenverhiltnis p=2 


Abb. 14 bzw. 15 zeigen diese Koordinatennetze fiir Rechteckprofile mit den Seitenverhaltnissen 


pe = 1/2, bzw. w = 2, fiir die man a und 6 der Tabelle 1 entnehmen kann. Die Kurven € = konstant 
sind immer zur s-Achse parallele Geraden, die zwar nicht aequidistant sind, aber das interessante 
Intervall 0 < € <1 als Ganzes nicht verandern. Die Kurven s — konstant sind alle zucinander 
affin und in bezug auf die &-Achse fallend, bzw. steigend, je nachdem 6 > 8/13 oder 6 < 8/13 ist. 

Denkt man sich nun diese Netze in Abb. 9 bis 12 eingetragen, so erkennt man, da®B zwar die 
einzelnen Knickwerte numerisch stark mit der Querschnittsform variieren kénnen, da sich aber 
am ganzen Diagramm qualitativ nur wenig andert. So zerfallen die Knickkurven weiterhin in die 
zwei Typen, die beide ihre charakteristischen Eigenschaften beibehalten. 
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Die GréBe des Bereichs, in dem eine Erhéhung des Torsionsmomentes auch eine VergréRerung 
der kritischen Druckkraft bewirkt, hangt allerdings stark von der Querschnittsform ab. Unter- 
sucht man den Bereich, wo die Kurven s = konstant die Knickkurven in mehr als einem Punkte 
schneiden, so stellt man fiir alle im (s, €)-Koordinatensystem dargestellten Knickfalle fest, daf 
dieser Bereich fiir ~ = 2 (Abb. 15) ganz verschwunden ist, da er aber fiir 4. = 1/2 (Abb. 14) ge- 
geniiber dem Kreisprofil sich stark vergréBert hat, (ja im Knickfall [Vs tiberhaupt erst auftritt), so 
da8 er fiir ein solches Rechteckprofil auch praktisch eine Rolle spielen diirfte. 

Allgemein ergibt sich daraus, da in allen diskutierten Knickfallen bei einem in Richtung der 
Schraubenachse geniigend langgestreckten Rechtecksprofil eine Erhéhung des Torsionsmomentes 
immer eine Verminderung der kritischen Druckkraft bewirkt und umgekehrt, wahrend bei Federn 
mit hochkant gewickeltem Rechtecksprofil immer ein Bereich fiir den Schlankheitsgrad existiert, 


wo eine Erhéhung der Druckkraft eine Erhéhung des kritischen Torsionsmomentes und umgekehrt 
nach sich zieht. 


(Eingegangen am 26. Juni 1956) 
Anschrift des Verfassers: H. Maag, dipl. Math. ETH. Bern (Schweiz), Simonstr. 25 
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Ein Sonderfall des senkrechten Profilgitters bei beliebigen 
instationdaren Bewegungen 


Von K. Nickel 


1. Einleitung. In einer kiirzlich hier erschienenen Arbeit! wurde u. a. das senkrechte Fligelgitter 
(Abb. 1) bei instationéren Bewegungen behandelt. Das strémende Medium war als volumenbe- 
standig und reibungsfrei angenommen, die Strémung sollte in der Ebene erfolgen, ihr instationares 
Verhalten war an allen Gitterschaufeln synchron und durch kleine Stérungen hervorgerufen. Unter 
diesen Voraussetzungen lief sich das Problem auf die Lésung von zwei Integralgleichungen zuriick- 
fiihren. Die erste, eine Fredholmsche Integralgleichung erster Art, beschreibt analog zur ,, Birn- 
baumschen Integralgleichung’® fiir das Hinzelprofil den stationdéren Anteil der Strémung und 
konnte vor einiger Zeit gelést werden®. Die instationare Komponente der Strémung, die sich physi- 
kalisch in der abschwimmenden ,,Wirbelschleppe“ (Abb. 1) bemerkbar macht, gehorcht mathema- 
tisch einer Volterraschen Integralgleichung erster Art entsprechend zur ,,Wagnerschen Integral- 
gleichung“* des Einzelprofils. Die Lésung dieser Integralgleichung wurde vor kurzem fiir harmo- 


nische Schwingungen von Ch.-Ch. Chang und W. H. Chu® angegeben; ferner vergleiche man dazu 
H. Séhngen®. DanachlaBt sich derlésende 


Kern der Wagnerschen Integralgleichung 
S in diesem Fall aus der hypergeometri- 
schen Funktion und der Gammafunktion 


ce sere —*) aufbauen und geht beim Ubergang zum 
Te J Einzelprofil (t > oc) in den entsprechen- 
on Wirbelsoemen jen Ausdruck aus Hankelfunktionen 


4 Antahrwirbel, 
Anstrimung Ea Saree =) SS (siehe H. G. KiiBner? und H. Séhngen®) 
ca D a J uber. Der allgemeine Fall der beliebigen 


Bewegungen laBt sich daraus als Fourier- 
integral gewinnen. 


Ce ee fe : 

6 > In der folgenden Note soll gezeigt 
° werden, da sich die mathematische L6- 
° sung des betrachteten Problems erheb- 
Abb. 1. Senkrechtes Profilgitter. lich vereinfacht, wenn man zu dem Son- 


derfall unendlich tiefer Gitterschaufeln 
(l—> co) iibergeht (Abb. 2), bei dem also die Profilvorderkanten in Anstrémrichtung ins Unendliche 
geriickt sind. Man kommt auf diesen Fall, wenn man bei einem ,,nichtentarteten“ Gitter (J endlich) 
gemaB Abb. 1 die Betrachtung auf eine kleine Umgebung der Hinterkante beschrankt, wobei auBer- 
dem noch die instationare Bewegung erst kurze Zeit andauert (s < /in Abb. 1). Dieser Grenziiber- 
gang 1+ oo bei fester Gitterteilung t bedeutet rein geometrisch dasselbe wie der Séhngensche® 
Grenziibergang t-—> 0 bei konstanter Profiltiefe 1. Séhngens Betrachtungsweise hat den Vorzug, 
daB er bei der Ubertragung seiner Rechnungen auf endliche Werte des Teilungsverhaltnisses t/I 
nicht wie hier nur auf den ersten Augenblick der instationaren Bewegung angewiesen ist, und da- 
her seine Rechnungen z. B. fiir Gitterschwingungen nutzbar machen kann. Als Nachteil muf da- 
fiir der Naherungscharakter der Lésungen in Kauf genommen werden. Im Gegensatz dazu sind 
die vorliegenden Untersuchungen dem Bewegungsbeginn angepaBt. Damit bleibt zwar die An- 
wendung auf das Gitter der Abb. 1 nur auf kleine Zeitintervalle beschrankt, jedoch werden sich 
dafiir geschlossene Formeln aus elementaren Funktionen ergeben, so daf es sich doch lohnt, auch 
diesen Sonderfall zu beleuchten. 


1 K. Nickel, Ing.-Archiv. 23 (1955) S. 179. 

2 W. Birnbaum, Z. angew. Math. Mech. 3 (1923) S. 290. 

* J. Dorr, Ing.-Arch. 19 (1951), S. 66 und K. Nickel, Z. angew. Math. Mech. 31 (1951) S. 285. 
* H. Wagner, Z. angew. Math. Mech. 5 (1925) S. 17. 

° Ch.-Ch, Chang u. W. H. Chu, Journ. Appl. Mech. 22 (1955) S. 503. 

® H. Séhngen, Z, angew. Math. Mech. 35 (1955) S. 81. 

’ H.G. Kuffner, Luftfahrtforschung 17 (1940) S. 355. 

8 H. Sohngen, Luftfahrtforschung 17 (1940) S. 401. 

° H. Sohngen, Z. angew. Math. Phys. 4 (1953) S. 267; vgl. auch FuBnote °, 
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2. Mathematische Formulierung des Problems. In Abb. 2 ist das Koordinatensystem so gelegt 
worden, daf die x«-Achse in Verlangerung eines Gitterprofils liegt, mit dem Nullpunkt in der Schau- 
felhinterkante. Da das Gitter der Abb. 2 keine bevorzugte Linge mehr besitzt, kann die Teilung ¢ 
beliebig gesetzt werden, es geniigt also, etwa die beiden Falle t —27 (Profilgitter) und t— co 
(EKinzelprofil) zu betrachten. Der andere mégliche Grenzfall t—> 0 besitzt kaum praktisches Inter- 
esse und soll daher im folgenden nicht weiter verfolgt werden. Die Strecke s ist ein Ma® fiir die 
Zeit, die seit dem Beginn der instationéren Bewegung verstrichen ist. Die Normalkomponente 
w(x, s) der Strémungsgeschwindigkeit an 
einem Fliigelpunkt auf der Strecke —oco 
<x <0 bestimmt sich aus der Forderung, 
daB der Fliigel nicht durchstrémt werden 
darf, und kann daher im folgenden als 
bekannt angesehen werden. Ersetzt man, 
wie in Abb. 2 angedeutet, die Gitterprofile 
durch gebundene Wirbel der Starke /"(x, s) 
auf der negativen x-Achse und entsprechend 
die Wirbelschleppen durch freie Wirbel der 
Stirke e(x,s) auf der positiven x-Achse}, 
so bestimmen sich die beiden Funktionen /"(x, s) und ¢(x, s) aus der Gleichung 


Abb. 2. Senkrechtes Profilgitter fiir 1—+ co. 


0 s 
i Rees i ies(Sors ms 
we) = 52 Gerd ty5 | Tiedt fir too (1a) 
—oo 0 
bzw. 
0 ; s 
1 =e en eG ; 
wees) =, ETE) C0*s dB +7; | os) Cig Tue fir t= 2e, (1b) 
—oo 0 
Weiter sind noch die folgenden Relationen zu beachten?: 
0 8 
J TGs) db + f elG, 3) de =0, (2) 
IO, 8) = £(0, 8), (3) 
E(x, 8) = £(0, s— x) (Zn) (4) 
SchlieBlich werde noch durch 
W(x, 8) = I(x, s\-}- = f L(&, s) dé (x < 0) (5) 


die értliche Auftriebsverteilung (x, s) am Profil eingefihrt. 
Berechnet man analog wie beim endlichen Gitter? /’ aus den Gleichungen (1) bis (4) und setat 


in (5) ein, so findet man die Auftriebsverteilung tiber die Fliigeltiefe 


0 
dy caer H(é, s) d fir t—> 6 

ema)? eof ms Sian i 

0 

cea H€, s) ? 
x, 8) = —j1—e’ —— d§ -fir t=2a (6b) 
y( ) mali Ga aia 

mit der Abkiirzung | 

H(x, 8) = w(x, 8) +2 f wlé, s) dé. (7) 


Die Formeln (6) haben sich gegeniiber den entsprechenden Gleichungen des Gitters aus endlich 
langen Schaufeln® sehr vereinfacht, insbesondere kommt in (6) die Starke «(&, s) der freien Wirbel 
1 Die genaueren Voraussetzungen dieser Setzung sind in der o.a. Arbeit (FuBnote! von Seite 134) angegeben. 


2 Vel. FuBnote | von Seite 134. 
3 Vel. FuBnote ! von Seite 134, Formel (28). 
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gar nicht mehr vor. Dieser merkwiirdige Sachverhalt bedeutet physikalisch, daB die drtliche Auf- 
triebsverteilung am Gitterprofil bei dem betrachteten speziellen Schaufelgitter nur vom augen- 
perme es instationdren Bewegungszustand abhangt, nicht aber von der Woreeschigiee der Be- 
vegung, wie sie sich in den abschwimmenden Anfahrwirbeln ausdriickt. Die Formeln fiir die Auf- 
vieleat envoiltite tragen also quasistationaren Charakter, nur ist gut zu beachten, daB das im sta- 
tiondren Fall allein rauftrevents Glied w(x) hier nicht nur durch ae s), sondern durch H(x, s) nach 
(7) zu ersetzen ist. 
Fir den Auftrieb eines Schaufelstiicks a < x <0 ergeben sich die geschlossenen Formeln 
(0 = Dichte des stromenden Mediums, u.. = Anblasgeschwindigkeit) 


maleic 


0 ae 
SHG, s) | == TA dx fir t—oo, (8a) 
A= 0 Uo [ y(%, 8) dx =— 0 Ucoy 9 ee : B ; 
a : ; aretg | «—@ ales ae e—* —1 + /e— a2 q 
Ae (ei 108 ——— 1 


fir t=27. (8b) 


Im Falle des Einzelprofils (¢ > oo) fiihrt auch die Berechnung des Moments dieses Schaufelstiicks 
um einen beliebigen Drehpunkt 6 (linksdrehend positiv) noch auf einen geschlossenen Ausdruck 


M=o eo J y(x, s) (x — b) dx 
s) pe [: b+ a (x — b) log a 


Als Beispiel werde die geknickte Platte mit dem Ruderausschlagwinkel 7 bei stationdrer 
Strémung berechnet (Abb. 3). Mit 


x = {0 fur)” “e-<a 
w(M5.S hae I (a, 8) == lites far axes 5) und alle s 


Abb. 3. Flosse mit Ruder. findet man 


dx fiir t—-oo. (9) 


A 0 
=) les if LMG 


4 ae ogy ak a ait 4 P 4 

Kn) == Uo | ages : se es mer cen/ Tie hc M =—1 Q Ue a b— Fa). 
Dieses Ergebnis ist nach den aiieeh sre: auch z. B. fiir die plétzlich einsetzende Bewe- 
gung, den ,,Anfahrvorgang™ giiltig. 


1 1 

3. Die Wagnersche Integralgleichung. Multipliziert man (la) mit (— x) ®, (1b) mit (c= 1) 23 
integriert iiber (—co,0) und beachtet (4), so kommt die Wagnersche Integralgleichung des 
Problems! 


s 0 
" (0, 0) w(x, s) an 
| =e 12 ee dc — 0% fir t->o, (9a) 
: , 
| es do 4 2 | = dx = 0 firs h— 40 (9b) 
= eo—s ex ek 


Gleichung (9a) besitzt schon die Gestalt der wohlbekannten Abelschen Integralgleichung?, aber 
auch (9b) wird durch die Setzung 


S=e—1, J=e—1, «0,c) =(14+ 5) 6(0,), w(x,s) = + S w(x, S) 


in die Integralgleichung 


i ‘Gant 1 von Seite 134. 
2 Vel. FuBnote 1 von S. 137. 
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von demselben Typ ibergefiihrt. Mit der bekannten Lésung der Abelschen Integralgleichung? 1aBt 
sich dann sofort die Starke ¢ der freien Wirbel in geschlossener Form angeben zu 


s 0 
14 Be Gh ff il w(x, 0) & 
€(0, s) oars | ae ae dx do fir t->oo, (10a) 
0 —oo 
ed s 0 
cf ! A 1 * w(x, o 
(0,9) =— = £ ea — dx do. fir 4=2n. (10b) 
0 —oo 


4. Beispiele: Es gelte 
ae nO fiir c= 0 
thts 8) > w(x) qs) fir s>0. 
Dann ergibt sich folgender Zusammenhang zwischen einigen einfachen Funktionen q(s) und den 
zugehérigen Wirbelstarken ¢(0, s): 


Oe | Pages 
ele le | 4 
eee q(s) | 1 | — | el? —] ite 
| eee 
Dabei wurde zur Abkiirzung gesetzt 
22 EC eral cae ec 


5. Folgerungen aus den abgeleiteten Formeln. a) Ein besonders einfacher Fall einer instationaren 
Bewegung ist der Anfahrvorgang, wenn also das Fliigelgitter aus der Ruhe heraus plotzlich mit 
einer konstanten Geschwindigkeit u.. bewegt wird. Obwohl der Fliigel seine Lage nach dem An- 


{ fr (8) 


40 


G5 


0 20 40 60 80 100 720 140 
Abb. 4, Abminderungsfaktor k,(s) fiir den Auftrieb des endlich langen Einzelprofils beim Anfahrvorgang. 


fahren nicht mehr andert, stellt sich beim endlich langen Einzelfliigel oder beim endlichen senk- 
rechten Fligelgitter der endgiiltige Auftrieb nicht sofort ein, sondern wird erst asymptotisch 
erreicht. Man erhalt die Auftriebsverteilung und damit auch den Auftrieb aus dem stationaren 
Wert durch Multiplikation mit einem Abminderungsfaktor? k,(s) (Abb.4). Da man den hier voll- 
zogenen Grenziibergang zum unendlich langen Fliigel (J > oo) auch als Lésung fiir den ersten Beginn 
der Bewegung (s < /) ansehen kann, folgt schon rein anschaulich, daB beim Anfahrvorgang jetzt 
k,(s) = k,(0) = konst. wird (Abb. 4). Dieses Ergebnis bestatigt man sofort aus (6) und (8), da ja 
beim Anfahrvorgang 0 H (x, s)/ds = 0 wird fiir s > 0. Ganz entsprechend gilt offensichtlich auch 


1 Vel. W. Schmeidler, Integralgleichungen mit Anwendungen in Physik und Technik. Seite 212/213, Leip- 


zig 1950. 
2 Vel. FuBnote ! von Seite 134 und FuRnote ® von Seite 134. 
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die Verallgemeinerung: Hat sich das betrachtete Fligelgitter in beliebiger Weise bewegt und endet 
diese instationére Bewegung plétzlich, so haben im gleichen Augenblick Auftriebsverteilung und 
Auftrieb schon ihren endgiiltigen Wert angenommen, namlich denjenigen des stationaren Zustands. 

b) Da der Auftrieb eines Fliigelprofils mit dem Impuls der hinter dem Fliigel abwarts bewegten 
Strémung zusammenhdngt, vermutet man weiterhin: Wenn sich der Fliigelauftrieb nicht mehr 
aindert, so induziert der Fliigel auch keine weiteren Anderungen in den Abwindgeschwindigkeiten 
hinter dem Fliigel, obwohl die Wirbelschleppe evtl. noch nicht abbricht, also noch weitere freie 
Wirbel abschwimmen. Dieses Ergebnis laBt sich in der Tat aus den angegebenen Gleichungen 
erhalten. Die dazu notwendigen Umformungen sind zwar einfach, jedoch etwas umfangreich und 
mégen daher iibergangen werden. Es gilt also: Hat die instationdre Bewegung des betrachteten 
senkrechten Fliigelgitters zum Zeitpunkt s, aufgehért, so hangt die induzierte Abwindgeschwindig- 
keit w(x, s) in Verlangerung der Fliigelprofile im Bereich 0 < « < s, nicht mehr von s ab. Dieses 
Ergebnis gilt nicht mehr auBerhalb der Fliigelsehnen, wie ein Blick auf das Beispiel in Abb. 5 zeigt. 


Abb. 5. Anfahrvorgang am senkrechten Profilgitter. Profiltiefe ]—> co, Rudertiefe 1p = a, Teilung t = 27, s = 4 
Eingetragen sind die Abwindgeschwindigkeiten auf der Profilsehne und in der Mitte zwischen zwei Profilen, 


c) Fir das endlich lange Einzelprofil gilt bekanntlich das interessante Ergebnis, daB Fliigel- 
flattern schon bei nur einem Freiheitsgrad méglich ist. Z. B. gibt es nach Séhngen' eine reine 
Ruderschwingung eines reibungs- und kraftefrei angelenkten Ruders. Séhngen bemerkt dort, daB 
allein fiir unendlich kleines Ruder (oder: unendliche Flossenlange bei endlichem Ruder, also in dem 
hier vorliegenden Sonderfall), keine Anfachung mehr vorhanden ist. Setzt man namlich (Abb. 3) 


| " { 0 fir —coo<x<a, 
ro(% 8) = Hoo | n(s) + («—a)y'(s) fir, a Sx SO, 
so wird 


0 fiir —09 <*% <4, 


EL (as )\ezesUhas 1 

n(s) + 2(% — a) 7'(s) + a —— a) ni (Ss) fin ars oe Oe 
Wegen der Linearitaét von (9) hat dann das Rudermoment um den Drehpunkt x = a die Gestalt 
M=4Ay7+ By’ + Cx" mit Konstanten A, B und C. Wirken nun auf das Ruder auBer diesem 
Luftkraftmoment M noch andere Momente ein, die nur von 7 (wie das Riickstellmoment) und 7’ 
(Tragheit) abhangen, so sind ungedampfte Schwingungen nur méglich, wenn der Koeffizient B 


1 FuBnote ® von Seite 134. 
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von 7’ verschwindet. Aus (9) ergibt sich jedoch 


pee io, 
=~ 978 oo w+ 0, 


womit Séhngens Bemerkung noch einmal bestatigt ist. Dieses Ergebnis gilt jedoch auch noch fiir 
ein senkrechtes Gitter aus festen Flossen, die sich unendlich weit in Anstrémrichtung erstrecken 
und an die reibungsfrei (keine Momente proportional zu 7’) im Punkte a Ruder angebracht sind 


(Abb. 3). Der wesentliche Term 2 u., (x — a) 7’ von H(x, s) liefert hierbei nach (6b) eine Zirku- 
lationsverteilung 


y(%, 8) = “Het 
nt ne E—x PE ae , 
2 Vi a e Ve s | 


‘ 0 rs — 
— | Bemis dé + 4 (x — a) log eats : ap eat pa 
Sin 


die in a < x <0 das Vorzeichen nicht wechselt. Damit verschwindet aber auch das zugehérige 
Glied B 7’ des Luftkraftmoments nicht und dies bedeutet ja gerade die Unméglichkeit einer unge- 
dampften Schwingung. 

d) Die angegebenen geschlossenen Ausdriicke erlauben die einfache Berechnung und das Stu- 
dium vieler Probleme, fiir die im Falle des endlich langen Gitters der Formelapparat zu unhandlich 
werden wiirde. Nur darf diese Leichtigkeit nicht dariiber hinwegtauschen, da8 eine Ubertragung 
der gewonnenen Ergebnisse auf das allgemeine Gitter der Abb. 1 haufig nicht zulassig ist, wie die 
vorangegangenen Beispiele a) bis c) zeigen! 


(Eingegangen am 6. Juli 1956.) 


Anschrift des Verfassers: Dr. Karl Nickel, Karlsruhe, Institut fiir angewandte Mathematik der Technischen 
Hochschule. 
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Die Extrema der Geschwindigkeiten an Kurbeltrieben 


Von W. Meyer zur Capellen 


1. Einleitung. Bei Kurbeltrieben, d. h. den Formen der Viergelenkkette, die Kurbeln enthalten, 
also bei der Kurbelschwinge (Abb. 1) und bei der Doppelkurbel (Abb. 2) nebst ihren Sonderformen, 
wird immer wieder die Frage nach der gréBten Abtriebsgeschwindigkeit bei gleichférmigem Antrieb 
gestellt. Fiir spezielle Falle, wie fiir die zentrische! und fiir die geschrankte Schubkurbel ist die 
Frage zusammenfassend auf rechnerischem Wege beantwortet worden”, wobei die Rechenarbeit 
im letzteren Fall nicht geringfiigig war. Es soll nun hier nach den allgemeinen Bedingungen gefragt 
werden, die erfiillt sein miissen, damit das Abtriebsglied, d. h. also die Schwinge baw. die Abtriebs- 
kurbel, oder aber auch die Koppelebene eine extremale Winkelgeschwindigkeit im Verlaufe eines 
Bewegungsspieles hat. Hierzu miissen vorher kurz einige im 
Folgenden bendtigte, allgemeine kinematische Beziehungen 
zusammengestellt werden. 


2. Benutzte Ortskurven. a) Die Bresseschen Kreise. 
Alle Punkte einer komplan bewegten Ebene, welche im 
Augenblick einen Wende- oder Undulationspunkt ihrer Bahn 
durchlaufen und damit nur Tangential-, aber keine Normal- 
beschleunigung haben, erfiillen den Wendekreis K,,, der die 


Abb. 1. Kurbelschwinge, 
Abybeoiy Ieee loy Alea, Zhe aah 


aN 


g 


Abb. 2. Doppelkurbel. Abb. 3. Wendekreis und Gleichenkreis. 


Polbahntangente t im Momentanpol P beriihrt; sein Durchmesser sei PW = D (Abb. 3). Dem- 
gegentiber liegen alle Punkte, welche im Augenblick nur Normal-, also keine Tangentialgeschleu- 
nigung haben, auf dem Gleichenkreis K,, welcher die Polbahnnormale n in P beriihrt; sein Durch- 
messer sei PG = g (Abb. 3). Beide Kreise schneiden sich (unter rechtem Winkel) im Beschleu- 
nigungspol J. Wahrend der Wendekreis durch die Geometrie des Getriebes festgelegt ist, bestimmt 
der Beschleunigungszustand die Lage von G: Aus Winkelgeschwindigkeit w und Winkelbeschleu- 
nigung € der bewegten Ebene folgt 


2 
g= PE =D— =Detgs=Dtgqy, (1) 


wenn Winkel PGW = 06 und Winkel JPG = gy gesetzt ist. Die Beschleunigungsvektoren bilden 
mit den jeweilig zum Beschleunigungspol J gezogenen Strahlen den konstanten Winkel 6 =2/2— qo, 
so daB z. B. die Polbeschleunigung von P nach W gerichtet ist. 


1 W. Meyer zur Capellen, Getriebetechnik 5 (1937) S. 529. 

* W. Meyer zur Capellen u. Mitarbeiter, Bewegungsverhaltnisse an der geschrankten Schubkurbel. Er- 
scheint demnachst in den vom Wirtschafts- u. Verkehrsministerium von Nordrhein- Westfalen herausgegebenen 
Forschungsberichten. 
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Durch den Wendekreis sind auch die Kriimmungsverhaltnisse der bewegten Ebene vermége 
der Euler-Savaryschen Formel 


Sat ir (2) 


gegeben, worin (Abb. 4) r = PA, ry) = PA,, w = PAy = D sin y und A, der Kriimmungsmittel- 
punkt der von 4 beschriebenen Bahnstelle ist. 


b) Kreisungs- und Angelpunktkurve. 1. Allgemeine Darstellung. Bekanntlich erfiillen die- 
jenigen Punkte der bewegten Ebene, die im Augenblick eine Bahnstelle mit vier- oder mehrpunktig 
berithrendem Kriimmungskreis durchlaufen, die Kreisungspunktkurve k,, und die dazugehérigen 
Kriimmungsmittelpunkte die Angelpunktkurve! k,. Es handelt sich um zirkulare Kurven dritter 
Ordnung, welche in Polarkoordinaten r, p baw. ro, ~ 


Ag 

die folgende auf P als Pol und t als Polarachse be- 
zogenen Gleichungen haben: 
k,: r(mcosy + lsingy)—mlsing cosy =0, (3a) 
oder 

1 1 1 

— i 

r lsing ' mcosq~’ 8) 
ferner nach (2) 
k,: r9(mcosy +1, sing~)—ml sing cosg = 0, (4a) i 
oder 

tomer e*.1 

ro sing ' mcosp’ =>) 
Hierbei sind — wesentlich nach Miiller! — die fol- ji 
genden Abkiirzungen benutzt worden: 

3 sy 
Fe ay? ES re | | t 
3 Pp 
PE = IP = mee > Abb. 4, Erlauterung zur Euler-Savaryschen Formel, 


ds = Bogenelement der Rastpolbahn, di) = Elementardrehung der bewegten Ebene, d. h. d)/dt=w, 
dr = Kontingenzwinkel der Rastpolbahn, Striche Ableitungen nach s, d. h. z. B. } = dd/ds = 1/D, 
t’ = dt/ds = 1/R, = Kriimmung der Rastpolbahn. 

SchlieBlich sind / und J, 
die Durchmesser der Kriim- 
mungskreise von k, bzw. k, 
in P, deren Mittelpunkte auf 
n liegen und ist m der Durch- 
messer des Kriimmungskrei- 
ses vonk, undk, in P, dessen 
Mittelpunkt auf ¢ liegt. Zu 
beachten ist noch, da nach 
(3) einerseits 

Hl if 1 

aes oD (5a) 
und andererseits 

1 1 1 

aD (5b) 
mit R als Kriimmungsradius 
der Gangpolbahn gilt. 


c) Sonderfille. Fiir die vorliegende Untersuchung werden von den vielen an sich médglichen 
Sonderfallen nur die folgenden benétigt: 


Abb. 5. ,,Kardanlage‘‘: PH senkrecht A,B. 


itung v ii infu in di i i i lin 1931 (dort werden 
1 Herleitung vgl. z. B. Miller, R., Einfiihrung in die theoretische Kinematik, Ber ( 
k,, und See Repke und Mittelpunktkurve bezeichnet), ferner R. Beyer, Kinematische Getriebesynthese, 
u a 


Berlin/Gottingen/Heidelberg 1953. 
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a) Die Kardanlage?: Die Mittelpunktkurve k,, entartet in die Polbahntangente t und den die 
Polbahnnormale n in P berihrenden Kreis k,* vom Durchmesser PM = m, wobei k,, nicht ent- 
artet. Diese Entartung tritt auf, wenn J, = oo, d. h. 1 =D und 7’ = oO =0, dh; Ry = D, also 
nach (5b) R = D/2 wird. Der Kriimmungsradius der Rastpolbahn ist doppelt so groB als jener 
der Gangpolbahn, wobei die erstere die letztere umschlieBt, und da der Kriimmungskreis der Gang- 
polbahn mit dem Wendekreis iibereinstimmt, liegt momentan die Kardanbewegung vor, wonach 
der Name ,,Kardanlage“ geprigt wurde. Beim ebenen Gelenkviereck hat man diesen Fall in der 
Totlage oder wenn die Verbindungsgrade des Momentanpols P = 20 mit dem Relativpol H = 31 
(Abb. 5) auf dem Steg A, By senkrecht steht. 


Abb. 6. Kurbeltriebe in den Steglagen. a,b) Kurbelschwinge; c, d) Doppelkurbel 


6) Die umgekehrte Kardanlage: Bei Umkehr der soeben betrachteten Lage entartet die 
Kreisungspunktkurve k, in die Polbahntangente t und den die Polbahnnormale in P beriihrenden 
Kreis k,,* vom Durchmesser PM = m, wobeik, nicht entartet. Es ist dann] = cound |, = — D, 
der Kreis vom Durchmesser /, ist identisch mit dem Riickkehrkreis, kurz gesagt, es liegt momentan 
die umgekehrte Kardanbewegung vor, weswegen hier diese neue Bezeichnung ,,umgekehrte 
Kardanlage“ gewahlt wurde. Beim Kurbeltrieb erhalten wir diesen Fall in den Steglagen, (Abb. 6a/d) 
oder, wenn die Verbindungsgerade der Punkte P = 20 und H = 31 senkrecht zur Koppel steht 
(Abb. 7) oder beides zutrifft (Abb. 8a/f). 

y) Stationarer Wendekreis: Hat der Wendekreisdurchmesser D einen stationaren Wert, 
d. h. ist 0’ = 0, so geht m nach unendlich und entartet k,, in die Polbahnnormale n und den Kreis 
k,,** vom Durchmesser PL, = I, ebenso k, in die Polbahnnormale n und den Kreis k,* * vom Durch- 
messer PL, = |), wobei diese Kreise die Polbahntangente in P beriihren. Dieser Fall tritt beim 
Gelenkviereck dann auf, wenn Koppel und Steg parallel? sind oder wenn die Kurbel 4, A zur Pol- 
bahnnormale wird. ; 

1 W. Meyer zur Capellen, Ing.-Arch. 18 (1949) S. 308; O. Bottema, Nederlandsche Ak. Wetenschaapen, 
Vol, I—H, Nr. 6 (1949). 


5 Uber eine ausfiihrliche Anwendung dieser Stellung vgl.: W. Meyer zur Capellen, Konstruktion 8 (1956) 
Ou Os 
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; 3. Die Extrema der Abtriebsgeschwindigkeit. Nach diesen vorbereitenden Bemerkungen kommen 
wir nun zum eigentlichen Problem: Die Kurbel .4,.4 der Getriebe gema®B Abb. 1 und 2 werde mit 
konstanter Winkelgeschwindigkeit @, = Ww ) angetrieben und bewirke dadurch in Abb. 1 eine 
schwingende, in Abb. 2 eine umlaufende, aber ungleichformige Bewegung im Abtriebsglied B, B =3. 
Gefragt ist nach denjenigen Stellungen, in denen die Geschwindigkeit vg des Punktes B oder auch 
die Winkelgeschwindigkeit «, = Ws) dieses Gliedes ein Extremum hat. 


a) Bedingungen. Die Punkte 4 und B der allgemein bewegten Koppelebene beschreiben 
Kreise bzw. Kreisbégen, d. h. durchlaufen Bahnstellen mit (mindestens) vierpunktig beriihrendem 
Kriimmungskreis, gehéren, also der Kreisungspunktkurve k, an. Wenn weiter v, ein Extremum 
haben soll, so muf dv,/dt, d. h. die Tangentialbeschleunigung b, = dvz/dt gleich Null sein. Da A 
von vornherein nur Normal-, also keine Tangentialbeschleunigung hat und B in dem betrachteten 
Augenblick ebenfalls keine Tangentialbeschleunigung hat, miissen beide Punkte auf dem Gleichen- 


8 


Abb, 7. ,,Umgekehrte Kardanlage‘‘, PH senkrecht AB. 


a) Kurbelschwinge; 
b) Doppelkurbel, 


b 


kreis K, liegen. In beliebiger Lage beriihrt K, die Kurve k, in P langs der gemeinsamen Tan- 
gente n und hat dort mit k, zwei unendlich benachbarte Punkte gemeinsam, trifft also k,, noch in 
einem weiteren Punkt. Wenn aber k, und K, aufBer in P noch zwei Punkte, nimlich A und B 
gemeinsam haben sollen, so miissen K, und k,, ganz iibereinstimmen, d.h. aber die Kreisungs- 
punktkurve muB hier in den oben genannten Kreis kj vom Durchmesser m und die Polbahn- 
tangente t entarten, wobei ky und K, identisch sind. 

Ergebnis: Die Abtriebsgeschwindigkeit vg bzw. die Abtriebswinkelgeschwindigkeit w, = Wzo 
beim Kurbeltrieb hat dann ein Extremum, wenn sich das Getriebe in der umgekehrten Kardanlage 
befindet, d. h. wenn die Gerade PH senkrecht zur Koppel steht (Abb. 7), also der Abstand & des 
Poles P von der Senkrechten zu AB in H gleich Null wird (vgl. u.), hierzu kann in Sonderfallen 
auch die Steglage gerechnet werden. 

Die Steglage erfiillt diese Bedingung nur dann, wenn der Punkt G mit dem Punkt B iiberein- 
stimmt, d. h. wenn B auf der Polbahntangente t liegt (Abb. 8a/d). Dieser Sonderfall kann aber 
nur dann auftreten, wenn gewisse MaBbeziehungen zwischen den Gliedlangen erfiillt werden, d.h. 
wenn 

c2 + (d—a)?— 8? =0 (Abb. 8a, c) (6a) 


oder 


2 + (d+ a)?—b? =0 (Abb. 8b, d) (6b) 
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ist. Hierbei beziehen sich Abb. 8c und 8d auf die Doppelkurbel, wenn die Antriebskurbel a auf dem 
Steg liegt. Es kann aber auch die Abtriebskurbel auf dem Steg liegen (Abb. 8e und 8f), voraus- 
gesetzt, dah 
ce? + (b—d)?-—- a? =0 (Abb. 8e) (6c) 
oder 
co - (bd)? a7 = 0 (Abb. 8f) (6d) 
ist. 
b) Praktische Anwendung. Ist ein Getriebe mit bekannten Abmessungen gegeben, so 
ermittelt man die Geschwindigkeit des Abtriebes, d.h. vg, in bekannter Weise mit Hilfe der ge- 
drehten Geschwindigkeiten, wobei vy “a cm gezeichnet wird (w = ,,1“) (Abb. 9a und 10a) oder 


Abb, 8. Kurbeltriebe, bei denen in der Steglage die Koppel senkrecht zum Steg steht bzw. senkrecht zu PH. 
a) und b) Kurbelschwingen; c), d), e) und £) Doppelkurbeln. 
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mit Hilfe des Satzes von Euler, wonach 0g = 04+ Vp 4 ist und bg , senkrecht zur Koppel AB steht, 
beilaufig auch den Betrag cw = c Wy hat mit W = Wy als Winkelgeschwindigkeit der Koppel- 
ebene gegentiber dem Gestell (Abb. 9b und 10b). Tragt man vz tiber dem Kurbelwinkel a und 
damit uber der Zeit auf mit ~ = 0 in der rechten Steglage, da « ja proportional der Zeit t ist, und 
schreibt an v, den Wert eins an, so stellt die erhaltene Kurve mafstablich das Ubersetzungsver- 


Abb. 9, 10. Geschwindigkeitsermittlung. 
a) Methode der gedrehten Geschwindigkeiten; 
b) Satz von Euler; 

c) Hodograph. 

9: Kurbelschwinge; 10: Doppelkurbel. 


haltnis vp/v, dar (Abb. 11 und 12) baw. auch das Ubersetzungsverhaltnis w,/w, = a/b-(vg/vy)- 
Beim Zeichnen des Hodographen der gedrehten Geschwindigkeiten (Abb. 9c, 10c) erkennt man 
bereits etwa die Stellung, in der vz einen Héchst- oder Tiefstwert hat. Man zeichnet dann in einigen 
benachbarten Stellungen die Pole P und H, errichtet in H auf AB die Senkrechte und ermittelt 
den Abstand PF = & des Punktes P von dieser Senkrechten (Abb. 9a und 10a) oder fallt von P 
das Lot auf die Koppel, FuBpunkt X, so daB auch HX = & wird. Dann tragt man diese Strecke & 
iiber w auf (Abb. 11 und 12) und verbindet die so erhaltenen Punkte. Jetzt liefert der Schnittpunkt 
der &-Kurve mit der a-Achse denjenigen Wert a, in welchem v4 ein Extremum hat. 
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Diese Konstruktion la8t sich rasch und einfach durchfiihren.! Sie ist vor allem einfacher als die 
Ermittlung der Nullstelle von b,, indem man fiir einige benachbarte Stellungen die Beschleunigung des 
Punktes B ermittelt und die Tangentialkomponente b, iiber « auftragt. Abb. 13a zeigt die allgemeine 


& J 


Abb. 12. Doppelkurbel. 


a’ 


Abb. 11. Kurbelschwinge, 


b) Winkelgeschwindigkeit w = ® der Koppel bwz. w/w 


Abb, 11,12. Verlauf der Geschwindigkeiten uber dem Kurbelwinkel, 
a) Geschwindigkeit vp des Punktes B, bzw. vB/e43 


SI3 


SS y Sg 


g % 3 


Beschleunigungskonstruktion an der Kurbelschwinge, Abb. ihveine Stellung, in welcher b= 0 ists 
wahrend Bild 14a,b das Gleiche fiir die Doppelkurbel darstellt [vgl. Abs. d2)]. Wenn beilaufig die 
Antriebswinkelgeschwindigkeit ©,» also v, nicht konstant ist, so bleibt das Ubersetzungsverhailtnis 


* Der Vollstandigkeit halber ist die gesamte &-Kurve gezeichnet. 
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poate (,/w,, das ja nur geometrisch bedingt ist, erhalten, und die vorhin bestimmte Stellung 
ur ein Extremum von vp stellt dann ein Extremum eben dieses Ubersetzungsverhaltnisses dar. 


Abb. 13. Kurbelschwinge. 


Abb, 14. Doppelkurbel, 


Abb. 13, 14, Konstruktion der Beschleunigung bp des Punktes B vermige 
bp =b4 + bapa t Ba = np t Op. 
a) in beliebiger Lage; 
b) in der Stellung eines Extremums von vp, dh. bsp = by = 0; 


awe 


c) in der Stellung eines Extremums von v p 4 bzw. w.,d.h. bs By = 0- 
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c) Allgemeines Bild. Gibt man sich Wendekreis und Gleichenkreis vor (Abb. 15), also 
neben D noch den Quotienten w2/e, [Abs. 2a)], so erfiillen alle Gelenkvierecke, deren Koppel- 
punkte 4 und B auf dem Gleichenkreis bzw. der damit identischen Kreisungspunktkurve k* liegen, 
die betrachtete Bedingung, wobei A, und By nach der Euler-Savaryschen Formel (3), rechnerisch 
oder zeichnerisch zu ermitteln sind (Kontrolle PH | AB). Gibt man sich die Koppellange 


AB =c <g vor, so erhalt man unendlich viele Gelenkvierecke mit der betrachteten Bedingung; 
da man aber noch ce variieren 


kann, schlieBlich aber auch noch 
g, so erhalt man insgesamt oo? 
mogliche Gelenkvierecke, wobei 
dann im Einzelfall untersucht 
werden muB, ob das Getriebe dreh- 
fahig ist, also die Grashofsche Be- 
dingungen erfiillt, und ob es sich 
um eine Kurbelschwinge oder eine 
Doppelkurbel handelt. Da ein Ge- 
lenkviereck durch die vier GréBen 
a, b, c,d bestimmt ist, infolge der 
geometrischen Ahnlichkeit aber 
nur drei Parameter, z. B. a/c, b/c, 
d/c wesentlich sind, so fiihrt die 
rechnerische Ermittlung der Stel- 
lungen, in denen ein Extemum 
Abb, 15. Kurbeltrieb bei gegebenem Wendekreis und Gleichenkreis = Kreisungspunkt- won rs auftritt, auf Abhangigkeiten 
kurve k,,*, A und B beliebig auf k,,* angenommen (vgl. Text). von eben diesen drei variablen 
Parametern, wonach wieder eine 

dreifach unendliche Mannigfaltigkeit der Méglichkeiten gegeben ist. Die rechnerische Lésung ware 
allerdings sehr weitlaufig. Denn bereits beim exzentrischen Schubkurbelbetrieb! mit zwei Para- 
metern erhalt man fiir sina, worin « der Kurbelwinkel, eine Gleichung sechsten Grades, wahrend 
beim zentrischen Schubkurbeltrieb? mit einem Parameter nur eine Gleichung dritten Grades auftritt. 


Abb. 16. Schubkurbeltrieb in einer Stellung extremaler Geschwindigkeit des Punktes B. 


d) Sonderfalle. 1. Die Schubkurbel. Bei der soeben erwahnten geschrankten oder auch zen- 
trischen Schubkurbel als Sonderfall der Kurbelschwinge entartet die Schwinge 3 zu dem Schub- 
glied 3, d. h. es ist Ws) = 0, und es fallt der Gelenkpunkt B, da gerade gefiihrt, mit dem Beschleu- 
nigungspol J selbst zusammen, (Abb. 16). Da PH senkrecht AB stehen muB, kann nicht gleich- 
zeitig die Koppel AB auch senkrecht zum Kurbelradius stehen — womit die auch heute noch im 
Schrifttum zu findende Behauptung, vp hatte ein Extremum, wenn die Koppel Tangente an den 
Kurbelkreis ist, sofort widerlegt ist. 


1 Vgl. FuBnote 2 von S. 140. 2 Vel. FuBnote 1 von S. 140. 
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Eine Ausnahme hiervon macht nur die Steglage: Da B, im Unendlichen liegt, ist diese hier die 
Stellung, in der auch die Kurbel 4) A senkrecht zur Fithrung von B steht. Da eine Schubkurbel 
aus der Kurbelschwinge fiir b > co und d—> ov, aber d — b = e = endlich, entsteht, folgt aus (6a), 
wenn dort 6 = d—e eingesetzt wird, (c? + a? — e?)/d — 2(a — e) = 0 und fiir d—> 0 die Bedin- 
gung a =e, d.h. die ge- 
rade Bahn vonB muB den 
Kurbelkreis _ beriihren. 
Die Koppel kann aber 
dann den Kurbelkreis 
nur beriihren, wenn 

a= 3 7/2 

ist. Beilaufig liegt dann 
P unendlich fern [vel. 
folgenden Abs. 2)] und 
ist )p = 0D, Wenn also 
a =e oder A=—€e mit 
A =al/c und € = e/c ist, 
so tritt das eine Extre- 
mum in dieser Stellung 
auf 1, 

AufdieVerzweigungs- 
lagen, dieja letztenEndes 
auch Steglagen sind, wird 
unten noch kurz einge- 
gangen. 

2. Die Parallel-Lagen. 
Da fir «o = 3 7/2 bei der 
Schubkurbel der Momen- 
tanpol ins Unendliche 
riickt, sind die urspriing- 
lichen Voraussetzungen 
eigentlich nicht mehr er- 
fiillt, und es bedarf da- 


her einer besonderen Be- Abb. 17. Parallel-Lagen der Kurbelschwinge mit Beschleunigungskonstruktion. 


a) Vierecklage und b) zugehérige Beschleunigungskonstruktion; 
c) Uberkreuzlage und d) zugehérige Beschleunigungskonstruktion. 


Ro 


trachtung: Bei der Kur- 
belschwinge liegt der 
Momentanpol im Unendlichen, wenn Kurbel und Schwinge parallel sind (Abb. 17a—d); dieser Fall 
kann also bei einer nicht durchschlagenden Doppelkurbel nicht auftreten. Ohne nun auf die Ent- 
artung der Ortskurven einzugehen, sei hier unmittelbar von der Beschleunigung ausgegangen: 
Da vz = vy ist, wird vg, = 0 und damit m = Wy) = 0. Im Allgemeinen gilt fir die Beschleuni- 


Foo 


Abb. 18. Parallel-Lagen zweier Kurbelschwingen, wenn Koppel senkrecht Kurbel und damit vp ein Extremum. 


gung des Punktes B (Abb. 13 und 14) bg = by +6, 24+ 5,84, worin b, p4 = vBy/e = ¢ w” ist, 
die Richtung BA hat, und hier gleich Null wird, und worin b,,4 senkrecht zu AB verlauft. 
Andererseits gilt 6, = 6, , + b,g, worin b,, = v3/b, also hier gleich v4/b ist, und die Tangential- 


1 Vel. FuBnote 2 von S. 140. 
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komponente b,,, oben kurz mit b, bezeichnet, senkrecht zu By, B verlauft, wobei aber b, , parallel 
b, 4 = 64 ist. Nach Abb. 17b, d ist 34 = b,, und diese Strecke wird gleich Null, wenn in der 
Parallel-Lage zudem noch die Koppel senkrecht zur Kurbel baw. zur Schwinge steht (Abb. 18a, b). 
Das stimmt auch mit der friiheren Bedingung ,,PH senkrecht Koppel“ iiberein. Die Stellungen 
nach Abb. 18a,b sind aber nur bei bestimmten 

MaBbeziehungen méglich, d. h. wenn 


ce + (b—a)?—d?=0 (Abb. 18a) (7a) 


oder 


2 +(b +a)?—d2=0 (Abb. 18b) (7b) 


ist. Setzt man beilaufig in die letzte Gleichung 
d—b=e ein und 1]a8t dannd nach unendlich 
gehen, so folgt wieder a =e bei der Schubkurbel 
und die dort angegebene Kurbelstellung. 


4, Das Extremum der Winkelgeschwindigkeit der 
Koppelebene. a) Allgemeine Bedingung. Die 
Winkelgeschwindigkeit jw[= W.)| der Koppelebene 


gegeniiber dem Steg hat dann ein Extremum, wenn 


b 


Abb, 19. Stellungen fiir extremale 
Winkelgeschwindigkeit w,, der Koppel- 
ebene, PH senkrecht B,P. 


a) Kurbelschwinge, 
b) Doppelkurbel. 


ihre Ableitung, d.h. die Winkelbeschleunigung ¢ der Koppelebene durch Null geht. Da der 
Durchmesser des Gleichenkreises durch g = D w?/e gegeben war, muB g—> co gehen, d.h. der 
Gleichenkreis entartet in die Polbahnnormale n. Da 
nun aber die Punkte A und B der Koppelebene gleich- 
zeitig auf der Kreisungspunktkurve k, liegen miissen, 
A aber auf der Polbahnnormale als dem entarteten 
Gleichenkreis, mu auch die Kurve k, entarten, 
und zwar in die Polbahnnormale n und in einen, 
die Polbahntangente t in P berithrenden Kreis k** 
vom Durchmesser PL =1, auf dem auch der Punkt B 
liegen muB. Es handelt sich also um den oben unter 
Abs. 2c) y hervorgehobenen Sonderfall. Nun ist nach 
der Bobillierschen Konstruktion (Abb. 19, a, b) die 
Antriebskurbel A,A dann auf der Polbahnnormale 
gelegen, wenn HP senkrecht zu PB steht: 

Ergebnis: Die Winkelgeschwindigkeit jm der 
Koppelebene hat dann ein Extremum, wenn die 
Kurbel bzw. die Antriebskurbel bei der Doppel- 
kurbel auf der Polbahnnormale liegt, d.h. wenn 
die Verbindungsgerade zwischen Pol P=20 und 
Relativpol H=31 senkrecht zur Schwinge bzw. zur 
Abtriebskurbel 3 steht. 

Denkt man sich einmal A, A festgehalten und betrachtet die Bewegung des Gliedes AB gegen- 
iiber diesem, wobei das Glied B, Ay gleichférmig angetrieben sei (d und a bilden ja den Winkel a), 


23 


Abb. 20. Relativpole gemaB Abb. 19b (Doppelkurbel). 
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so ist B,B die Koppel, und diese befindet sich in der oben erérterten umgekehrten Kardanlage, 
da die Verbindungsgerade vom Momentanpol 13 dieser Koppelebene und dem Hilfspol 20 seukrecht 
zur ,,Koppel B,B steht. Dies folgt aber auch aus den allgemeinen Satzen iiber die Drehung von 
drei Ebenen gegeneinander: In Abb. 20 sind alle Relativpole eingetragen, und fiir die Winkelge- 
schwindigkeiten (49, W29, W2, der drei Ebenen 1, 2,0 gegeneinander um die Achsen (Pole) 10, 20 
21 gilt doch 


Woy = Wiy + Wo, Oder Woy = Woy — Wy, 


Abb. 22. Kurbeltriebe fiir extremale Winkelgeschwindigkeit ,, 
der Koppelebene bei gegebenem Wendekreis und gegebener ent- 
arteter Kreisungspunktkurve k,,**. Es stellt k,** die zugehérige 
Angelpunktkurve dar. 
a) L, L, in der oberen, b) L, L, in der unteren Halbebene. L in 
der oberen Halbebene auBerhalb Ky liefert L, in der unteren 
Halbebene. 


L 


wobei Wi) = M, = constant, Wy = W und w,, die Winkelgeschwindigkeit des Gliedes 2 gegen | 
bedeutet. Somit hat auch w,, ein Extremum, wenn W,) ein Extremum hat, und W., hat ein Ex- 
tremum, wenn eben 13—20 senkrecht 30—20 ist. 

b) Koppel parallel Steg. Die stationdre Lage ist auch in dieser Stellung verwirklicht, aber 
sie liefert nur dann ein Extremum in w, wenn gleichzeitig AA) Polbahnnormale ist und PH senk- 


ll 


152 Meyer zur Capellen: Die Extrema der Geschwindigkeiten an Kurbeltrieben Ingenieur-Archiv 


recht BB, steht, wobei H unendlich fern, also PH || AB und || Ay By ist (Abb. 21a, b). Eine solche 
Lage kann aber nie bei der Kurbelschwinge, wohl aber bei der Doppelkurbel auftreten. Denkt man 
sich wieder 4A, festgehalten, so liegt die oben betrachtete Parallel-Lage vor [Abs.3d) 2)]. Bei der 


Abb. 23. Schubkurbeltrieb in der Stellung extremaler Winkelgeschwindigkeit w,, der Koppelebene: Winkel BPH = PHB, = PHA, = 2/2. 


& Doppelkurbel ist die genannte Stellung nur bei be- 
stimmten MaBbeziehungen erfillt: 
s b+ (c—d)2}—a2?=0, Abb.2la, (8a) 
oder 
Ky b? + (c + d)?— a*—0 Abb. 21b. (8b) 


c) Praktische Anwendung. Die Winkel- 
geschwindigkeit @ der Koppelebene kann zeich- 


| nerisch aus v, = PA w oder vg = PBw oder noch 


besser aus vp 4 = ABw = cw (vgl. Abs. 3b) leicht 

Boo gewonnen werden: Tragt man vz 4 iber dem Kurbel- 

winkel « auf, so stellt die erhaltene Kurve bereits 

Boo maBstablich den Verlauf von w dar. In der Nahe 

der Extrema fallt man nun 4hnlich wie oben in 

Abs. 3 b) vom Relativpol H = 13 das Lot auf den 

Strahl B)BP, FuBpunkt Y, Abb. 9a, und tragt die 

Strecke PY =7 tiber « auf. Der Schnittpunkt der 

Kurve 7(«) mit der «-Achse bestimmt dann die ge- 

suchten Stellungen der Extrema (Abb. 11b und 12b). 

Auch dieser Weg ist einfacher als der der Beschleuni- 

gungskonstruktion, bei welcher man b, p , = c € tiber 
a auftragen miiBte: b, 5 4 = 0 bedeutet ¢ = 0.1 

d) Allgemeines Bild. Gibt man sich den 

Ky Wendekreis K,, und den Kreis ki* (die entartete 

Kreisungspunktkurve k,) vom Durchmesser PL =I 

vor (Abb. 22a, b), so kann A auf n und B auf k** 

gewahlt werden. A, und B, folgen aus der Euler- 

Savaryschen Formel, wobei B, auf dem Kreis k** 

Abb. 24. Kurbelschleife (B — oo) in einer Stellung extremaler vom Durchmesser ly hegt [vgl. Gl. Sb]. In ahnlicher 

Winkelgeschwindigkeit 9 der Koppel = wy des Abtriebs- Weise wie oben ergibt sich dann wieder eine drei- 

Ce er aeeenee s fach unendliche Mannigfaltigkeit. 

e) Sonderfalle. «) Schubkurbel. Hier fallt der geradgefiihrte Punkt B auf den Wendekreis 

selbst, und da B gleichzeitig dem Kreis ki* angehért, sind Wendekreis und Kreisungspunktkurve 

ki* identisch: Alle Punkte des Wendekreises beschreiben Undulationspunkte! Da 1, co geht, 


4 Vell Abb: 13¢ und 14. 
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ist auch hier t’ = 0, d.h. Ry = D und R = D/2: Es liegt gleichzeitig die Kardanlage vor. Dies 
zeigt sich auch in Abb. 23, wonach nicht nur HP senkrecht PBB, (B)— ©) sondern auch PH 
senkrecht zum Steg A,B, steht. Wie an anderer Stelle gezeigt wurde!, ergibt sich fiir sina eine 
quadratische Gleichung und daraus eine einfache Bestimmungsgleichung fiir den Kurbelwinkel 
a =a**. Dann liegt ein Extremum in w fiir 1 = 0** und fir « = 180° — q** vor. 


B) Kurbelschleife. Die Kurbelschleife wurde in Ziff. 3 noch nicht erwahnt, da B-> oo geht. 
Dann entartet der Kreis k** in die unendlich ferne Gerade, und aus dem entsprechenden Teil der 
Angelpunktkurve k** wird der Riickkehrkreis K,, vom Durchmesser |, = (—) D (Abb. 24)2. Es 


Abb. 25. Verzweigungslage eines durchschlagenden Getriebes, MaBe. 


kann also A, beliebig auf n gewahlt werden und B, auf K,; man erhalt im allgemeinen Fall die 
geschrankte Kurbelschleife?: Die in A mit der Kurbel 1 gelenkig verbundene Koppelebene 2 gleitet 
durch die um B, drehbare entartete Schwinge 3. Die Winkelgeschwindigkeit j.,) ist immer gleich 
Wo, und es steht hier sowohl HP senkrecht zu B)PB wie auch PH senkrecht zur Koppel AB! 


Wird der Punkt B, im Riickkehrpol ¥Y selbst gewahlt, so erhalt man die zentrische Kurbel- 
schleife mit e = AS = 0, und es bestatigt sich, daB das Abtriebsglied in der Steglage ein Maximum 
bzw. Minimum seiner Winkelgeschwindigkeit durchlauft. 4 


y) Totlage. Auch die Totlage einer Kurbelschwinge gehért zu der hier betrachteten Stellung 
dann, wenn Schwinge und Kurbel dabei aufeinander senkrecht stehen: Die Schwinge liegt auf der 
Polbahntangente und die Kurbel einschlieBlich der Koppel, wie erértert, auf der Polbahnnormale. 


gd 


Abb. 26. Wie Abb. 25, aber Konstruktion der beiden Geschwindigkeiten vB und Bo 
gemiB den beiden Momentanpolen P, und P, (vgl. Text). ts 


5. Verzweigungslagen. Wenn die Punkte A, Ay und B, By eines Kurbeltriebes in einer Richtung 
und zwar auf der Polbahnnormale liegen, so ist die Bewegung des Getriebes nicht mehr eindeutig. 
Es kann das Getriebe durchschlagen, es kann aber auch die Geschwindigkeit um einen endlichen 
Wert springen. Diese Stellung gehért als Grenzfall in die Stellung der extremalen Abtriebsgeschwin- 
digkeit wie auch der extremalen Koppelwinkelgeschwindigkeit — hat dariiber hinaus aber keine 
praktische Bedeutung. Beim Schubkurbeltrieb® ist dieser Grenzfall durch a + e = oderA+e=1 
gegeben [vgl. Abs. 3 d)], d. h. alle untersuchten Schubkurbeltriebe miissen die Bedingung A -+é <1 
erfiillen, wobei 2 +- € = 1 die Grenzkurven liefert. 

Es sei hier nur kurz erlautert, wie die Geschwindigkeiten in der Verzweigungslage ermittelt 
werden kénnen (Abb. 25). Der Momentanpol der Koppelebene ist zunachst nicht eindeutig angeb- 
bar, er sei im Abstande x vom Kurbelpunkt 4 angenommen. Dann gilt nach der Euler-Savaryschen 


1 Vel. FuBnote 2 von S. 140. 

2 Beilaufig wird auch R, = — D/2 und R = — D: umgekehrte Kardanlage. 

3 Uber eine praktische Anwendung vgl. W. Meyer zur Capellen, Werkstatt und Betrieb 89 (1956) S. 67. 
4 Vel. W. Meyer zur Capellen, Werkstatt und Betrieb 89 (1956) hy Oxy le 

5 Vel. FuBnote 2 von S. 140. 
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Formel (2) 
| 1 1 1 | 


x BG BSE és x+e+ b D 


oder geordnet 
x? (a—b)+2acx-+ac(ec+b)=0 


mit den Lésungen 


ac ac /od— 
4226+), = OCS) eee 


b—a 


wobei hier a + b =c + d. 


(9) 


(10) 


(11) 


Es gibt also zwei Momentanpole! und damit zwei Abtriebsgeschwindigkeiten, wie in Abb. 26 
angegeben, und auch zwei Wendekreise. Es ist nicht die Absicht hier im Einzelnen auf die ver- 
schiedenen méglichen Faille einzugehen, zu denen auch der Antiparallelkurbeltrieb oder das Gallo- 


waygetriebe gehéren, sondern es sollte nur kurz auf den Grenzfall hingewiesen werden. 


Beim Schubkurbeltrieb? ergeben sich mit b—> co und d -> oo die einfachen Werte x, = Vac 


Tee ee sore 
und x4. jac. 


(Eingegangen am 13. Juli 1956) 


Anschrift des Verfassers: Dr.-Ing. W. Meyer zur Capellen, Aachen, Ronheiderweg 46 


* Die beiden Pole sind die Doppelpunkte der durch die Paare A, By, B, Ay bestimmten Inyolution, Vel. 


R. Miller (FuBnote 1 von S, 141), 
2 Vel. FuBnote 2 von S. 140. 


i aa Wahrend der Korrekturen wurde mir die Arbeit von F. Freudenstein »On the Maximum and Minimum 
Velocities and the Accelerations in Four-Link Mechanisms‘ Trans. of ASME (1956) S. 779/787 -ASME Paper 
55-SA-20 bekannt, in welcher auf anderem Wege die gleichen Bedingungen wie in Abs. 3a iiber das Extremum 


der Abtriebsgeschwindigkeit gefunden wurden. 


q 


Die selbsttiitige Regelung. Theoretische Grundlagen 


mit praktischen Beispielen,. Von Professor Dr.-Ing. A. Leonbard, Stuttgart. Zweite neubear- 


beitete Auflage. Mit 319 Abbildungen~ XII, 376 Seiten Gr.-8°. 1957. : Ganzleinen DM 39,— 


/ ‘ 

Die Regelungstechnik hat als verhaltnismaBig junges Teilgebiet der Technik in den letzten Jahren einen 
ungeheuren Aufschwung erlebt. An dieser Aufwartsentwicklung sind Techniker, Naturwissenschaftler 
und Mathematiker aller zivilisierten Lander beteiligt. In das nunmehr in 2. Auflage vorliegende Werk 
sind neben den bereits als klassisch anzusehenden Methoden miglichst weitgehend auch die in, der Zwischen- 


zeit international erarbeiteten neuesten Verfahren fiir die Behandlung von Regelproblemen eingearbeitet 


worden. 
( 


Seit Erscheinen des Buches ,,Die selbsttatige Regelung in der Elektrotechnik“ im Jahre 1940 im gleichen 
Verlag und vom gleichen Verfasser, das als eines der ersten auch neuere Verfahren der Regelungstechnik 
behandelt hat und als Vorlaufer des jetzigen Werkes biedechaa ist, sind zwar Daliveishie Junaminmaessaente 
Werke in verschiedenen Sprachen uber das Gebiet der Regelung eatstandens teotedéns wird das vorliegende 
Werk dank seiner Eigenart, theoretische Methoden immer nur im Zusammenhang mit durchgerechneten 


Beispielen zu behandeln, seinen Platz in der Regler-Literatur behaupten. - 


Auch in der 2. Auflage, bei der die Benennungen und Begriffe weitgehend dem inzwischen erschienenen 
Normblatt angeglichen sind, Neldon zunichst in einem ersten Abschnitt die Grundlagen der Regelungs- 
technik nnd anschlieSend die Einzelglieder des Regelkreises, also im wesentlichen Regler und Regelstrecke, 
hekasdelt. Ein zweiter Abschnitt beschaftigt sich mit den verschiedenen vachuseeicker und graphischen 


Methoden zur Ermittlung des Regelvorganges, ein dritter mit dem Stabilitatsproblem. SchlieSlich wird 


_im vierten Teil ausfiihrlich auf die Synthese des Regelkreises, also vor allem auf die zweckmaBige Bestim- 


o 


mung frei wahlbarer Reglerkonstanten eingegangen. 


Inhaltsiibersicht: Grandingcets Allgemeines tiber Regelung. — Verhalten von MeBwerken baw. Reg- 
lern ‘ahihé Hilfsenergie. — Verhalten von Einzelgliedern (Regelgliedern) des Regelkreises. — Versuchs- 
etiniichs Feststellung des Verhaltens von Einzelgliedern des Regelkreises. — Ermitilung des Regelvor- 
ganges: Klassische Verfahren zur Ermittlung des Regelvorganges mit Hilfe der Differentialgleichung. — 
Die Ermittlung des Regelvorganges mit Hilfe der Laplace-Transformation. — Ermittlung der charakteristi- 
schen Gleichung des Regelvorganges mit Hilfe des Frequenzganges. — Ermittlung des Regelvorganges mit 
Hilfe des Frequenzganges. — Bauitiang des Regelvorganges bei. nichtlinearen Regelgliedern. — Die 
Stabilitat der ie le Stabilitatsuntersuchungen. — Verbesserung der Stabilitat. — Festlegung frei 
wahlbarer Regelkonstauten. re Versctiedone Methoden fiir die aweckmaBige Wahl der Regelkonstanten. — 
Kurvenblatter fiir die zweckmaBige Bestimmung frei wahlbarer Regelkonstanten, — Anhang: Die Lésung 
von Chidwnden ‘ohana Grades mit konjugiert komplexen Wurzeln, — Literaturverzeichnis, — Sach- 


yerzeichnis. 


SPRINGER-VERLAG > BERLIN / GOUTTINGEN/HEIDELBERG 


¢e 


Elektr essung mechanischer 
ektrische M o h her 
pds Von Dr.-Ing. Paul M. Pflier, Nirnberg. Vi erte neubearbeitete Auflage. Mit 
Gro i} en ® 349 Abbildungen. VIII, 276 Seiten Gr.-8°. 1956. Ganzleinen DM 33,— 
Inhaltsiibersicht; A. Grundlagen der elektrischen’ Messung. I. Vorziige elektrischer MeBver- 
fahren. Il. Die Eigenschaften der elektrischen MeSgerate.. Allgemeines. Grundbegriffe der 
MeBtechnik. VDE-Vorschriften. Symbole und Schaltzeichen. MeSfehbler, Empfindlichkeit. Spannungs- 
* einfluB. ‘Spannungsk onstanthalter, Spannungsunabhangige MeBgeraite. Frequenzeinflub. Temperatur- 
- einfluB. Schwingungseigenschaften des beweglichen Organs. Ferniibertragung. Til. MeBwerke und 
MeByverfahren aaa Durchfiihren von Rechenoperationen. Addition und Subtraktion mit dem 
DrehspulmeBwerk. Multiplikation mit dem elektrodynamischen MeBwerk. Multiplikation mit dem Hall- 
effekt. Division mit dem Kreuzspulinstrument. IV. Elektrische Abbildung mathemati cher B e- 
ziehungen. Addition zweier GréBen. Subtraktionsschaltungen. Addition und Subtraktion von Wechsel- 
spannungen und Strémen. Divisionsschaltungen. Multiplikationsschaltungen, Differentiation. Integration. 
Briickenschaltungen. — B. Die Méglichkeiten, mechanische GréBen elektrisch zu messen. I. E nergieum- 
former. Allgemeines. Induktionsprinzip.’ Piezoelektrizitat. Stroémung selektrizitat: Thermoelektrische 
Umformer. Lichtelektrische Umformer. Induktive Umformer, I. Physikalische Zusammenhiange 
zwischen nichtelektrischen und elektrischen GréBen. Allgbinethes. Beeinflussung des Widerstandes 
eines Leiters. Steuerung der Kapazitat eines Kondensators. Mechanische Steuerung der Induktivitat einer 
Spule. Zeit und Elektrizitatsmenge. Elastische Beanspruchung und Eigenfrequenz. III. Wellen aus- 
breitung und Bestrahlung. /Allgemeines. Radioaktive und Rontgenstrabler. StrahlungsmeBgerate. 
Ultraschall-Bestrahlung. Lichtdurchstrahlung (Photcelastisches Verfahren. IV. MechanischeSteuerung 
eines elektrischen Kreises, Allgemeines. Kontinuierlich arbeitende Verfahren. Impulsverfahren. cae 
_¢. MeBverfahren. I. Lan genmessun g. Stellungsanzeige und Wegmessung. Elektrische Mikrometer. 
Dickenmessung. II. Kraftmessung. Dehnungsmessung.  KraftmeBdosen. Elektrische Waagen. Band- 


zugmesser. Schnittkraftmesser. Gasdruckmesser. Drehmomentmesser. III. Geschwindigkeitsmes- 


sung. Drehzahlmesser. Drehzahldifferenz- und Schlupfmessung. Messung yon Fahrzeyggeschwindigkeiten. 


Strémungsmesser. IV. Messung von Beschleunigung, Schwingungen, Erschitterungen. All- 


gemeines. Drehbeschleunigungsmesser. Beschleunigungsmesser. V. Zeitmessung. Langzeitmesser. 


Zeitschreiber. Hlektromache Zeitmesser. Zeitwaagen. — Literaturverzeichnis. Namenverzeichnis. — 


Sachverzeichnis. 


Die nich aniachn GréBen an Ort und Stelle mit mechanischen oder optischen Mitteln zu messen ist auch 
heute noch in vieley#f allen iiblich und ausreichend. Uberall aber, wo die MeBstelle schlecht zuganglich ist, 
zwischen MeBstelle und Ableseort erafere Entfernungen liegen, mebrere GrofSen zueinander in Beziehung 
gebracht werden sollen oder wo die MeBgréBe irgendwelche Vorgiinge steuern soll, haben sich die elek- 


trischen MeBverfahren als zweckmaBig erwiesen und in standig wachsendem Umfang eingeschaltet. 


Die 4. Auflage dieses allgemein anerkannten Buches behalt im wesentlichen die friihere Einteilung bei, 


versuchi jedoch Grundsatzliches vom Speziellen scharfer zu trennen und die Systematik zu verbessern; 
sie befaBt sich auch in stirkerem MaBe mit den neuerdings in den Vordergrund geriickten MeBverfahren. 


“Das Buch behandelt zunachst die Eigenschaften der elektrischen MeBgerate und die Vorziige der elek- ‘ 


trischen Messung, anschlieBend werden die Méglichkeiten der mechanisch-elektrischen Umwandlung und 
die Beziehungen zwischen mechanischen und elektrischen GréSen gezeigt, im letzten Abschnitt werden 


die zweckmaBigen MeBverfahren fiir die einzelnen mechanischen GréBen besprochen. 
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